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Recommended Books and Resources 


• M. Peskin and D. Schroeder, An Introduction to Quantum Field Theory 

This is a very clear and comprehensive book, covering everything in this course at the 
right level. It will also cover everything in the “Advanced Quantum Field Theory” 
course, much of the “Standard Model” course, and will serve you well if you go on to 
do research. To a large extent, our course will follow the first section of this book. 


There is a vast array of further Quantum Field Theory texts, many of them with 
redeeming features. Here I mention a few very different ones. 

• S. Weinberg, The Quantum Theory of Fields, Vol 1 

This is the first in a three volume series by one of the masters of quantum field theory. 
It takes a unique route to through the subject, focussing initially on particles rather 
than fields. The second volume covers material lectured in “AQFT”. 

• L. Ryder, Quantum Field Theory 

This elementary text has a nice discussion of much of the material in this course. 

• A. Zee, Quantum Field Theory in a Nutshell 

This is charming book, where emphasis is placed on physical understanding and the 
author isn’t afraid to hide the ugly truth when necessary. It contains many gems. 

• M Srednicki, Quantum l Field Theory 

A very clear and well written introduction to the subject. Both this book and Zee’s 
focus on the path integral approach, rather than canonical quantization that we develop 
in this course. 


There are also resources available on the web. Some particularly good ones arc listed 
on the course webpage: http://www.damtp.cam.ac.uk/user/tong/qft.html 
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0. Introduction 


“There are no real one-particle systems in nature, not even few-particle 
systems. The existence of virtual pairs and of pair fluctuations shows that 
the days of fixed particle numbers are over.” 

Viki Weisskopf 

The concept of wave-particle duality tells us that the properties of electrons and 
photons are fundamentally very similar. Despite obvious differences in their mass and 
charge, under the right circumstances both suffer wave-like diffraction and both can 
pack a particle-like punch. 

Yet the appearance of these objects in classical physics is very different. Electrons 
and other matter particles are postulated to be elementary constituents of Nature. In 
contrast, light is a derived concept: it arises as a ripple of the electromagnetic field. If 
photons and particles are trucly to be placed on equal footing, how should we reconcile 
this difference in the quantum world? Should we view the particle as fundamental, 
with the electromagnetic held arising only in some classical limit from a collection of 
quantum photons? Or should we instead view the held as fundamental, with the photon 
appearing only when we correctly treat the held in a manner consistent with quantum 
theory? And, if this latter view is correct, should we also introduce an “electron held”, 
whose ripples give rise to particles with mass and charge? But why then didn’t Faraday, 
Maxwell and other classical physicists hnd it useful to introduce the concept of matter 
helds, analogous to the electromagnetic held? 

The purpose of this course is to answer these questions. We shall see that the second 
viewpoint above is the most useful: the held is primary and particles are derived 
concepts, appearing only after quantization. We will show how photons arise from the 
quantization of the electromagnetic held and how massive, charged particles such as 
electrons arise from the quantization of matter helds. We will learn that in order to 
describe the fundamental laws of Nature, we must not only introduce electron helds, 
but also quark helds, neutrino helds, gluon helds, W and Z-boson helds, Higgs helds 
and a whole slew of others. There is a held associated to each type of fundamental 
particle that appears in Nature. 

Why Quantum Field Theory? 

In classical physics, the primary reason for introducing the concept of the held is to 
construct laws of Nature that are local. The old laws of Coulomb and Newton involve 
“action at a distance”. This means that the force felt by an electron (or planet) changes 
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immediately if a distant proton (or star) moves. This situation is philosophically un¬ 
satisfactory. More importantly, it is also experimentally wrong. The held theories of 
Maxwell and Einstein remedy the situation, with all interactions mediated in a local 
fashion by the held. 

The requirement of locality remains a strong motivation for studying held theories 
in the quantum world. However, there are further reasons for treating the quantum 
held as fundamental 1 . Here I'll give two answers to the question: Why quantum held 
theory? 

Answer 1: Because the combination of quantum mechanics and special relativity 
implies that particle number is not conserved. 

Particles are not indestructible objects, made at the 
beginning of the universe and here for good. They can be 
created and destroyed. They are, in fact, mostly ephemeral 
and fleeting. This experimentally verihed fact was hrst 
predicted by Dirac who understood how relativity implies 
the necessity of anti-particles. An extreme demonstra¬ 
tion of particle creation is shown in the picture, which 
comes from the Relativistic Heavy Ion Collider (RHIC) at 
Brookhaven, Long Island. This machine crashes gold nu¬ 
clei together, each containing 197 nucleons. The resulting 
explosion contains up to 10,000 particles, captured here in 
all their beauty by the STAR detector. 

We will review Dirac’s argument for anti-particles later in this course, together with 
the better understanding that we get from viewing particles in the framework of quan¬ 
tum held theory. For now, we’ll quickly sketch the circumstances in which we expect 
the number of particles to change. Consider a particle of mass m trapped in a box 
of size L. Heisenberg tells us that the uncertainty in the momentum is A p > h/L. 
In a relativistic setting, momentum and energy are on an equivalent footing, so we 
should also have an uncertainty in the energy of order A E > hc/L. However, when 
the uncertainty in the energy exceeds A E = 2me 2 , then we cross the barrier to pop 
particle anti-particle pairs out of the vacuum. We learn that particle-anti-particle pairs 
are expected to be important when a particle of mass m is localized within a distance 
of order 



me 


A concise review of the underlying principles and major successes of quantum field theory can be 
found in the article by Frank Wilczek, http://arxiv.org/abs/hep-th/9803075 
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At distances shorter than this, there is a high probability that we will detect particlc- 
anti-particle pairs swarming around the original particle that we put in. The distance A 
is called the Compton wavelength. It is always smaller than the de Broglie wavelength 
A dB = h/\p\. If you like, the de Broglie wavelength is the distance at which the wavelike 
nature of particles becomes apparent; the Compton wavelength is the distance at which 
the concept of a single pointlike particle breaks down completely. 

The presence of a multitude of particles and antiparticles at short distances tells us 
that any attempt to write down a relativistic version of the one-particle Schrodinger 
equation (or, indeed, an equation for any fixed number of particles) is doomed to failure. 
There is no mechanism in standard non-relativistic quantum mechanics to deal with 
changes in the particle number. Indeed, any attempt to naively construct a relativistic 
version of the one-particle Schrodinger equation meets with serious problems. (Negative 
probabilities, infinite towers of negative energy states, or a breakdown in causality are 
the common issues that arise). In each case, this failure is telling us that once we 
enter the relativistic regime we need a new formalism in order to treat states with an 
unspecified number of particles. This formalism is quantum held theory (QFT). 

Answer 2: Because all particles of the same type are the same 

This sound rather dumb. But it’s not! What I mean by this is that two electrons 
are identical in every way, regardless of where they came from and what they’ve been 
through. The same is true of every other fundamental particle. Let me illustrate this 
through a rather prosaic story. Suppose we capture a proton from a cosmic ray which 
we identify as coming from a supernova lying 8 billion lightyears away. We compare 
this proton with one freshly minted in a particle accelerator here on Earth. And the 
two are exactly the same! How is this possible? Why aren’t there errors in proton 
production? How can two objects, manufactured so far apart in space and time, be 
identical in all respects? One explanation that might be offered is that there’s a sea 
of proton “stuff” filling the universe and when we make a proton we somehow dip our 
hand into this stuff and from it mould a proton. Then it’s not surprising that protons 
produced in different parts of the universe are identical: they’re made of the same stuff. 
It turns out that this is roughly what happens. The “stuff” is the proton held or, if 
you look closely enough, the quark held. 

In fact, there’s more to this tale. Being the “same” in the quantum world is not 
like being the “same” in the classical world: quantum particles that are the same are 
truely indistinguishable. Swapping two particles around leaves the state completely 
unchanged — apart from a possible minus sign. This minus sign determines the statis¬ 
tics of the particle. In quantum mechanics you have to put these statistics in by hand 
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and, to agree with experiment, should choose Bose statistics (no minus sign) for integer 
spin particles, and Fermi statistics (yes minus sign) for half-integer spin particles. In 
quantum held theory, this relationship between spin and statistics is not something 
that you have to put in by hand. Rather, it is a consequence of the framework. 

What is Quantum Field Theory? 

Having told you why QFT is necessary, I should really tell you what it is. The clue is in 
the name: it is the quantization of a classical held, the most familiar example of which 
is the electromagnetic held. In standard quantum mechanics, we’re taught to take the 
classical degrees of freedom and promote them to operators acting on a Hilbert space. 
The rules for quantizing a field are no different. Thus the basic degrees of freedom in 
quantum held theory are operator valued functions of space and time. This means that 
we are dealing with an infinite number of degrees of freedom — at least one for every 
point in space. This infinity will come back to bite on several occasions. 

It will turn out that the possible interactions in quantum held theory are governed 
by a few basic principles: locality, symmetry and renormalization group how (the 
decoupling of short distance phenomena from physics at larger scales). These ideas 
make QFT a very robust framework: given a set of helds there is very often an almost 
unique way to couple them together. 

What is Quantum Field Theory Good For? 

The answer is: almost everything. As I have stressed above, for any relativistic system 
it is a necessity. But it is also a very useful tool in non-relativistic systems with many 
particles. Quantum held theory has had a major impact in condensed matter, high- 
energy physics, cosmology, quantum gravity and pure mathematics. It is literally the 
language in which the laws of Nature are written. 

0.1 Units and Scales 

Nature presents us with three fundamental dimensionful constants; the speed of light c, 
Planck’s constant (divided by 27r) h and Newton’s constant G. They have dimensions 

[c] = LT - 1 

[h] = l 2 mt ~ 1 

[G] = L^M-'T” 2 

Throughout this course we will work with “natural” units, defined by 

c = h= 1 (0.1) 
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which allows us to express all dimensionful quantities in terms of a single scale which 
we choose to be mass or, equivalently, energy (since E = me 2 has become E = m). 
The usual choice of energy unit is eV , the electron volt or, more often GeV = 10 9 eU or 
TeV = 10 12 eU. To convert the unit of energy back to a unit of length or time, we need 
to insert the relevant powers of c and h. For example, the length scale A associated to 
a mass m is the Compton wavelength 



me 


With this conversion factor, the electron mass m e = 10 6 eV translates to a length scale 
A e ~ 10 _12 m. (The Compton wavelength is also defined with an extra factor of 27r: 
A = 2nh/mc.) 


Throughout this course we will refer to the dimension of a quantity, meaning the 
mass dimension. If A" has dimensions of (mass) d we will write [A"] = d. In particular, 
the surviving natural quantity G has dimensions [G] = —2 and defines a mass scale, 


G = 


he 

M2 


1 

M2 


( 0 . 2 ) 


where M p « 10 19 GeU is the Planck scale. It corresponds to a length l p ps 10~ 33 cm. The 
Planck scale is thought to be the smallest length scale that makes sense: beyond this 
quantum gravity effects become important and it’s no longer clear that the concept 
of spacetime makes sense. The largest length scale we can talk of is the size of the 
cosmological horizon, roughly lO 60 ^. 


Observable 

Universe 

~ 20 billion light years 

y 

length _ 

10‘ 33 eV 


Earth 10 cm 


Cosmological 

Constant 


+ 


Atoms Nuclei 
—8 -13 

10 cm 10 cm 


— LHC 


Planck Scale 
10 33 cm 


Energy 


10 eV 


10 n -10 12 eV 
= 1 TeV 


10 28 eV 
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Figure 2: Energy and Distance Scales in the Universe 


Some useful scales in the universe are shown in the figure. This is a logarithmic plot, 
with energy increasing to the right and, correspondingly, length increasing to the left. 
The smallest and largest scales known are shown on the figure, together with other 
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relevant energy scales. The standard model of particle physics is expected to hold np 
to about the TeV. This is precisely the regime that is currently being probed by the 
Large Hadron Collider (LHC) at CERN. There is a general belief that the framework 
of quantum held theory will continue to hold to energy scales only slightly below the 
Planck scale — for example, there are experimental hints that the coupling constants 
of electromagnetism, and the weak and strong forces unify at around 10 18 GeV. 

For comparison, the rough masses of some elementary (and not so elementary) par¬ 
ticles are shown in the table, 


Particle 

Mass 

neutrinos 

~ 10" 2 eV 

electron 

0.5 MeV 

Muon 

100 MeV 

Pions 

140 MeV 

Proton, Neutron 

1 GeV 

Tau 

2 GeV 

W,Z Bosons 

80-90 GeV 

Higgs Boson 

125 GeV 




1. Classical Field Theory 


In this first section we will discuss various aspects of classical fields. We will cover only 
the bare minimum ground necessary before turning to the quantum theory, and will 
return to classical held theory at several later stages in the course when we need to 
introduce new ideas. 

1.1 The Dynamics of Fields 

A field is a quantity defined at every point of space and time (x,t). While classical 
particle mechanics deals with a finite number of generalized coordinates q a (t), indexed 
by a label a, in held theory we are interested in the dynamics of helds 

( 1 . 1 ) 

where both a and x are considered as labels. Thus we are dealing with a system with an 
infinite number of degrees of freedom — at least one for each point x in space. Notice 
that the concept of position has been relegated from a dynamical variable in particle 
mechanics to a mere label in held theory. 

An Example: The Electromagnetic Field 

The most familiar examples of helds from classical physics are the electric and magnetic 
helds, E(x,t) and B(x,t). Both of these helds are spatial 3-vectors. In a more sophis¬ 
ticated treatement of electromagnetism, we derive these two 3-vectors from a single 
4-component held A^^x, t ) = (0, A) where jx — 0,1, 2, 3 shows that this held is a vector 
in spacetime. The electric and magnetic helds are given by 

_ g 4 

E = —V0 — — and B = V x A (1.2) 

which ensure that two of Maxwell’s equations, V • B — 0 and dB/dt = —V x E, hold 
immediately as identities. 

The Lagrangian 

The dynamics of the held is governed by a Lagrangian which is a function of 0(T, t), 
4>(x,t) and V0(T, t). In all the systems we study in this course, the Lagrangian is of 
the form, 

L{t) = I dfix £(0 a ,<9 M 0 a ) (1.3) 
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where the official name for C is the Lagrangian density, although everyone simply calls 
it the Lagrangian. The action is, 

rt 2 

S= dt 

Jt i 

Recall that in particle mechanics L depends on q and q, but not q. In held theory 
we similarly restrict to Lagrangians C depending on 0 and 0, and not 0. In principle, 
there’s nothing to stop C depending on V0, V 2 0, V 3 0, etc. However, with an eye to 
later Lorentz invariance, we will only consider Lagrangians depending on V0 and not 
higher derivatives. Also we will not consider Lagrangians with explicit dependence on 
aft; all such dependence only comes through 0 and its derivatives. 



We can determine the equations of motion by the principle of least action. We vary 
the path, keeping the end points fixed and require SS = 0, 


5S = J d A x 
= [ d A x 


dC 


d(pa + 


dC 


dC 

9(p a 


d„ 


<9(<9 M 0a) 

dC 


S(d„<f) a ) 
d(pa + C0 


DC 


d 0 C 


(1.5) 


\d(d^(p a )J\ ' J ' ra ' ^ \d(du,(/) a ) 

The last term is a total derivative and vanishes for any 5(p a {x,t) that decays at spatial 
infinity and obeys 50 a (x, H) = 8 <p a (x,t 2 ) = 0. Requiring 5S = 0 for all such paths 
yields the Euler-Lagrange equations of motion for the fields 0 a , 

dC \ dC 


d„ 


= 0 


X d{d^ a )J d(p a 

1.1.1 An Example: The Klein-Gordon Equation 

Consider the Lagrangian for a real scalar held 0(x, t), 

C = \ if v d^cpdvcf) - \m 2 4> 2 


( 1 . 6 ) 


1 


(1.7) 


- -(V0) - \m 


1^2 ±2 


2V 2 

where we are using the Minkowski space metric 


+i 


r r = vu* = 


( 1 . 8 ) 


Comparing (1.7) to the usual expression for the Lagrangian L = T — V, we identify the 
kinetic energy of the held as 


'3^ 1 12 


T = ft) 


(1.9) 













and the potential energy of the held as 


1/ = j d 3 x l(V</>) 2 + \m 2 <p 2 (1.10) 

The first term in this expression is called the gradient energy, while the phrase “poten¬ 
tial energy”, or just “potential”, is usually reserved for the last term. 

To determine the equations of motion arising from (1.7), we compute 
fjr nr 

9t = - m 4, and (L11) 

The Euler-Lagrange equation is then 

4 > — V 2 0 + m 2 (j) = 0 (1.12) 

which we can write in relativistic form as 

d^d^cj) + m 2 (j) = 0 (1.13) 

This is the Klein-Gordon Equation. The Laplacian in Minkowski space is sometimes 
denoted by □. In this notation, the Klein-Gordon equation reads D0 + m 2 (f) = 0. 


An obvious generalization of the Klein-Gordon equation comes from considering the 
Lagrangian with arbitrary potential V(<p), 


C = \d^d^ - V(<P) 




dV 

d<p 


0 


(1.14) 


1.1.2 Another Example: First Order Lagrangians 


We could also consider a Lagrangian that is linear in time derivatives, rather than 
quadratic. Take a complex scalar held ip whose dynamics is defined by the real La¬ 
grangian 


£ = ^(ip*ip — , 0*'0) — • W — mil)*ip 


(1.15) 


We can determine the equations of motion by treating ip and ip* as independent objects, 
so that 


d£ 

dip * 


-ip — mip and 


d£ 
dip* 


—ip and 
2 


d£ 

d\7ip* 


-Wip 


(1.16) 


This gives us the equation of motion 

dib 

i— = — X7 2 ip + mip (1.17) 

This looks very much like the Schrodinger equation. Except it isn’t! Or, at least, the 
interpretation of this equation is very different: the held ip is a classical held with none 
of the probability interpretation of the wavefunction. We’ll come back to this point in 
Section 2.8. 
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The initial data required on a Cauchy surface differs for the two examples above. 
When £ ~ </> 2 , both (p and <p must be specified to determine the future evolution; 
however when £ ~ if;*ip, only ip and ip* are needed. 


1.1.3 A Final Example: Maxwell’s Equations 

We may derive Maxwell’s equations in the vacuum from the Lagrangian, 

C = -\ (<VU (3M") + (1.18) 

Notice the funny minus signs! This is to ensure that the kinetic terms for A* are positive 
using the Minkowski space metric (1.8), so £ ~ -(A 2 . The Lagrangian (1.18) has no 
kinetic term Aq for A 0 . We will see the consequences of this in Section 6. To see that 
Maxwell’s equations indeed follow from (1.18), we compute 

d&Aj) = -SM- + (df.An IT (1.19) 

from which we may derive the equations of motion, 

d, + d u (d P A p ) = - FA*) = -d,F^ (1.20) 

where the field strength is defined by F^ u = d tl A v — d u A /1 . You can check using (1.2) 
that this reproduces the remaining two Maxwell’s equations in a vacuum: V ■ E = 0 
and dE/dt = V x B. Using the notation of the field strength, we may rewrite the 
Maxwell Lagrangian (up to an integration by parts) in the compact form 

C = -\ F ijm F>^ ( 1 . 21 ) 


1.1.4 Locality, Locality, Locality 

In each of the examples above, the Lagrangian is local. This means that there are no 
terms in the Lagrangian coupling cp(x, t ) directly to cp(y, t ) with x j- y. For example, 
there are no terms that look like 


L 


d 3 xd 3 y (p(x)(p(y) 


( 1 . 22 ) 


A priori, there’s no reason for this. After all, x is merely a label, and we’re quite 
happy to couple other labels together (for example, the term c^Aq <9oA 3 in the Maxwell 
Lagrangian couples the y = 0 field to the y = 3 field). But the closest we get for the 
x label is a coupling between (p{x) and <p(x + Sx) through the gradient term (V0) 2 . 
This property of locality is, as far as we know, a key feature of all theories of Nature. 
Indeed, one of the main reasons for introducing field theories in classical physics is to 
implement locality. In this course, we will only consider local Lagrangians. 
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1.2 Lorentz Invariance 


The laws of Nature are relativistic, and one of the main motivations to develop quantum 
field theory is to reconcile quantum mechanics with special relativity. To this end, we 
want to construct field theories in which space and time are placed on an equal footing 
and the theory is invariant under Lorentz transformations, 

x» —► (x'Y = A %x v (1.23) 


where A^ v satisfies 


A M CT rf T A u t = r t 


(1.24) 


For example, a rotation by 9 about the a; 3 -axis, and a boost by v < 1 along the ad-axis 
are respectively described by the Lorentz transformations 


/l 0 0 o\ 

0 cos 6 — sin 9 0 
0 sin 6 cos 6 0 
\0 0 01 / 


and A'd 


/ /y 0 0 \ 

—jv 7 0 0 

0 0 10 

0 0 0 1 / 


(1.25) 


with 7 = 1 /v/1 — v 2 . The Lorentz transformations form a Lie group under matrix 
multiplication. You’ll learn more about this in the “Symmetries and Particle Physics” 
course. 


The Lorentz transformations have a representation on the fields. The simplest ex¬ 
ample is the scalar field which, under the Lorentz transformation x —y Ax, transforms 
as 


4 >{x) —» <p'(x) = (f ){A 1 x) (1-26) 

The inverse A -1 appears in the argument because we are dealing with an active trans¬ 
formation in which the field is truly shifted. To see why this means that the inverse 
appears, it will suffice to consider a non-relativistic example such as a temperature field. 
Suppose we start with an initial field 4>(x) which has a hotspot at, say, x = (1,0,0). 
After a rotation x —y Rx about the z-axis, the new field (j)'{x) will have the hotspot at 
x — (0,1, 0). If we want to express </>'(x) in terms of the old field q i, we need to place 
ourselves at x — ( 0 , 1 , 0 ) and ask what the old field looked like where we’ve come from 
at f? _ 1 (0,1,0) = (1,0,0). This R is the origin of the inverse transformation. (If we 
were instead dealing with a passive transformation in which we relabel our choice of 
coordinates, we would have instead <j)(x) —>■ (f>'(x) = <j)(Ax)). 
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The definition of a Lorentz invariant theory is that if <p(x) solves the equations of 
motion then (j)(A~ 1 x ) also solves the equations of motion. We can ensure that this 
property holds by requiring that the action is Lorentz invariant. Let’s look at our 
examples: 


Example 1: The Klein-Gordon Equation 

For a real scalar field we have 0(x) —y cj>'(x) = 0(A~ 1 a:). The derivative of the scalar 
field transforms as a vector, meaning 

(d^)(x) -)■ (A 

where y = A _1 :r. This means that the derivative terms in the Lagrangian density 
transform as 

Cdenv(x) = d^(x)d u (j)(x)rf v —> (A {K~ x Y v {d a <j>){y) rf v 

= {d P <t>){y) (d*(f))(y)v pa 

= C deriv {y) (1.27) 

The potential terms transform in the same way, with (p 2 (x ) —> 4> 2 (y). Putting this all 
together, we find that the action is indeed invariant under Lorentz transformations, 

S — j d 4 x C{x) —> J (Yx £{y) = j d A y C{y) = S (1-28) 

where, in the last step, we need the fact that we don’t pick up a Jacobian factor when 
we change integration variables from f d 4 x to f d 4 y. This follows because det A = 1. 
(At least for Lorentz transformation connected to the identity which, for now, is all we 
deal with). 


Example 2: First Order Dynamics 

In the first-order Lagrangian (1.15), space and time are not on the same footing. (£ 
is linear in time derivatives, but quadratic in spatial derivatives). The theory is not 
Lorentz invariant. 

In practice, it’s easy to see if the action is Lorentz invariant: just make sure all 
the Lorentz indices /i = 0,1, 2, 3 are contracted with Lorentz invariant objects, such 
as the metric rj llu . Other Lorentz invariant objects you can use include the totally 
antisymmetric tensor e^ vprj and the matrices that we will introduce when we come 
to discuss spinors in Section 4. 
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Example 3: Maxwell’s Equations 

Under a Lorentz transformation A^(x) —y A IJ U A u ( y A~ 1 x). Yon can check that Maxwell’s 
Lagrangian (1.21) is indeed invariant. Of course, historically electrodynamics was the 
first Lorentz invariant theory to be discovered: it was found even before the concept of 
Lorentz invariance. 

1.3 Symmetries 

The role of symmetries in held theory is possibly even more important than in particle 
mechanics. There are Lorentz symmetries, internal symmetries, gauge symmetries, 
supersymmetries.... We start here by recasting Noether’s theorem in a held theoretic 
framework. 

1.3.1 Noether’s Theorem 

Every continuous symmetry of the Lagrangian gives rise to a conserved current 
such that the equations of motion imply 

<9^ = 0 (1.29) 

or, in other words, dj°/dt + V ■ j — 0. 

A Comment: A conserved current implies a conserved charge (), defined as 


Q = [ d 3 x j 0 

J R3 

which one can immediately see by taking the time derivative, 


dQ = f af 

n J R 3 at 



(1.30) 


(1.31) 


assuming that j —> 0 sufficiently quickly as |x| —» oo. However, the existence of a 
current is a much stronger statement than the existence of a conserved charge because 
it implies that charge is conserved locally. To see this, we can define the charge in a 
finite volume V, 


Qv = [ d 3 x j 0 

Jv 

Repeating the analysis above, we find that 

d -w = -j v d 3 xVj= -Jj - i3 


(1.32) 


(1.33) 
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where A is the area bounding V and we have used Stokes’ theorem. This equation 
means that any charge leaving V must be accounted for by a flow of the current 3- 
vector j out of the volume. This kind of local conservation of charge holds in any local 
held theory. 

Proof of Noether’s Theorem: We’ll prove the theorem by working infinitesimally. 
We may always do this if we have a continuous symmetry. We say that the transfor¬ 
mation 


S(p a (x) = X a ((j )) 

is a symmetry if the Lagrangian changes by a total derivative, 


(1.34) 


SC = d„F“ 


(1.35) 


for some set of functions E M (0). To derive Noether’s theorem, we first consider making 
an arbitrary transformation of the fields S(p a . Then 


SC 



dC 

<9(<9/A) 


dp{6<l> a ) 


dC _ dC 
d(p a M d(dM 


S(f) a + dp 



(1.36) 


When the equations of motion are satisfied, the term in square brackets vanishes. So 
we’re left with 


SC 


dp 



(1.37) 


But for the symmetry transformation Sc/) a = X a (0), we have by definition SC = dpF C 
Equating this expression with (1.37) gives us the result 


dpj» = 0 with ^ = __x a (0)-^(0) (1.38) 

1.3.2 An Example: Translations and the Energy-Momentum Tensor 

Recall that in classical particle mechanics, invariance under spatial translations gives 
rise to the conservation of momentum, while invariance under time translations is 
responsible for the conservation of energy. We will now see something similar in held 
theories. Consider the infinitesimal translation 


x v y x v C =>■ (j) a {x) ->■ 4 > a (x) + e u d v (f> a (x) 


(1.39) 
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(where the sign in the field transformation is plus, instead of minus, because we’re doing 
an active, as opposed to passive, transformation). Similarly, once we substitute a spe¬ 
cific field configuration (j>{x) into the Lagrangian, the Lagrangian itself also transforms 
as 


C(x) —> jC(x) + e v d v C(x) 


(1.40) 


Since the change in the Lagrangian is a total derivative, we may invoke Noether’s 
theorem which gives us four conserved currents one for each of the translations 

e v with u — 0,1,2, 3, 

dC 

= oToNTT Wa ~ 6ZC = T1 (1.41) 

T p v is called the energy-momentum tensor. It satisfies 


d.T'l = 0 


The four conserved quantities are given by 


E = d 3 x T 


oo 


and P l = / d 3 x T 


0z 


(1.42) 

(1.43) 


where E is the total energy of the field configuration, while P l is the total momentum 
of the field configuration. 


An Example of the Energy-Momentum Tensor 

Consider the simplest scalar field theory with Lagrangian (1.7). From the above dis¬ 
cussion, we can compute 

T 1 " = &*<!> (1.44) 

One can verify using the equation of motion for (j) that this expression indeed satisfies 
d/jT 111 ' = 0. For this example, the conserved energy and momentum are given by 

E — I d 3 x l<j) 2 + i(V</>) 2 + \m 2 ct) 2 (1.45) 

P l = J d 3 x (1-46) 

Notice that for this example, T^ came out symmetric, so that T^ u = This 

won’t always be the case. Nevertheless, there is typically a way to massage the energy 
momentum tensor of any theory into a symmetric form by adding an extra term 

& jy = T IJy + d p V^ v (1.47) 

where is some function of the fields that is anti-symmetric in the first two indices so 
Y ppv = — Y ppy . This guarantees that d p d p Y pp,l/ = 0 so that the new energy-momentum 
tensor is also a conserved current. 
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A Cute Trick 


One reason that you may want a symmetric energy-momentum tensor is to make con¬ 
tact with general relativity: such an object sits on the right-hand side of Einstein’s 
held equations. In fact this observation provides a quick and easy way to determine a 
symmetric energy-momentum tensor. Firstly consider coupling the theory to a curved 
background spacetime, introducing an arbitrary metric g lw {x) in place of r/ /w , and re¬ 
placing the kinetic terms with suitable covariant derivatives using “minimal coupling”. 
Then a symmetric energy momentum tensor in the hat space theory is given by 

0 »„ 2 au^gC) 

V~9 '%» 

It should be noted however that this trick requires a little more care when working 
with spinors. 

1.3.3 Another Example: Lorentz Transformations and Angular Momentum 

In classical particle mechanics, rotational invariance gave rise to conservation of angular 
momentum. What is the analogy in held theory? Moreover, we now have further 
Lorentz transformations, namely boosts. What conserved quantity do they correspond 
to? To answer these questions, we hrst need the inhnitesimal form of the Lorentz 
transformations 


(1.48) 


A“„ = K + < (1-49) 

where okj, is inhnitesimal. The condition (1.24) for A to be a Lorentz transformation 
becomes 


{K + ^aWr + u v T ) rf T = tT 

=> + u Uf * = 0 (1.50) 

So the inhnitesimal form of the Lorentz transformation must be an anti-symmetric 
matrix. As a check, the number of different 4x4 anti-symmetric matrices is 4 x 3/2 = 6, 
which agrees with the number of different Lorentz transformations (3 rotations + 3 
boosts). Now the transformation on a scalar held is given by 

4>(x) —» (f>\x) = (f>( A" 1 ^) 

= — oZ/aA) 

= 0(x M ) — Cck/ x u d^(j)(x) (1.51) 
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from which we see that 


Sep = -a/ „x v d 

By the same argument, the Lagrangian density transforms as 


(1.52) 


SC = (1.53) 

where the last equality follows because cU), = 0 due to anti-symmetry. Once again, 
the Lagrangian changes by a total derivative so we may apply Noether’s theorem (now 
with F M = —t o p v x v C) to find the conserved current 


r = 


dC 


d(d^) 


u p v x v d p (f) + uj p v x v C 


= - u p 


dC 


[d(dp(j)) 


x v d p (j> - 8 P x v C 


-u p v ry 


(1.54) 


Unlike in the previous example, I’ve left the infinitesimal choice of u p u in the expression 
for this current. But really, we should strip it out to give six different currents, i.e. one 
for each choice of u p u . We can write them as 


(J p ) p,T = x p T pa - x a T pp 


(1.55) 


which satisfy d ll (J IJ ') p(T = 0 and give rise to 6 conserved charges. For p,a = 1,2,3, 
the Lorentz transformation is a rotation and the three conserved charges give the total 
angular momentum of the field. 

Q ij = I d 3 x (x ,: T 0 l - x j T 0i ) (1.56) 

But what about the boosts? In this case, the conserved charges are 

Q 0i = I d 3 x (x°T 0i - FT 00 ) (1.57) 

The fact that these are conserved tells us that 

0 = ^= [ d 3 xT 0l +t [ d 3 x d ^-^- [ d 3 x FT 00 

dt J J dt dt J 

= P l + t - 4- f d 3 x i ! T 00 (1.58) 

dt dt J y J 

But we know that P l is conserved, so dP l /dt = 0, leaving us with the following conse¬ 
quence of invariance under boosts: 

7 j* 

— / d 3 x x l T m = constant (1.59) 

This is the statement that the center of energy of the field travels with a constant 
velocity. It’s kind of like a field theoretic version of Newton’s first law but, rather 
surprisingly, appearing here as a conservation law. 
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1.3.4 Internal Symmetries 


The above two examples involved transformations of spacetime, as well as transforma¬ 
tions of the field. An internal symmetry is one that only involves a transformation of 
the fields and acts the same at every point in spacetime. The simplest example occurs 
for a complex scalar field ip(x) = (4>i(x) + i(p 2 (x)) / \/ 2. We can build a real Lagrangian 
by 


C = d fi rd^-V( |^| 2 ) 


(1.60) 


where the potential is a general polynomial in \ip\ 2 = ip*^- To find the equations of 
motion, we could expand ip in terms of (p i and (p 2 and work as before. However, it’s 
easier (and equivalent) to treat ip and ip* as independent variables and vary the action 
with respect to both of them. For example, varying with respect to ip* leads to the 
equation of motion 


dnd^i> + 


dV (ip*' l l J ) 


= 0 


(1.61) 


m r 1 

The Lagrangian has a continuous symmetry which rotates (p\ and (p 2 or, equivalently, 
rotates the phase of ip: 


ip —> e za ip or Sip = iaip 


(1.62) 


where the latter equation holds with a infinitesimal. The Lagrangian remains invariant 
under this change: SC = 0. The associated conserved current is 

= i(d»ip*)ip - iV^ip) (1.63) 


We will later see that the conserved charges arising from currents of this type have 
the interpretation of electric charge or particle number (for example, baryon or lepton 
number). 


Non-Abelian Internal Symmetries 

Consider a theory involving N scalar fields (p a , all with the same mass and the La¬ 
grangian 

2 i2 

m <p a -9 

a =1 a =1 

In this case the Lagrangian is invariant under the non-Abelian symmetry group G = 
SO(N). (Actually O(N) in this case). One can construct theories from complex fields 
in a similar manner that are invariant under an SU(N) symmetry group. Non-Abelian 
symmetries of this type are often referred to as global symmetries to distinguish them 
from the “local gauge” symmetries that you will meet later. Isospin is an example of 
such a symmetry, albeit realized only approximately in Nature. 


C = \ £ - \ J2 



(1.64) 
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Another Cute Trick 

There is a quick method to determine the conserved current associated to an internal 
symmetry Sep = ot(p for which the Lagrangian is invariant. Here, a is a constant real 
number. (We may generalize the discussion easily to a non-Abelian internal symmetry 
for which a becomes a matrix). Now consider performing the transformation but where 
a depends on spacetime: a = a(x). The action is no longer invariant. However, the 
change must be of the form 

SC = (<V*) h^cP) (1.65) 

since we know that SC = 0 when a is constant. The change in the action is therefore 

SS = j d A x SC = — J d A x a(x) d^lP (1.66) 

which means that when the equations of motion are satisfied (so SS = 0 for all varia¬ 
tions, including 5<p = a(x)<p) we have 

d t Jd l = 0 (1.67) 

We see that we can identify the function h fl = as the conserved current. This way 
of viewing things emphasizes that it is the derivative terms, not the potential terms, 
in the action that contribute to the current. (The potential terms are invariant even 
when a = a(x)). 


1.4 The Hamiltonian Formalism 


The link between the Lagrangian formalism and the quantum theory goes via the path 
integral. In this course we will not discuss path integral methods, and focus instead 
on canonical quantization. For this we need the Hamiltonian formalism of field theory. 
We start by defining the momentum 7T a (x) conjugate to <p a (x ), 


7r a (a:) = 



( 1 . 68 ) 


The conjugate momentum 7 r a (x) is a function of x, just like the field <p a {x) itself. It 
is not to be confused with the total momentum P l defined in (1.43) which is a single 
number characterizing the whole field configuration. The Hamiltonian density is given 

by 


U = 7 T a (x)<p a (x) - C(x) 


(1.69) 


where, as in classical mechanics, we eliminate <p a (x) in favour of n a (x) everywhere in 
H. The Hamiltonian is then simply 


H = 



(1.70) 
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An Example: A Real Scalar Field 

For the Lagrangian 



1 

2 


W ) 2 - v[<t>) 


(1.71) 


the momentum is given by 7 t — (j), which gives us the Hamiltonian, 


H 


I 


d?x lir 2 + i(V^.) 2 + V(<t>) 


(1.72) 


Notice that the Hamiltonian agrees with the definition of the total energy (1.45) that 


we get from applying Noether’s theorem for time translation invariance. 

In the Lagrangian formalism, Lorentz invariance is clear for all to see since the action 
is invariant under Lorentz transformations. In contrast, the Hamiltonian formalism is 
not manifestly Lorentz invariant: we have picked a preferred time. For example, the 
equations of motion for 4>(x) = t ) arise from Hamilton’s equations, 



(1.73) 


which, unlike the Euler-Lagrange equations (1.6), do not look Lorentz invariant. Nev¬ 
ertheless, even though the Hamiltonian framework doesn’t look Lorentz invariant, the 
physics must remain unchanged. If we start from a relativistic theory, all final answers 


must be Lorentz invariant even if it’s not manifest at intermediate steps. We will pause 


at several points along the quantum route to check that this is indeed the case. 
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2. Free Fields 


“The career of a young theoretical physicist consists of treating the harmonic 

oscillator in ever-increasing levels of abstraction.” n . , „ , 

bidney Coleman 


2.1 Canonical Quantization 

In quantum mechanics, canonical quantization is a recipe that takes us from the Hamil¬ 
tonian formalism of classical dynamics to the quantum theory. The recipe tells us to 
take the generalized coordinates q a and their conjugate momenta p a and promote them 
to operators. The Poisson bracket structure of classical mechanics morphs into the 
structure of commutation relations between operators, so that, in units with h — 1, 

[Qa,Qb] = \p a ,P b } = 0 

[Qa,p b ] = iS b a (2.1) 

In field theory we do the same, now for the field f> a {x) and its momentum conjugate 
n b (x). Thus a quantum field is an operator valued function of space obeying the com¬ 
mutation relations 


[</>»(£), <M$] = [n a (x),n b (y)\ = 0 

[fia(x),TT b (y)]=i5 (3 \x-y)5 b ( 2 . 2 ) 


Note that we’ve lost all track of Lorentz invariance since we have separated space x 
and time t. We are working in the Schrodinger picture so that the operators <f> a {x) and 
n a (x) do not depend on time at all — only on space. All time dependence sits in the 
states | if) which evolve by the usual Schrodinger equation 


■ d\ip) 

l -r- 

dt 


H\ip) 


(2.3) 


We aren’t doing anything different from usual quantum mechanics; we’re merely apply¬ 
ing the old formalism to fields. Be warned however that the notation \if) for the state 
is deceptively simple: if you were to write the wavefunction in quantum field theory, it 
would be a functional , that is a function of every possible configuration of the field qb. 


The typical information we want to know about a quantum theory is the spectrum of 
the Hamiltonian H. In quantum field theories, this is usually very hard. One reason for 
this is that we have an infinite number of degrees of freedom — at least one for every 
point x in space. However, for certain theories — known as free theories — we can find 
a way to write the dynamics such that each degree of freedom evolves independently 


21 



from all the others. Free field theories typically have Lagrangians which are quadratic 
in the fields, so that the equations of motion are linear. For example, the simplest 
relativistic free theory is the classical Klcin-Gordon (KG) equation for a real scalar 
field (j)(x, t ), 

d^d^cj) + m 2 (j) = 0 (2.4) 

To exhibit the coordinates in which the degrees of freedom decouple from each other, 
we need only take the Fourier transform, 

$&t) = J (2.5) 

Then cp(p, t) satisfies 

+ ( P 2 + 0(P> t) = 0 ( 2 -6) 

Thus, for each value of p, (p(p, t ) solves the equation of a harmonic oscillator vibrating 
at frequency 

up = + \/p 2 + m 2 (2.7) 

We learn that the most general solution to the KG equation is a linear superposition of 
simple harmonic oscillators, each vibrating at a different frequency with a different am¬ 
plitude. To quantize (p(x, t ) we must simply quantize this infinite number of harmonic 
oscillators. Let’s recall how to do this. 

2.1.1 The Simple Harmonic Oscillator 

Consider the quantum mechanical Hamiltonian 

H = \p 2 + \u 2 q 2 (2.8) 

with the canonical commutation relations [q,p\ — %■ To find the spectrum we define 
the creation and annihilation operators (also known as raising/lowering operators, or 
sometimes ladder operators) 



which can be easily inverted to give 

9=^^(a + a f ) , p= a f ) (2.10) 
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Substituting into the above expressions we find 

[a, a'] = 1 


( 2 . 11 ) 


while the Hamiltonian is given by 

H = \u(aal + ada) 

= w(a t a+§) (2.12) 

One can easily confirm that the commutators between the Hamiltonian and the creation 
and annihilation operators are given by 

[H,al]=uja) and [H J a\ = —ua (2-13) 

These relations ensure that a and at take us between energy eigenstates. Let | E) be 
an eigenstate with energy E, so that H\E) = E\E). Then we can construct more 
eigenstates by acting with a and ad, 

Hat \E) = (E + u)a! \E) , Ha \E) = (E - u)a\E) (2.14) 

So we find that the system has a ladder of states with energies 

. . ., E — uj, E, E + uj, E + 2cn, ... (2-15) 

If the energy is bounded below, there must be a ground state |0) which satisfies a |0) = 0. 
This has ground state energy (also known as zero point energy), 

H |0) = \uj |0) (2.16) 

Excited states then arise from repeated application of ad, 

| n) = (cd) n |0) with H \n) = (n + \)uj \n) (2-17) 

where I’ve ignored the normalization of these states so, (n\ n) ^ 1. 


2.2 The Free Scalar Field 


We now apply the quantization of the harmonic oscillator to the free scalar field. We 
write (f) and n as a linear sum of an infinite number of creation and annihilation oper¬ 
ators at and ap, indexed by the 3-momentum p, 


cj)(x) 


7r(x) 


d 3 p 


(2vr) 3 
d 3 p 


(2tt) 



a f e?* 3 + ale-v-* 


a^P * - ale-^* 


(2.18) 

(2.19) 
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Claim: The commutation relations for 0 and n are equivalent to the following com¬ 
mutation relations for ap and at 


[0(f), 0(y)] = [ir(x),7r (£)] = 0 
[<t>{x),n{y)\ = iS^(x-y) 


[a?, ap] = [ajf, aj] = 0 

[ap,a t] = (2vr ) 3 6^ 3 \p- q) 


( 2 . 20 ) 


Proof: We’ll show this just one way. Assume that [a-, at] = (27T ) 3 8^ 3 \p — q). Then 
d 3 p d 3 q (—i) fuTq 




(2tt) 6 2 

d 3 p (— i) 


— K a S-] ***_** + [4’ 4 e~ ifS+i ™ 


J (27r) 3 2 

= id^ 3 \x — y) 


_ e ip-(x-y) _ e ip-{y-x)\ 


( 2 . 21 ) 

□ 


The Hamiltonian 

Let’s now compute the Hamiltonian in terms of ap and at. We have 


H — - d 3 x 7r 2 + (V0) 2 + m 2 0 2 


1 /“ d 3 x d 3 p d 3 q 
2. 


(27r) 6 
1 


4 


UJfPJ, 


p^q 


a-e 


„-ip-x 


ale~^)(a^ g - ate -4 **) 


V 




(ipa f e ifg - ip ale~^ S ) • {iqa.e^ - iqale “ 4 * 2 ) 


2 -y/ LOpUJ q 

i r d 3 p i 

4 J (27r) 3 a; ; 


(a^e®^ -2 + ate-**)^ 4 ** + ate"" 4 * 2 ) 


(—+ p 2 + m 2 )(apd-p + atal-) + (cut + p 2 + m 2 )(ap-at + at a^) 


,t , „t 


V 


where in the second line we’ve used the expressions for 0 and 7r given in (2.18) and 
(2.19); to get to the third line we’ve integrated over d 3 x to get delta-functions 5^(p±q) 
which, in turn, allow us to perform the d 3 q integral. Now using the expression for the 
frequency cat = p 2 + m 2 , the first term vanishes and we’re left with 


H = 1 


f d 3 p 
(2^ 
d 3 p 


«p4+ 4 a P 


(2n 


i.3 U P 


4 a P+ |(2vr) 3 5 (3) (0) 


( 2 . 22 ) 
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Hmmmm . We’ve found a delta-function, evaluated at zero where it has its infinite 
spike. Moreover, the integral over up diverges at large p. What to do? Let’s start by 
looking at the ground state where this infinity first becomes apparent. 


2.3 The Vacuum 


Following onr procedure for the harmonic oscillator, let’s define the vacuum |0) by 
insisting that it is annihilated by all a~, 

|0) = 0 Vp (2.23) 


With this definition, the energy Eq of the ground state comes from the second term in 

( 2 . 22 ), 


H\0) = E 0 \0) 



d 3 p \upS^( 0) 


0 ) 


oo |0) 


(2.24) 


The subject of quantum held theory is rife with infinities. Each tells us something 
important, usually that we’re doing something wrong, or asking the wrong question. 
Let’s take some time to explore where this infinity comes from and how we should deal 
with it. 


In fact there are two different oo’s lurking in the expression (2.24). The first arises 
because space is infinitely large. (Infinities of this type are often referred to as infra-red 
divergences although in this case the oo is so simple that it barely deserves this name). 
To extract out this infinity, let’s consider putting the theory in a box with sides of 
length L. We impose periodic boundary conditions on the held. Then, taking the limit 
where L —)■ oo, we get 

r L l 2 pL /2 

(2tt) 3 5 (3) (0) = lim / d 3 x e™* L_= lim / d 3 x = V (2.25) 

L^oo J_ L j 2 p L- s-oo J _ L / 2 

where V is the volume of the box. So the 5(0) divergence arises because we’re computing 
the total energy, rather than the energy density Sq. To find £ 0 we can simply divide by 
the volume, 


£ - Eo - 
£ °-v~ 


d 3 p 1 


(2tt) 3 2 




(2.26) 


which is still infinite. We recognize it as the sum of ground state energies for each 
harmonic oscillator. But So —> oo due to the \p\ —> oo limit of the integral. This is 
a high frequency — or short distance — infinity known as an ultra-violet divergence. 
This divergence arises because of our hubris. We’ve assumed that our theory is valid 
to arbitrarily short distance scales, corresponding to arbitrarily high energies. This is 
clearly absurd. The integral should be cut-off at high momentum in order to reflect 
the fact that our theory is likely to break down in some way. 
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We can deal with the infinity in (2.24) in a more practical way. In physics we’re 
only interested in energy differences. There’s no way to measure Eq directly, so we can 
simply redefine the Hamiltonian by subtracting off this infinity, 


H 


d 3 p f 

(2tt) 3 UJ P a p a v 


(2.27) 


so that, with this new definition, H |0) = 0. In fact, the difference between this Hamilto¬ 
nian and the previous one is merely an ordering ambiguity in moving from the classical 
theory to the quantum theory. For example, if we defined the Hamiltonian of the har¬ 
monic oscillator to be H = (l/2)(ojq — ip)(ojq + ip), which is classically the same as our 
original choice, then upon quantization it would naturally give H = oja^a as in (2.27). 
This type of ordering ambiguity arises a lot in field theories. We’ll come across a number 
of ways of dealing with it. The method that we’ve used above is called normal ordering. 


Definition: We write the normal ordered string of operators ... 0 n (x n ) as 


: 0i(xi)...0 n (x„) : 


(2.28) 


It is defined to be the usual product with all annihilation operators a- placed to the 
right. So, for the Hamiltonian, we could write (2.27) as 

rr [ f 

:H Wf * f r 

In the remainder of this section, we will normal order all operators in this manner. 


(2.29) 


2.3.1 The Cosmological Constant 

Above I wrote “there’s no way to measure E 0 directly”. There is a BIG caveat here: 
gravity is supposed to see everything! The sum of all the zero point energies should 
contribute to the stress-energy tensor that appears on the right-hand side of Einstein’s 
equations. We expect them to appear as a cosmological constant A = E 0 /V, 

Ruv - \Rg^v = -8 ttGT\ v + Ag llu (2.30) 

Current observation suggests that 70% of the energy density in the universe has the 
properties of a cosmological constant with A ~ (10 -3 eH) 4 . This is much smaller than 
other scales in particle physics. In particular, the Standard Model is valid at least up 
to 10 12 eV. Why don’t the zero point energies of these fields contribute to A? Or, if 
they do, what cancels them to such high accuracy? This is the cosmological constant 
problem. No one knows the answer! 
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2.3.2 The Casimir Effect 


“I mentioned my results to Niels Bohr, during a walk. That is nice, he 
said, that is something new... and he mumbled something about zero-point 

energy. Hendrik Casimir 


Using the normal ordering prescription we can happily set E 0 = 0, while chanting 
the mantra that only energy differences can be measured. But we should be careful, for 
there is a situation where differences in the energy of vacuum fluctuations themselves 
can be measured. 

To regulate the infra-red divergences, we’ll make the x 1 direction periodic, with size 
L, and impose periodic boundary conditions such that 


0(x) = (f>{x + Ln) 


(2.31) 


with n = (1,0,0). We’ll leave y and z alone, but remember 
that we should compute all physical quantities per unit area 
A. We insert two reflecting plates, separated by a distance 
d <C L in the x 1 direction. The plates are such that they 
impose (f>(x) = 0 at the position of the plates. The presence of 
these plates affects the Fourier decomposition of the field and, 
in particular, means that the momentum of the field inside the 
plates is quantized as 

/ 77/7T \ 

V=y-j,Py,Pz) ne z + (2.32) 


d 

<—> 



Figure 3: 


For a massless scalar field, the ground state energy between the plates is 


m 

A 



dpydp z 1 / /rvn\ 

(2vr) 2 2\\d) 


+ P 2 y +P 2 z 


(2.33) 


while the energy outside the plates is E(L — d). The total energy is therefore 


E = E(d ) + E(L - d) 


(2.34) 


which - at least naively - depends on d. If this naive guess is true, it would mean 
that there is a force on the plates due to the fluctuations of the vacuum. This is the 
Casimir force, first predicted in 1948 and observed 10 years later. In the real world, the 
effect is due to the vacuum fluctuations of the electromagnetic field, with the boundary 
conditions imposed by conducting plates. Here we model this effect with a scalar. 
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But there’s a problem. E is infinite! What to do? The problem comes from the 
arbitrarily high momentum modes. We could regulate this in a number of different 
ways. Physically one could argue that any real plate cannot reflect waves of arbitrarily 
high frequency: at some point, things begin to leak. Mathematically, we want to find 
a way to neglect modes of momentum p a -1 for some distance scale a<d, known 
as the ultra-violet (UV) cut-off. One way to do this is to change the integral (2.33) to, 



which has the property that as a —» 0, we regain the full, infinite, expression (2.33). 
However (2.35) is finite, and gives us something we can easily work with. Of course, 
we made it finite in a rather ad-hoc manner and we better make sure that any physical 
quantity we calculate doesn’t depend on the UV cut-off a, otherwise it’s not something 
we can really trust. 


The integral (2.35) is do-able, but a little complicated. It’s a lot simpler if we look 
at the problem in d = 1 + 1 dimensions, rather than d — 3 + 1 dimensions. We’ll find 
that all the same physics is at play. Now the energy is given by 

OO 

E i+i( d ) = n ( 2 - 36 ) 

n= 1 

We now regulate this sum by introducing the UV cutoff a introduced above. This 
renders the expression finite, allowing us to start manipulating it thus, 


E 1+1 (d) ne~ an ^ d 



1 d 1 

2 da 1 - e~ aTC / d 

7r e an l d 

2d ( e an / d - l) 2 


where, in the 
energy, 


d 7T , „, 

= " 245 + ° {a ] 

last line, we’ve used the fact that We can now compute 


(2.37) 
the full 


E i+i 


E i+i (d) + Ei + i(L — d) 


L 

27ra 2 


- 28 - 


7T /l 1 \ 

24 \d + L-d) 


+ 0(a 2 ) 


(2.38) 











This is still infinite in the limit a —> 0, which is to be expected. However, the force is 
given by 


dE 1+1 _ Ti 
dd ~ 24 d? + ‘ ‘' 


(2.39) 


where the ... include terms of size d/L and a/d. The key point is that as we remove 
both the regulators, and take a —> 0 and L —> oo, the force between the plates remains 
finite. This is the Casimir force 2 . 


If we ploughed through the analogous calculation in d — 3 + 1 dimensions, and 
performed the integral (2.35), we would find the result 


IdE _ 7r 2 

A~dd ~ 480# 


(2.40) 


The true Casimir force is twice as large as this, due to the two polarization states of 
the photon. 


2.4 Particles 

Having dealt with the vacuum, we can now turn to the excitations of the field. It’s 
easy to verify that 

[H, at] = ujpai and [H, a^} = (2.41) 

which means that, just as for the harmonic oscillator, we can construct energy eigen¬ 
states by acting on the vacuum |0) with at. Let 

|p) = 4l°) (2.42) 

This state has energy 

H \p) = oop\p) with c Op = p 2 + m 2 (2.43) 

But we recognize this as the relativistic dispersion relation for a particle of mass m and 
3-momentum p. 


El = p 2 + m 2 (2.44) 

2 The number 24 that appears in the denominator of the one-dimensional Casimir force plays a more 
famous role in string theory: the same calculation in that context is the reason the bosonic string lives 
in 26 = 24 + 2 spacetime dimensions. (The +2 comes from the fact the string itself is extended in one 
space and one time dimension). You will need to attend next term’s “String Theory” course to see 
what on earth this has to do with the Casimir force. 
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We interpret the state | p) as the momentum eigenstate of a single particle of mass m. 
To stress this, from now on we’ll write Ep everywhere instead of c op. Let’s check this 
particle interpretation by studying the other quantum numbers of | p). We may take 
the classical total momentum P given in (1.46) and turn it into an operator. After 
normal ordering, it becomes 

P = - J d 3 x 7T V0 = f 4A p a ^ af (2.45) 

Acting on our state | p) with P, we learn that it is indeed an eigenstate, 

P\p) = p \p) (2.46) 

telling us that the state | p) has momentum p. Another property of | p) that we can 
study is its angular momentum. Once again, we may take the classical expression for 
the total angular momentum of the held (1.55) and turn it into an operator, 

(2.47) 

It’s not hard to show that acting on the one-particle state with zero momentum, 
J l \p = 0) = 0, which we interpret as telling us that the particle carries no internal 
angular momentum. In other words, quantizing a scalar held gives rise to a spin 0 
particle. 

Multi-Particle States, Bosonic Statistics and Fock Space 

We can create multi-particle states by acting multiple times with ad’s. We interpret 
the state in which n ad’s act on the vacuum as an n-particle state, 

\Pi,---,Pn)= 4i • • • a k 1°) ( 2 - 48 ) 

Because all the ad’s commute among themselves, the state is symmetric under exchange 
of any two particles. For example, 


\p,Q) = \q,P) (2.49) 

This means that the particles are bosons. 

The full Hilbert space of our theory is spanned by acting on the vacuum with all 
possible combinations of ad’s, 

lo), 4-1°). 44i°> ■ 4°H-io}(2.50) 
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This space is known as a Fock space. The Fock space is simply the sum of the n-particle 
Hilbert spaces, for all n > 0. There is a useful operator which counts the number of 
particles in a given state in the Fock space. It is called the number operator N 

N = f (^3 4 a P ( 2 - 51 ) 

and satisfies N \pi,... ,p n ) = n\pi,... ,p n ). The number operator commutes with the 
Hamiltonian, [N, H] = 0, ensuring that particle number is conserved. This means that 
we can place ourselves in the n-particlc sector, and stay there. This is a property of 
free theories, but will no longer be true when we consider interactions: interactions 
create and destroy particles, taking us between the different sectors in the Fock space. 

Operator Valued Distributions 

Although we’re referring to the states | p) as “particles”, they’re not localized in space 
in any way — they are momentum eigenstates. Recall that in quantum mechanics the 
position and momentum eigenstates are not good elements of the Hilbert space since 
they are not normalizable (they normalize to delta-functions). Similarly, in quantum 
field theory neither the operators 4>(x), nor a - are good operators acting on the Fock 
space. This is because they don’t produce normalizable states. For example, 

(0|a^at|0) — {p\p) — (27t) 3 «5(0) and (O|0(T) 4>(x) |0) = (x\x) = 5(0) (2.52) 

They are operator valued distributions, rather than functions. This means that al¬ 
though (j>{x) has a well defined vacuum expectation value, (0| <j)(x) |0) = 0, the fluc¬ 
tuations of the operator at a fixed point are infinite, {0\4>(x)4>(x) |0) = oo. We can 
construct well defined operators by smearing these distributions over space. For exam¬ 
ple, we can create a wavepacket 

J-0y> e ~‘ f! v(p)\P) (2.53) 

which is partially localized in both position and momentum space. (A typical state 
might be described by the Gaussian <p(p) — exp(— p 2 /2m 2 )). 

2.4.1 Relativistic Normalization 

We have defined the vacuum |0) which we normalize as (0| 0) = 1. The one-particle 
states | p) = at|0) then satisfy 

(p\q) = (2n) 3 5 {3 \p - q) (2.54) 
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But is this Lorentz invariant? It’s not obvious because we only have 3-vectors. What 
conld go wrong? Suppose we have a Lorentz transformation 

p* 0 p'Y = A \p v (2.55) 

such that the 3-vector transforms as p —» p'. In the quantum theory, it would be 
preferable if the two states are related by a unitary transformation, 

Ip) Ip) = u( A) | p) (2.56) 

This would mean that the normalizations of | p) and |p') are the same whenever p and 
p' are related by a Lorentz transformation. But we haven’t been at all careful with 
normalizations. In general, we could get 

|p) A(p,p') |p') ( 2 . 57 ) 


for some unknown function A(p, p'). How do we figure this out? The trick is to look at 
an object which we know is Lorentz invariant. One such object is the identity operator 
on one-particle states (which is really the projection operator onto one-particle states). 
With the normalization (2.54) we know this is given by 


1 



Ip) (pI 


(2.58) 


This operator is Lorentz invariant, but it consists of two terms: the measure f d 3 p and 
the projector |p)(p|. Are these individually Lorentz invariant? In fact the answer is no. 


Claim The Lorentz invariant measure is, 


/ 


d 3 p 
2 E~p 


(2.59) 


Proof: f d^p is obviously Lorentz invariant. And the relativistic dispersion relation 
for a massive particle, 

p M p M = m 2 =>- Pq = = p 2 + m 2 (2.60) 

is also Lorentz invariant. Solving for po, there are two branches of solutions: po = ±Ep. 
But the choice of branch is another Lorentz invariant concept. So piecing everything 
together, the following combination must be Lorentz invariant, 

J d 4 p 5(po — p 2 — rn 2 ) 

which completes the proof. □ 


p o>0 


d 3 p 

2p 0 


(2.61) 


PO Erf 
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From this result we can figure out everything else. For example, the Lorentz invariant 
5-function for 3-vectors is 

2 E f 6^ 3 \p-q} (2.62) 

which follows because 

/^2Sp'5 (3, (p-9T = 1 (2.63) 

So finally we learn that the relativistically normalized momentum states are given by 

Ip) = Ip) = V^4 1°> ( 2 - 64 ) 

Notice that our notation is rather subtle: the relativistically normalized momentum 
state |p) differs from \p) by the factor ^2Ep. These states now satisfy 

(p\q) = (2t r) 3 2 E^S (3 \p- q) 

Finally, we can rewrite the identity on one-particle states as 

1 = /(S ^> ><p| 

Some texts also define relativistically normalized creation operators by a){p) = s /2EpaL. 
We won’t make use of this notation here. 


(2.65) 

( 2 . 66 ) 


2.5 Complex Scalar Fields 

Consider a complex scalar field ip(x) with Lagrangian 

C = 8 V d»ip - M 2 W (2.67) 

Notice that, in contrast to the Lagrangian (1.7) for a real scalar field, there is no factor 
of 1/2 in front of the Lagrangian for a complex scalar field. If we write ip in terms 
of real scalar fields by ip = (0i + icp 2 )/x/2, we get the factor of 1/2 coming from the 
l/\/2’s. The equations of motion are 

d^d 11 ip + M 2 ip = 0 

d„d»ip* + M 2 ip* = 0 (2.68) 


where the second equation is the complex conjugate of the first. We expand the complex 
field operator as a sum of plane waves as 

d 3 p 1 


ip = 
ip^ = 


(2tt)3 v /2 Ep 
d 3 p 1 

VP y/2E p 


t e -ip-x 


b f e +ip ' x + cte 


e~ ip ' g + c 0 e +ip - s 


(2.69) 
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Since the classical field 0 is not real, the corresponding quantum field if; is not hermitian. 
This is the reason that we have different operators b and c 1 appearing in the positive 
and negative frequency parts. The classical field momentum is 7r = dC /<90 = 0*. We 
also turn this into a quantum operator field which we write as, 



The commutation relations between fields and momenta are given by 

[0(f), tt(^)] = iS {3) (x-y) and [0(f), n\y)] = 0 (2.71) 

together with others related by complex conjugation, as well as the usual [0(f), 0(0)] = 
[0(f), 0' (0)] = 0, etc. One can easily check that these field commutation relations are 
equivalent to the commutation relations for the operators bp and cp, 

[V,$ = (2ir) 3 i <3) (?-9) 

[c p -,ct] = (2 nfS<-V(p-q) (2.72) 

and 

[bp, bp] = [cp, cp] = [bp, cp] = [bp, ct] = 0 (2.73) 

In summary, quantizing a complex scalar field gives rise to two creation operators, 
and ct. These have the interpretation of creating two types of particle, both of mass M 
and both spin zero. They are interpreted as particles and anti-particles. In contrast, 
for a real scalar field there is only a single type of particle: for a real scalar field, the 
particle is its own antiparticle. 


Recall that the theory (2.67) has a classical conserved charge 

Q=i S Sx w-**«*) 

After normal ordering, this becomes the quantum operator 

d 3 p 


Q = 


(2 fk 


(,c'cp-blbp) = N c -N b 


(2.74) 


(2.75) 


so Q counts the number of anti-particles (created by c') minus the number of particles 
(created by 0). We have [H,Q] = 0, ensuring that Q is conserved quantity in the 
quantum theory. Of course, in our free field theory this isn’t such a big deal because 
both N c and 7V& are separately conserved. However, we’ll soon see that in interacting 
theories Q survives as a conserved quantity, while N c and N b individually do not. 
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2.6 The Heisenberg Picture 

Although we started with a Lorentz invariant Lagrangian, we slowly butchered it as we 
quantized, introducing a preferred time coordinate t. It’s not at all obvious that the 
theory is still Lorentz invariant after quantization. For example, the operators 0(x) 
depend on space, but not on time. Meanwhile, the one-particle states evolve in time 
by Schrodinger’s equation, 

= H \iKt)) =* \p(t)) = e-‘ % ‘ \P) (2.76) 

Things start to look better in the ffeisenberg picture where time dependence is assigned 
to the operators O, 


O h = e iHt O s e- iHt (2.77) 

so that 

d ^ = i[H,0„} (2.78) 

where the subscripts S and H tell us whether the operator is in the Schrodinger or 
ffeisenberg picture. In field theory, we drop these subscripts and we will denote the 
picture by specifying whether the fields depend on space 0(5?) (the Schrodinger picture) 
or spacetime 0(5?, t) = (p(x) (the Heisenberg picture). 

The operators in the two pictures agree at a fixed time, say, t — 0. The commutation 
relations (2.2) become equal time commutation relations in the Heisenberg picture, 

[(j)(x,t),(f)(y,t)] = [ir(x,t),Tr{y,t)] =0 

[(f)(x,t),7r(y,t)] =iS^\x-y) (2.79) 


Now that the operator 0(x) = 0(x, t) depends on time, we can start to study how it 
evolves. For example, we have 


0 = i[H, 0] = -[ j d 3 y n(y) 2 + V00/) 2 + m 2 (j){y) 2 , 0(x)] 
= i d 3 y 7 r(y) (—i) S^ 3 \y — x) = ti(x) 


(2.80) 


Meanwhile, the equation of motion for n reads, 


7T = i[H, 7r] = - [ / d 3 y 7T0/) 2 + V0(|/) 2 + m 2 (j){y) 2 , 7r(x)] 
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= \ f d 3 y (V v [(p(y),7T(x)})V(j)(y)+ V(p{y)V y [(p{y),7T{x)] 
+2im 2 (f)(y ) 5^ 3 \x — y) 

= ~ (yf d 3 y (\/yS {3 \x-y)) V y <f>(y)^ -m 2 ^) 

= V 2 (p — m 2 (p (2.81) 


where we’ve included the subscript y on V y when there may be some confusion about 
which argument the derivative is acting on. To reach the last line, we’ve simply inte¬ 
grated by parts. Putting (2.80) and (2.81) together we find that the field operator (p 
satisfies the Klcin-Gordon equation 

dydy + m 2 (j) = 0 (2.82) 


Things are beginning to look more relativistic. We can write the Fourier expansion of 
c p(x ) by using the definition (2.77) and noting, 


iHt - " lHt — lEpt and e lHt ale lHt = e +lEpt al 


e""“ a^e lHt = e 


(2.83) 


which follows from the commutation relations [H,ap\ = —Ep a-and [H,al] = +Epal. 
This then gives, 


^ i) = /(05y=(^* + 4 e ' 


which looks very similar to the previous expansion (2.18) 
except that the exponent is now written in terms of 4- 
vectors, p ■ x = Ept — p ■ x. (Note also that a sign has 
flipped in the exponent due to our Minkowski metric con¬ 
traction). ft’s simple to check that (2.84) indeed satisfies 
the Klein-Gordon equation (2.82). 


t e +ip-~ 


(2.84) 



2.6.1 Causality 

We’re approaching something Lorentz invariant in the 
Heisenberg picture, where 4>(x) now satisfies the Klein- 

Gordon equation. But there’s still a hint of non-Lorentz invariance because (p an d vr 
satisfy equal time commutation relations, 


(x,t),ir(y,t)] = i5 {3) (x-y ) 


(2.85) 
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But what about arbitrary spacetime separations? In particular, for our theory to be 
causal, we must require that all spacelike separated operators commute, 

[Oi(x),0 2 {y)\= 0 V (x-y) 2 < 0 (2.86) 

This ensures that a measurement at x cannot affect a measurement at y when x and y 
are not causally connected. Does our theory satisfy this crucial property? Let’s define 

A (x-y) = [<l>(x),<j>(y)\ (2.87) 

The objects on the right-hand side of this expression are operators. However, it’s easy 
to check by direct substitution that the left-hand side is simply a c-number function 
with the integral expression 

A(x ~ v) = j Wji Wp ( e ~‘ r<X ~ V> - e"" 1 ’" 5 ’) (2.88) 

What do we know about this function? 

• It’s Lorentz invariant, thanks to the appearance of the Lorentz invariant measure 
f d 3 p/2Ep that we introduced in (2.59). 

• It doesn’t vanish for timelike separation. For example, taking x — y = (t, 0, 0, 0) 
gives [4>(x,0),(j)(x,t)] ~ e~ imt — e +imt . 

• It vanishes for space-like separations. This follows by noting that A(x — y) — 0 
at equal times for all (x — y) 2 = — (x — y) 2 < 0, which we can see explicitly by 
writing 

Mf, t), t)\ = \f 7^3 *- (2.89) 

J (27T) d 2 + m 2 

and noticing that we can flip the sign of p in the last exponent as it is an inte¬ 
gration variable. But since A(x — y) Lorentz invariant, it can only depend on 
(x — y) 2 and must therefore vanish for all (x — y) 2 < 0. 

We therefore learn that our theory is indeed causal with commutators vanishing 
outside the lightcone. This property will continue to hold in interacting theories; indeed, 
it is usually given as one of the axioms of local quantum held theories. I should mention 
however that the fact that [<p{x), 4>(y)] is a c-number function, rather than an operator, 
is a property of free fields only. 
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2.7 Propagators 


We could ask a different question to probe the causal structure of the theory. Prepare 
a particle at spacetime point y. What is the amplitude to find it at point x? We can 
calculate this: 


( 0 | <j>{x)<j>{y) | 0 ) 


d 3 pd 3 p' 1 
(27r) 6 E^, 


a p a p‘ 


d3 P 1 c -ip-(x-v) 

(2 t r) 3 2E P 


D{x — 


0) e 


-ip-x-\-ip'-y 


y ) 


(2.90) 


The function D{x — y ) is called the propagator. For spacelike separations, (x — y) 2 < 0, 
one can show that D{x — y) decays like 


D(x-y)~ e~ mM 


(2.91) 


So it decays exponentially quickly outside the lightcone but, nonetheless, is non-vanishing! 
The quantum held appears to leak out of the lightcone. Yet we’ve just seen that space¬ 
like measurements commute and the theory is causal. How do we reconcile these two 
facts? We can rewrite the calculation (2.89) as 

[<f>(x),<f>(y)] = D(x-y)~ D(y~x) = 0 if (x - y) 2 < 0 (2.92) 

There are words you can drape around this calculation. When (x — y) 2 < 0, there 
is no Lorentz invariant way to order events. If a particle can travel in a spacelike 
direction from x —> y, it can just as easily travel from y —» x. In any measurement, the 
amplitudes for these two events cancel. 

With a complex scalar held, it is more interesting. We can look at the equation 
[i/j( x), (y)\ = 0 outside the lightcone. The interpretation now is that the amplitude 
for the particle to propagate from x —> y cancels the amplitude for the antiparticle 
to travel from y —y x. In fact, this interpretation is also there for a real scalar held 
because the particle is its own antiparticle. 


2.7.1 The Feynman Propagator 

As we will see shortly, one of the most important quantities in interacting held theory 
is the Feynman propagator , 


A F (x - y) 


(0| T(j)(x)(f>(y) |0) = 


D(x - y) 
D(y - x) 


x° > y° 

yO > x o 


(2.93) 
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where T stands for time ordering, placing all operators evaluated at later times to the 
left so, 


T<j>{x)<j>{y) 


4>(x)(j)(y) x° > y° 

4>(y)(f>(x) y° > x° 


(2.94) 


Claim: There is a useful way of writing the Feynman propagator in terms of a 4- 
momentum integral. 

A f(x - y) = f - j -5 — - e-*<— '> (2.95) 

J (27r) 4 p 2 — m 2 

Notice that this is the hrst time in this course that we’ve integrated over 4-momentum. 
Until now, we integrated only over 3-momentum, with p° fixed by the mass-shell con¬ 
dition to be p° = Ep. In the expression (2.95) for A p, we have no such condition on 
p°. However, as it stands this integral is ill-defined because, for each value of p, the 
denominator p 2 — m 2 = (p 0 ) 2 — p 2 — m 2 produces a pole when p° = EEp = ± \Jp 2 + m 2 . 
We need a prescription for avoiding these singularities in the po integral. To get the 
Feynman propagator, we must choose the contour to be 



Figure 5: The contour for the Feynman propagator. 


Proof: 


p 2 m 2 ( p°) 2 -E 2 (p« - E p )(pP + E p ) ' 

so the residue of the pole at p° = ±Ep is ±1/2 Ep. When x° > y°, we close the contour 
in the lower half plane, where p° — > —ioo, ensuring that the integrand vanishes since 
e -ip°(x°-y°) q xhe integral over p° then picks up the residue at p° = +Ep which 
is —2ni/2Ep where the minus sign arises because we took a clockwise contour. Hence 
when x° > y° we have 


- y) 


d3 P ~ 2m j p -iE^x°-y°)+ipix-y) 

(27r) 4 2 Ep 
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D(x - y) 


(2.97) 


d 3 p 1 
(2vr) 3 2Ep 


,-ip-(x-y) _ 


which is indeed the Feynman propagator for x° > y°. In contrast, when y° > x°, we 
close the contour in an anti-clockwise direction in the upper half plane to get, 


A f (x - y) 


d3 P 2711 j +iE^x°-y°)+ipix-y) 

(2vr) 4 (-2 Ep) 

d P 1 p -iEp(y°-X°)~ip-(y-x) 
(2vr)3 2 Ep 

^ P _ L_ e -ip-(y-x) — ]j(y _ x \ 

(2vr) 3 2 E p e } 


(2.98) 


where to go to from the second line to the third, we have flipped the sign of p which 
is valid since we integrate over d 3 p and all other quantities depend only on p 2 . Once 
again we reproduce the Feynman propagator. □ 


Instead of specifying the contour, it is standard to write 
the Feynman propagator as 


A F (x - y) 


d i p i e -ip-(x-y) 
(27t) 4 p 2 — m 2 + ie 


(2.99) 


with e > 0, and infinitesimal. This has the effect of shifting 
the poles slightly off the real axis, so the integral along the 
real p° axis is equivalent to the contour shown in Figure 5. 
This way of writing the propagator is, for obvious reasons, 
called the “ie prescription”. 


Im(p°) 


1 


+ ie 


> 


Re(p°) 
1 - i£ 


Figure 6: 


2.7.2 Green’s Functions 


There is another avatar of the propagator: it is a Green’s function for the Klcin-Gordon 
operator. If we stay away from the singularities, we have 


(dl — V 2 + m 2 ) A f (x — y) 


d*p 

(27t) 4 p 2 


m “ 


(-p 2 + m 2 ) e - ip ' {x ~ y) 


d P „-ip-(x-y) 


J (2vr) 4 

—i 8^ 4 \x — y) 


( 2 . 100 ) 


Note that we didn’t make use of the contour anywhere in this derivation. For some 
purposes it is also useful to pick other contours which also give rise to Green’s functions. 
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Figure 7: The retarded contour Figure 8: The advanced contour 


For example, the retarded Green’s function A r(x — y ) is defined by the contour shown 
in Figure 7 which has the property 


A R {x -y) = 


D{x — y) — D(y — x) x° > y° 


0 


yO > x 0 


( 2 . 101 ) 


The retarded Green’s function is useful in classical held theory if we know the initial 
value of some held conhguration and want to figure out what it evolves into in the 
presence of a source, meaning that we want to know the solution to the inhomogeneous 
Klcin-Gordon equation, 


+ m 2 (f> = J{x ) (2.102) 

for some hxed background function J(x). Similarly, one can define the advanced Green’s 
function A^(x — y ) which vanishes when y° < x°, which is useful if we know the end 
point of a held conhguration and want to figure out where it came from. Given that 
next term’s course is called “Advanced Quantum Field Theory”, there is an obvious 
name for the current course. But it got shot down in the staff meeting. In the quantum 
theory, we will see that the Feynman Green’s function is most relevant. 

2.8 Non-Relativistic Fields 

Let’s return to our classical complex scalar field obeying the Klcin-Gordon equation. 
We’ll decompose the held as 


^(Xjt) = e imt ip(x,t) 


(2.103) 


Then the KG-equation reads 


<9 t 2 ^ — V 2 '!/’ + rn 2r ip = e 


—imt 


"0 — 2 im^\) — V 2 ^ 


= 0 


(2.104) 


with the m 2 term cancelled by the time derivatives. The non-relativistic limit of a 
particle is \p\ <C m. Let’s look at what this does to our held. After a Fourier transform, 
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this is equivalent to saying that |^)| <C m|^|. In this limit, we drop the term with two 
time derivatives and the KG equation becomes, 

’w = (2 - i05) 

This looks very similar to the Schrodinger equation for a non-relativistic free particle of 
mass m. Except it doesn’t have any probability interpretation — it’s simply a classical 
field evolving through an equation that’s first order in time derivatives. 


We wrote down a Lagrangian in section 1.1.2 which gives rise to field equations 
which are first order in time derivatives. In fact, we can derive this from the relativistic 
Lagrangian for a scalar field by again taking the limit (9*0 -C m0. After losing the 
tilde, so ip —y 0, the non-relativistic Lagrangian becomes 

£ =+iip*ip -V0* V0 (2.106) 

2m 

where we’ve divided by 1/2 m. This Lagrangian has a conserved current arising from 
the internal symmetry 0 —>■ e*“0. The current has time and space components 


f= ( 


(2.107) 


To move to the Hamiltonian formalism we compute the momentum 


7T = 



(2.108) 


This means that the momentum conjugate to 0 is i0*. The momentum does not depend 
on time derivatives at all! This looks a little disconcerting but it’s fully consistent for a 
theory which is first order in time derivatives. In order to determine the full trajectory 
of the field, we need only specify 0 and 0* at time t = 0: no time derivatives on the 
initial slice are required. 


Since the Lagrangian already contains a “pg” term (instead of the more familiar \pq 
term), the time derivatives drop out when we compute the Hamiltonian. We get, 

n = — V0*V0 (2.109) 

2m 

To quantize we impose (in the Schrodinger picture) the canonical commutation relations 


[0(f), 0(0)] = [0t(f),0t(0)] = 0 

[0(f),000)] = <5 (3 ) (f-p) (2.110) 
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We may expand ^(x) as a Fourier transform 



( 2 . 111 ) 


where the commutation relations (2.110) require 


[«P’4] = (2vr ) 3 5 {3 \p-q) 


( 2 . 112 ) 


The vacuum satishes a-|0) = 0, and the excitations are at ... at j0). The one-particle 
states have energy 



(2.113) 


which is the non-relativistic dispersion relation. We conclude that quantizing the first 
order Lagrangian (2.106) gives rise to non-relativistic particles of mass m. Some com¬ 
ments: 

• We have a complex held but only a single type of particle. The anti-particle is 
not in the spectrum. The existence of anti-particles is a consequence of relativity. 

• A related fact is that the conserved charge Q = f d 3 x : : is the particle 

number. This remains conserved even if we include interactions in the Lagrangian 
of the form 


AC = y(v>*v>) 


(2.114) 


So in non-relativistic theories, particle number is conserved, ft is only with rela¬ 
tivity, and the appearance of anti-particles, that particle number can change. 

• There is no non-relativistic limit of a real scalar held. In the relativistic theory, 
the particles are their own anti-particles, and there can be no way to construct a 
multi-particle theory that conserves particle number. 

2.8.1 Recovering Quantum Mechanics 

In quantum mechanics, we talk about the position and momentum operators X and 
P. In quantum held theory, position is relegated to a label. How do we get back to 
quantum mechanics? We already have the operator for the total momentum of the held 



(2.115) 
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which, on one-particle states, gives P\p) = p\p). It’s also easy to construct the position 
operator. Let’s work in the non-relativistic limit. Then the operator 

d 3 p 


^d(f) = 


ate-** 


(2.116) 


, (2vr) 3 p 

creates a particle with 5-function localization at x. We write |f) = i]fi(x) |0). A natural 
position operator is then 

V = (2.117) 

so that X |f) — x\x). 

Let’s now construct a state \ip) by taking superpositions of one-particle states |f), 

| (p) = J d 3 x p(x) \x) (2.118) 

The function p(x) is what we would usually call the Schrodinger wavefunction (in the 
position representation). Let’s make sure that it indeed satisfies all the right properties. 
Firstly, it’s clear that acting with the position operator X has the right action of <p(x), 

X' \y) = J #: 

but what about the momentum operator P 

P i \<p) = j d 3 x 

which tells us that P l acts as the familiar derivative on wavefunctions \p). To see that 
this is the case, we write 


x 1 p(x) 

f) 

(2.119) 

We will 

now show that 


-a 

f) 

(2.120) 


p‘\f) = 


d j0^p'a t l p f <p(x)ili t (x) | 0 ) 

"Vat-VH) 10) 


j (2ir) 3 r (2.121) 

where we’ve used the relationship [a-, ^)'{x)\ = e~ ip ' x which can be easily checked. 
Proceeding with our calculation, we have 


P l W) = 


"" 4 ( 4 ^ 1 ^ 10 ) 


(27t) 3 p \ dx 

r 

-e 


d 3 xd 3 p r _ ifS f_,dp \ j 
(2tt) 


-i^ at |0) 

dx'J p 1 ' 


. dp 


= I d 3 x [ ) |f) 


( 2 . 122 ) 
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which confirms (2.120). So we learn that when acting on one-particle states, the oper¬ 
ators X and P act as position and momentum operators in quantum mechanics, with 
[X*,Pi] | cp) = i5 l i | <p). But what about dynamics? How does the wavefunction (p(x,t) 
change in time? The Hamiltonian (2.109) can be rewritten as 


H 





d 3 p p 2 | 
(27t) 3 2 m p ^ 


(2.123) 


so we find that 




(2.124) 


But this is the same equation obeyed by the original held (2.105)! Except this time, it 
really is the Schrodinger equation, complete with the usual probabilistic interpretation 
for the wavefunction ip. Note in particular that the conserved charge arising from the 
Noether current (2.107) is Q = f d 3 x \p(x)\ 2 which is the total probability. 


Historically, the fact that the equation for the classical held (2.105) and the one- 
particle wavefunction (2.124) coincide caused some confusion, ft was thought that 
perhaps we are quantizing the wavefunction itself and the resulting name “second quan¬ 
tization” is still sometimes used today to mean quantum held theory. It’s important to 
stress that, despite the name, we’re not quantizing anything twice! We simply quantize 
a classical held once. Nonetheless, in practice it’s useful to know that if we treat the 
one-particle Schrodinger equation as the equation for a quantum held then it will give 
the correct generalization to multi-particle states. 


Interactions 

Often in quantum mechanics, we’re interested in particles moving in some hxed back¬ 
ground potential V(x). This can be easily incorporated into held theory by working 
with a Lagrangian with explicit x dependence, 

C = (2.125) 

2m 

Note that this Lagrangian doesn’t respect translational symmetry and we won’t have 
the associated energy-momentum tensor. While such Lagrangians are useful in con¬ 
densed matter physics, we rarely (or never) come across them in high-energy physics, 
where all equations obey translational (and Lorentz) invariance. 
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One can also consider interactions between particles. Obviously these are only impor¬ 
tant for n particle states with n > 2. We therefore expect them to arise from additions 
to the Lagrangian of the form 

A£ = ip*(x) ip(x) (2.126) 

which, in the quantum theory, is an operator which destroys two particles before creat¬ 
ing two new ones. Such terms in the Lagrangian will indeed lead to inter-particle forces, 
both in the non-relativistic and relativistic setting. In the next section we explore these 
types of interaction in detail for relativistic theories. 
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3. Interacting Fields 


The free field theories that we’ve discussed so far are very special: we can determine 
their spectrum, but nothing interesting then happens. They have particle excitations, 
but these particles don’t interact with each other. 

Here we’ll start to examine more complicated theories that include interaction terms. 
These will take the form of higher order terms in the Lagrangian. We’ll start by asking 
what kind of small perturbations we can add to the theory. For example, consider the 
Lagrangian for a real scalar field, 

c = - imV I 3 - 1 ) 

n> 3 n ' 

The coefficients \ n are called coupling constants. What restrictions do we have on X n 
to ensure that the additional terms are small perturbations? You might think that we 
need simply make u X n <C 1”. But this isn’t quite right. To see why this is the case, let’s 
do some dimensional analysis. Firstly, note that the action has dimensions of angular 
momentum or, equivalently, the same dimensions as h. Since we’ve set h = 1, using 
the convention described in the introduction, we have [S'] = 0. With S = f d 4 x C, and 
[i d 4 x] = —4, the Lagrangian density must therefore have 

[£] = 4 (3.2) 

What does this mean for the Lagrangian (3.1)? Since [d^] = 1, we can read off the 
mass dimensions of all the factors to find, 

[<fi] = 1 , M = 1 , [A n ] =4-n (3.3) 

So now we see why we can’t simply say we need A n -C 1, because this statement only 
makes sense for dimensionless quantities. The various terms, parameterized by A n , fall 
into three different categories 

• [A 3 ] = 1: For this term, the dimensionless parameter is A 3 /E, where E has 
dimensions of mass. Typically in quantum field theory, E is the energy scale of 
the process of interest. This means that X^(j) 3 /3\ is a small perturbation at high 
energies E A 3 , but a large perturbation at low energies f? C A 3 . Terms that 
we add to the Lagrangian with this behavior are called relevant because they’re 
most relevant at low energies (which, after all, is where most of the physics we see 
lies). In a relativistic theory, E > m, so we can always make this perturbation 
small by taking A 3 <C m. 
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• [A 4 ] = 0: this term is small if A 4 1. Such perturbations are called marginal. 

• [A n ] < 0 for n > 5: The dimensionless parameter is (A n E n ~ 4 ), which is small at 
low-energies and large at high energies. Such perturbations are called irrelevant. 

As you’ll see later, it is typically impossible to avoid high energy processes in quantum 
held theory. (We’ve already seen a glimpse of this in computing the vacuum energy). 
This means that we might expect problems with irrelevant operators. Indeed, these lead 
to “non-renormalizable” field theories in which one cannot make sense of the infinities 
at arbitrarily high energies. This doesn’t necessarily mean that the theory is useless; 
just that it is incomplete at some energy scale. 

Let me note however that the naive assignment of relevant, marginal and irrelevant 
is not always fixed in stone: quantum corrections can sometimes change the character 
of an operator. 


An Important Aside: Why QFT is Simple 

Typically in a quantum held theory, only the relevant and marginal couplings are 
important. This is basically because, as we’ve seen above, the irrelevant couplings 
become small at low-energies. This is a huge help: of the infinite number of interaction 
terms that we could write down, only a handful are actually needed (just two in the 
case of the real scalar field described above). 


Let’s look at this a little more. Suppose that we some day discover the true su- 
perduper “theory of everything unimportant” that describes the world at very high 
energy scales, say the GUT scale, or the Planck scale. Whatever this scale is, let’s call 
it A. It is an energy scale, so [A] = 1. Now we want to understand the laws of physics 
down at our puny energy scale fcA. Let’s further suppose that down at the energy 
scale E, the laws of physics are described by a real scalar field. (They’re not of course: 
they’re described by non-Abelian gauge fields and fermions, but the same argument 
applies in that case so bear with me). This scalar field will have some complicated 
interaction terms (3.1), where the precise form is dictated by all the stuff that’s going 
on in the high energy superduper theory. What are these interactions? Well, we could 
write our dimensionful coupling constants X n in terms of dimensionless couplings g n , 
multiplied by a suitable power of the relevant scale A, 


A 


n 


9n 

An - 4 


(3.4) 


The exact values of dimensionless couplings g n depend on the details of the high-energy 
superduper theory, but typically one expects them to be of order 1: g n ~ 0(1). This 
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means that for experiments at small energies E •C A, the interaction terms of the 
form (f) 11 with n > 4 will be suppressed by powers of (E/A) n ~ 4 . This is usually a 
suppression by many orders of magnitude, (e.g for the energies E explored at the 
LHC, E/M p i ~ 1CT 16 ). It is this simple argument, based on dimensional analysis, that 
ensures that we need only focus on the first few terms in the interaction: those which 
are relevant and marginal. It also means that if we only have access to low-energy 
experiments (which we do!), it’s going to be very difficult to figure out the high energy 
theory (which it is!), because its effects are highly diluted except for the relevant and 
marginal interactions. The discussion given above is a poor man’s version of the ideas 
of effective field theory and Wilson’s renormalization group , about which you can learn 
more in the “Statistical Field Theory” course. 

Examples of Weakly Coupled Theories 

In this course we’ll study only weakly coupled field theories i.e. ones that can truly be 
considered as small perturbations of the free field theory at all energies. In this section, 
we’ll look at two types of interactions 

1) 0 4 theory: 

£ = - ^0 4 (3.5) 

with A 1. We can get a hint for what the effects of this extra term will be. Expanding 
out 0 4 in terms of a- and at, we see a sum of interactions that look like 

atatatat and atatata^ etc. (3.6) 

p p p p p p p P V/ 

These will create and destroy particles. This suggests that the 0 4 Lagrangian describes 
a theory in which particle number is not conserved. Indeed, we could check that the 
number operator N now satisfies [H, N] f 0. 

2) Scalar Yukawa Theory 

C = — ilf 2 0*0 — ^m 2 0 2 — (3.7) 

with g < M, m. This theory couples a complex scalar -0 to a real scalar 0. While 
the individual particle numbers of 0 and 0 are no longer conserved, we do still have 
a symmetry rotating the phase of 0, ensuring the existence of the charge Q defined 
in (2.75) such that [Q,H] = 0. This means that the number of 0 particles minus the 
number of anti-particles is conserved. It is common practice to denote the anti¬ 
particle as 0. 
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The scalar Yukawa theory has a slightly worrying aspect: the potential has a stable 
local minimum at 0 = ip — 0, but is unbounded below for large enough —g<j). This 
means we shouldn’t try to push this theory too far. 


A Comment on Strongly Coupled Field Theories 

In this course we restrict attention to weakly coupled held theories where we can use 
perturbative techniques. The study of strongly coupled held theories is much more 
difficult, and one of the major research areas in theoretical physics. For example, some 
of the amazing things that can happen include 

• Charge Fractionalization: Although electrons have electric charge 1, under 
the right conditions the elementary excitations in a solid have fractional charge 
1/A (where A e 2Z + 1). For example, this occurs in the fractional quantum 
Hall effect. 

• Confinement: The elementary excitations of quantum chromodynamics (QCD) 
are quarks. But they never appear on their own, only in groups of three (in a 
baryon) or with an anti-quark (in a meson). They are confined. 

• Emergent Space: There are held theories in four dimensions which at strong 
coupling become quantum gravity theories in ten dimensions! The strong cou¬ 
pling effects cause the excitations to act as if they’re gravitons moving in higher 
dimensions. This is quite extraordinary and still poorly understood. It’s called 
the AdS/CFT correspondence. 


3.1 The Interaction Picture 


There’s a useful viewpoint in quantum mechanics to describe situations where we have 
small perturbations to a well-understood Hamiltonian. Let’s return to the familiar 
ground of quantum mechanics with a finite number of degrees of freedom for a moment. 
In the Schrodinger picture, the states evolve as 


•4 0>s 

dtt 


H\^) s 


while the operators Os are independent of time. 


(3.8) 


In contrast, in the Heisenberg picture the states are fixed and the operators change 
in time 


0 H (t) = e iHt O s e- iHt 

Mh = e iHt Ms ( 3 - 9 ) 
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The interaction picture is a hybrid of the two. We split the Hamiltonian np as 

H = Hq + H mt (3.10) 

The time dependence of operators is governed by Ho, while the time dependence of 
states is governed by H int . Although the split into H 0 and H[ nt is arbitrary, it’s useful 
when Ho is soluble (for example, when H 0 is the Hamiltonian for a free field theory). 
The states and operators in the interaction picture will be denoted by a subscript / 
and are given by, 


\m), = e‘ H °‘ W))s 

£>,(t) = e iHot Ose~ iHot (3.11) 

This last equation also applies to H int , which is time dependent. The interaction 
Hamiltonian in the interaction picture is, 

Ht = ( Hint)i = e iHot (Hint)s e- iHot (3.12) 

The Schrodinger equation for states in the interaction picture can be derived starting 
from the Schrodinger picture 

i ^dt A = Hs ^ s ^ l Jt ^ € ~ iHOt ^ = ( H ° + H ' mt ^ S e ~ iHOt 

=► = e iHot (H mt )se- iHot Wj (3.13) 

So we learn that 

(3.14) 


3.1.1 Dyson’s Formula 

“Well, Birmingham has much the best theoretical physicist to work with, 
Peierls; Bristol has much the best experimental physicist, Powell; Cam¬ 
bridge has some excellent architecture. You can make your choice.” 

Oppenheimer’s advice to Dyson on which university position to accept. 
We want to solve (3.14). Let’s write the solution as 

I = U(t, t 0 ) \^(to)) r (3.15) 
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where U(t, t 0 ) is a unitary time evolution operator such that U(ti,t 2 )U(t 2 , t 3 ) = U(ti,t 3 ) 
and U(t,t ) = 1. Then the interaction picture Schrodinger equation (3.14) requires that 


dIJ 

%-r- 

dt 


H T (t)U 


(3.16) 


If Hi were a function, then we could simply solve this by 

U(t, t 0 ) = exp ^ —i £ Hi(t') dt'^ (3-17) 

But there’s a problem. Our Hamiltonian Hi is an operator, and we have ordering 
issues. Let’s see why this causes trouble. The exponential of an operator is defined in 
terms of the expansion, 

exp i Hi{t') dt'^j —l ~ i £ Hr(t') dt' + ^ ^ J ^ dt'^j + ... (3.18) 


But when we try to differentiate this with respect to t, we find that the quadratic term 
gives us 


1 

2 





Hi{t')dt' 


(3.19) 


Now the second term here looks good, since it will give part of the Hi(t)U that we 
need on the right-hand side of (3.16). But the first term is no good since the Hi(t) 
sits the wrong side of the integral term, and we can’t commute it through because 
[Hi(t'), Hi(t )] ^ 0 when t' ^ t. So what’s the way around this? 


Claim: The solution to (3.16) is given by Dyson’s Formula. (Essentially Erst figured 
out by Dirac, although the compact notation is due to Dyson). 

rt 


U(t,to) = T exp ( —i / Hi(t') dt' 


'to 


(3.20) 


where T stands for time ordering where operators evaluated at later times are placed 
to the left 


T {0 1 {t l )0 2 {t 2 )) = 


(h) 0 2 (t 2 ) t.\ > t 2 

0 2 (t 2 ) 01 (ti) t 2 > ti 

Expanding out the expression (3.20), we now have 


(3.21) 


U(t,t Q ) = 1 -i dt’ Hi[t’) + 


(-0 s 


2 r rt 


'to 


dt ' / dt" Hi(t")Hi(t') 


-J to 


+ / dt' dt" Hi{t')Hi(t") 


'to 


'to 


+ ... 
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Actually these last two terms double up since 



dt" Hj(O#/(0 



dt' 


'to 



dt" Hi(lf)Hi{t n 


'to 


(3.22) 


where the range of integration in the first expression is over t" > t', while in the second 
expression it is t' < t" which is, of course, the same thing. The final expression is the 
same as the second expression by a simple relabelling. This means that we can write 


U{t,t 0 ) = 1 


dt'H^t') + (-i) 2 / dt' / dt" Hi(t')Hi(t") + 


'to 


'to 


'to 


(3.23) 


Proof: The proof of Dyson’s formula is simpler than explaining what all the nota¬ 
tion means! Firstly observe that under the T sign, all operators commute (since their 
order is already fixed by the T sign). Thus 


i— T exp i ^ dt' Hi(t 1 )^ = T 


Hi(t) exp ( —i / dt' Hi (t') 


'to 


= H T (t) T exp ^ -i dt' H^t'^j (3.24) 


since t, being the upper limit of the integral, is the latest time so Hj(t) can be pulled 
out to the left. □ 


Before moving on, I should confess that Dyson’s formula is rather formal. It is 
typically very hard to compute time ordered exponentials in practice. The power of 
the formula comes from the expansion which is valid when Hi is small and is very easily 
computed. 


3.2 A First Look at Scattering 

Let us now apply the interaction picture to held theory, starting with the interaction 
Hamiltonian for our scalar Yukawa theory, 


Hint = 9 


I 


d 3 x 


(3.25) 


Unlike the free theories discussed in Section 2, this interaction doesn’t conserve particle 
number, allowing particles of one type to morph into others. To see why this is, we use 
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the interaction picture and follow the evolution of the state: \^{t)) = U(t,t 0 ) |^(t 0 )), 
where U(t,to) is given by Dyson’s formula (3.20) which is an expansion in powers of 
H mt . But Hi n t contains creation and annihilation operators for each type of particle. 
In particular, 

• (j) ~ a + ah This operator can create or destroy (j) particles. Let’s call them 
mesons. 

• 'ijj ~ b + ck This operator can destroy ijj particles through b, and create anti¬ 
particles through cb Let’s call these particles nucleons. Of course, in reality 
nucleons are spin 1/2 particles, and don’t arise from the quantization of a scalar 
field. But we’ll treat our scalar Yukawa theory as a toy model for nucleons 
interacting with mesons. 

• -0' ~ frt + c: This operator can create nucleons through , and destroy anti¬ 
nucleons through c. 

Importantly, Q = N c — Nb remains conserved in the presence of H- mt . At first order in 
perturbation theory, we find terms in H mi like cWa. This kills a meson, producing a 
nucleon-anti-nucleon pair. It will contribute to meson decay 0 —y 'ip'ip. 

At second order in perturbation theory, we’ll have more complicated terms in ( H int ) 2 , 
for example a)(cba)). This term will give contributions to scattering processes 
00 —> (j) —> 00 . The rest of this section is devoted to computing the quantum ampli¬ 
tudes for these processes to occur. 

To calculate amplitudes we make an important, and slightly dodgy, assumption: 

Initial and final states are eigenstates of the free theory 

This means that we take the initial state | i) at t —> — oo, and the final state j/) at 
t —>• Too, to be eigenstates of the free Hamiltonian H 0 . At some level, this sounds 
plausible: at t —> —oo, the particles in a scattering process are far separated and don’t 
feel the effects of each other. Furthermore, we intuitively expect these states to be 
eigenstates of the individual number operators N, which commute with Hq, but not 
H mt ■ As the particles approach each other, they interact briefly, before departing again, 
each going on its own merry way. The amplitude to go from |z) to |/) is 

lim (/| U(t+,t-) \i) = (f\S\i) (3.26) 

t ±—»±oo 

where the unitary operator S is known as the S-matrix. (S is for scattering). There 
are a number of reasons why the assumption of non-interacting initial and final states 
is shaky: 
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• Obviously we can’t cope with bound states. For example, this formalism can’t 
describe the scattering of an electron and proton which collide, bind, and leave 
as a Hydrogen atom. It’s possible to circumvent this objection since it turns out 
that bound states show up as poles in the S-matrix. 

• More importantly, a single particle, a long way from its neighbors, is never alone 
in held theory. This is true even in classical electrodynamics, where the electron 
sources the electromagnetic held from which it can never escape. In quantum 
electrodynamics (QED), a related fact is that there is a cloud of virtual photons 
surrounding the electron. This line of thought gets us into the issues of renormal¬ 
ization — more on this next term in the “AQFT” course. Nevertheless, motivated 
by this problem, after developing scattering theory using the assumption of non¬ 
interacting asymptotic states, we’ll mention a better way. 

3.2.1 An Example: Meson Decay 

Consider the relativistically normalized initial and hnal states, 

N> = 1°) 

l/) = 1°> < 3 - 27 ) 

The initial state contains a single meson of momentum p; the hnal state contains a 
nucleon-anti-nucleon pair of momentum q\ and g 2 - We may compute the amplitude for 
the decay of a meson to a nucleon-anti-nucleon pair. To leading order in g. it is 

(f\S\i) = -ig (f\ I d 4 x^(x)^(x)(j)(x) \i) (3.28) 

Let’s go slowly. We hrst expand out <fi ~ a + cd using (2.84). (Remember that the (j) in 
this formula is in the interaction picture, which is the same as the Heisenberg picture 
of the free theory). The a piece will turn |f) into something proportional to 10), while 
the cd piece will turn |z) into a two meson state. But the two meson state will have 
zero overlap with (/1, and there’s nothing in the i/j and ^ operators that lie between 
them to change this fact. So we have 

(/I S \i) = —ig (/| I d 4 x^\x)il:(x) J |°) 

= (f\ J d 4 x r ip\x)'ip(x)e~ ip ' x |0) (3.29) 

where, in the second line, we’ve commuted a% past at, picking up a 5( 3 \p — k) delta- 
function which kills the d 3 k integral. We now similarly expand out d ~ b + c' and 
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rjr ~ + c. To get non-zero overlap with (/|, only the /h and (d contribute, for they 

create the nucleon and anti-nucleon from |0). We then have 


(f\S\i) = -ig 



d 4 xd 3 kid 3 k 2 \/EqxE, 


O) 6 

= -ig (2j r) 4 S m (qi +t 2 ~p) 

and so we get our first quantum held theory amplitude. 


92 c 0 -6,ct bl |0) e *( fei + fc 2 -p)-a 
92 91 jfcj fc 2 1 ' 


(3.30) 


Notice that the <5-function puts constraints on the possible decays. In particular, the 
decay only happens at all if m > 2 M. To see this, we may always boost ourselves 
to a reference frame where the meson is stationary, so p = (m, 0,0,0). Then the 
delta function imposes momentum conservation, telling us that q\ = —q 2 and m = 
2 a/M 2 + \q\ 2 . 

Later you will learn how to turn this quantum amplitude into something more phys¬ 
ical, namely the lifetime of the meson. The reason this is a little tricky is that we must 
square the amplitude to get the probability for decay, which means we get the square 
of a ^-function. We’ll explain how to deal with this in Section 3.6 below, and again in 
next term’s “Standard Model” course. 


3.3 Wick’s Theorem 

From Dyson’s formula, we want to compute quantities like (/| T {Hj(x i) ... Hi(x n )} |i), 
where |i) and |/) are eigenstates of the free theory. The ordering of the operators is fixed 
by T, time ordering. However, since the L/j’s contain certain creation and annihilation 
operators, our life will be much simpler if we can start to move all annihilation operators 
to the right where they can start killing things in |i). Recall that this is the definition 
of normal ordering. Wick’s theorem tells us how to go from time ordered products to 
normal ordered products. 

3.3.1 An Example: Recovering the Propagator 

Let’s start simple. Consider a real scalar held which we decompose in the Heisenberg 
picture as 

<f>(x) = 

where 

(p + (x) = j 
<i>~(x) = j 


(x) + (j) (x) 


(3.31) 


d 3 p 1 


a-e 


-ip-x 


(2vr) 3 ^2E p 

d3 P 1 „t 

(2vr) 3 ^2E~ p 


aue 


(3.32) 
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where the ± signs on ejr 1 make little sense, but apparently you have Pauli and Heisenberg 
to blame. (They come about because (p + ~ e~ lEt , which is sometimes called the positive 
frequency piece, while </>“ ~ e +^Et j g ^.j ie ne g a ti ve frequency piece). Then choosing 
x° > y°, we have 

T (j){x)<t>{y) = (t>(x)(j)(y) 

= (0 + (x) + 0"(x))(0 + (?/) + 0"(y)) (3.33) 

= (p + (x)(j) + (y) + (j)~(x)(j) + (y) + <j>~{y) (p + {x) + [0 + (x), <f>~{y)\ + </>~(x) <j>~(y) 

where the last line is normal ordered, and for our troubles we have picked up the extra 
term D(x — y ) = [</> + (x), 4>~{y)\ which is the propagator we met in (2.90). So for x° > y° 
we have 


T (j>{x)(t>{y) =: <t>{x)(f>(y) : + D(x - y) (3.34) 

Meanwhile, for y° > x°, we may repeat the calculation to find 

T (j){x)(j)(y) =: (f>{x)(f>{y) : + D(y - x) (3.35) 

So putting this together, we have the final expression 

T <f>{x)<l>{y) =: <!>{x)(j>{y) : + A F (x - y) (3.36) 

where Ap(x — y) is the Feynman propagator defined in (2.93), for which we have the 
integral representation 


A F (x-y) = 


d A k 


: p ik-(x—y) 


ie 


(27r) 4 k 2 — m 2 + ie 


(3.37) 


Let me reiterate a comment from Section 2: although T (f)(x)(j>{y) and : (p(x)(p(y) : are 
both operators, the difference between them is a c-number function, A f{x — y). 


Definition: We define the contraction of a pair of fields in a string of operators 
... (j>(x i) .. . 4>(x2) ... to mean replacing those operators with the Feynman propaga¬ 
tor, leaving all other operators untouched. We use the notation, 

... 0(xi)... <f>(x 2 ) ■ ■ ■ (3.38) 

to denote contraction. So, for example, 

<t>{x)(f>{y) = A f (x - y) (3.39) 
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A similar discussion holds for complex scalar fields. We have 

T^{x)'ip\y) =: 0(x)0 t (y) : +A F (x - y) 


(3.40) 


prompting us to define the contraction 

0(x)0t(y) = A F (x — y ) and ip(x)ip(y) = ^ (x)^ (y) = 0 (3.41) 

3.3.2 Wick’s Theorem 

For any collection of fields 0i = 00ci), 02 = 00c 2 ), etc, we have 

T(0 1 ... <j> n ) =: 0!... 0 n : + : all possible contractions : (3.42) 

To see what the last part of this equation means, let’s look at an example. For n = 4, 
the equation reads 


r a r > 

T ( 01020304 ) = : 01020304 : +0102 : 0304 : +0103 : 0204 : + four similar terms 

+ 0102 0304 + 0103 0204 + 0104 0203 (3.43) 


Proof: The proof of Wick’s theorem proceeds by induction and a little thought. It’s 
true for n = 2. Suppose it’s true for 0 2 ... 0 n and now add 0i. We’ll take a;) 1 > 
for all k — 2, ... ,n. Then we can pull 0i out to the left of the time ordered product, 
writing 

T (0102 • • • 0n) = (0i~ + 0r) 0 02 • • • 0n : + : contractions :) (3.44) 

The 0]" term stays where it is since it is already normal ordered. But in order to write 
the right-hand side as a normal ordered product, the 0+ term has to make its way 
past the crowd of operators. Each time it moves past 0^, we pick up a factor of 

010A, = Aip(a;i — Xk) from the commutator. (Try it!) □ 


3.3.3 An Example: Nucleon Scattering 

Let’s look at 00 —y 00 scattering. We have the initial and final states 

N> = l°> = IP1.P2) 

I/) = 02^021% |0) = < 3 ' 45) 

We can then look at the expansion of (/|S'|i). In fact, we really want to calculate 
(/| S — 1 |f) since we’re not interested in situations where no scattering occurs. At 
order g 2 we have the term 

- y~ j d A xid A x 2 T (0 t (xi)0(a;i)0(a;i)0 t (a;2)0(a;2)0(a;2)) (3.46) 
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Now, using Wick’s theorem we see there is a piece in the string of operators which looks 
like 

: ip\xi)ip(xi)ip\x 2 )ip(x 2 ) : <p(x 1 )<p(x 2 ) (3.47) 


which will contribute to the scattering because the two ip fields annihilate the ip parti¬ 
cles, while the two ip' fields create ip particles. Any other way of ordering the ip and 
ip' fields will give zero contribution. This means that we have 


{p'i,P 2 1 ■'Ip\x 1 )ip(x 1 )ip\x 2 )ip(x 2 ) : \pi,p 2 ) 

= (p'l,P 2 \ V’ t (^l)'0 t (^2) |0) (0\ip(x 1 )ip(x 2 ) \pi,p 2 ) 


^ e ip' 1 -x 1 +ip' 2 -X2 _J_ e ip' 1 -X2+ip' 2 -x 1 ^ ^ e ~ipi-x 1 -ip 2 -X2 


_J_ e ~ipvx 2 -ip2-xi 


) 


— gixi-(pi-Pl)+“2-(P2-P2) _|_ e ixi-(p' 2 -pi)+iX 2 -(p' 1 -p 2 ) 


+ (xi O x 2 ) 


(3.48) 


where, in going to the third line, we’ve used the fact that for relativistically normalized 
states, 


(0| ip(x) \p) — e w ' x (3.49) 

Now let’s insert this into (3.46), to get the expression for (/1 S |i) at order g 2 , 

(— ig ) 2 f ,a r a a , si f d 4 k i e ik(x 1 -x 2 ) 


d A xid A x 2 \e 1 '" + e 1 '" + (xi <R- x 2 )] 


(27r) 4 k 2 — m 2 + ie 


(3.50) 


where the expression in square brackets is (3.48), while the final integral is the (p 
propagator which comes from the contraction in (3.47). Now the (aq -H- x 2 ) terms 
double up with the others to cancel the factor of 1/2 out front. Meanwhile, the X\ and 
x 2 integrals give delta-functions. We’re left with the expression 




d A k i{ 2ny 


(27r) 4 k 2 — m 2 + ie 


[<5 (4) G°i ~ Pi + k) h (4) (p' 2 -p2-k) 


+ <5 (4) (p' 2 —pi + k) S ( ' A \p[ - p 2 - k)] (3.51) 

Finally, we can trivially do the d A k integral using the delta-functions to get 




— + 


(pi ~ p[) 2 ~ m 2 + (pi ~ P 2 ) 2 ~ m2 + 


(27t) 4 5^ A \pi + p 2 — Pi — P 2 ) 


In fact, for this process we may drop the +ie terms since the denominator is never 
zero. To see this, we can go to the center of mass frame, where pi = — p 2 and, by 
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momentum conservation, |pj| = 1/7/1. This ensures that the 4-momentum of the meson 
is k = (0 ,p — p'), so k 2 < 0. We therefore have the end result, 


-wf 


+ 


(pi — p[) 2 — m 2 {pi — P 2) 2 — m 2 


(2vr) 4 5 {A \pi +P 2 - p'i ~ P 2 ) (3.52) 


We will see another, much simpler way to reproduce this result shortly using Feynman 
diagrams. This will also shed light on the physical interpretation. 


This calculation is also relevant for other scattering processes, such as 00 —> -0-0, 
'00 —>■ 00. Each of these comes from the term (3.48) in Wick’s theorem. However, we 
will never find a term that contributes to scattering 00 —>- 00, for this would violate 
the conservation of Q charge. 


Another Example: Meson-Nucleon Scattering 

If we want to compute 00 — Y ifjc/) scattering at order g 2 , we would need to pick out the 
term 

/-' N —-s 

: 0 t (xi)0(xi)0(x 2 )0(x 2 ) : 0(xi)0 t 0c 2 ) (3.53) 

and a similar term with 0 and 0' exchanged. Once more, this term also contributes to 
similar scattering processes, including 00 —)■ 00 and 00 —> 00. 


3.4 Feynman Diagrams 

“Like the silicon chips of more recent years, the Feynman diagram was 

bringing computation to the masses.” T n , . 

Julian bcnwmger 

As the above example demonstrates, to actually compute scattering amplitudes using 
Wick’s theorem is rather tedious. There’s a much better way. It requires drawing pretty 
pictures. These pictures represent the expansion of ( f\S |i) and we will learn how to 
associate numbers (or at least integrals) to them. These pictures are called Feynman 
diagrams. 

The object that we really want to compute is (/1 S' —1 |i), since we’re not interested in 
processes where no scattering occurs. The various terms in the perturbative expansion 
can be represented pictorially as follows 

• Draw an external line for each particle in the initial state |i) and each particle 
in the final state |/). We’ll choose dotted lines for mesons, and solid lines for 
nucleons. Assign a directed momentum p to each line. Further, add an arrow to 
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solid lines to denote its charge; we’ll choose an incoming (outgoing) arrow in the 
initial state for ^ (?/?). We choose the reverse convention for the final state, where 
an outgoing arrow denotes 


• Join the external lines together with trivalent vertices 



Each such diagram you can draw is in 1-1 correspondence with the terms in the 
expansion of (f \ S — 1 |i). 


3.4.1 Feynman Rules 

To each diagram we associate a number, using the Feynman rules 
• Add a momentum k to each internal line 


• To each vertex, write down a factor of 

(~ig) (2tt) 4 J (4) (^2 h) (3.54) 

i 

where )T) ki is the sum of all momenta flowing into the vertex. 


• For each internal dotted line, corresponding to a 0 particle with momentum k, 
we write down a factor of 

r d A k i 
J (27t) 4 k 2 — m 2 + ie 



We include the same factor for solid internal ^ lines, with m replaced by the 
nucleon mass M. 
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3.5 Examples of Scattering Amplitudes 

Let’s apply the Feynman rules to compute the amplitudes for various processes. We 
start with something familiar: 


Nucleon Scattering Revisited 

Let’s look at how this works for the ifnjj —> scattering at order g 2 . We can write 

down the two simplest diagrams contributing to this process. They are shown in Figure 
9. 



Figure 9: The two lowest order Feynman diagrams for nucleon scattering. 


Applying the Feynman rules to these diagrams, we get 
1 1 




+ 


(pi — p[) 2 — m 2 (pi — P 2) 2 ~ m2 


(2vr ) 4 <f (4) (pi +P 2 - p[ - p' 2 ) (3.56) 


which agrees with the calculation (3.51) that we performed earlier. There is a nice 
physical interpretation of these diagrams. We talk, rather loosely, of the nucleons 
exchanging a meson which, in the first diagram, has momentum k = (pi—ff) = G» 2 — ^ 2 )- 
This meson doesn’t satisfy the usual energy dispersion relation, because k 2 7 ^ m 2 : the 
meson is called a virtual particle and is said to be off-shell (or, sometimes, off mass- 
shell). Heuristically, it can’t, live long enough for its energy to be measured to great 
accuracy. In contrast, the momentum on the external, nucleon legs all satisfy p 2 = M 2 , 
the mass of the nucleon. They are on-shell. One final note: the addition of the two 
diagrams above ensures that the particles satisfy Bose statistics. 


There are also more complicated diagrams which will contribute to the scattering 
process at higher orders. For example, we have the two diagrams shown in Figures 
10 and 11, and similar diagrams with p\ and p' 2 exchanged. Using the Feynman rules, 
each of these diagrams translates into an integral that we will not attempt to calculate 
here. And so we go on, with increasingly complicated diagrams, all appearing at higher 
order in the coupling constant g. 
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Figure 10: A contribution at O(g^). Figure 11: A contribution at 0(g®) 


Amplitudes 

Our final result for the nucleon scattering amplitude (/| S — 1 |i) at order g 2 was 


-igf 


+ 


L(pi — p[) 2 — m 2 (p] — P 2) 2 — m 2 


(27r) 4 5 (4) (pi +p 2 -p[ -P 2 ) 


The 5-function follows from the conservation of 4-momentum which, in turn, follows 
from spacetime translational invariance. It is common to all S-matrix elements. We will 
define the amplitude Afi by stripping off this momentum-conserving delta-function, 


(/I s - 1 |i) = iAfi (2jr) 4 <5 (J, (p F - Pi) 


(3.57) 


where pi ( pp ) is the sum of the initial (final) 4-momenta, and the factor of i out front 
is a convention which is there to match non-rclativistic quantum mechanics. We can 
now refine our Feynman rules to compute the amplitude iAfi itself: 


• Draw all possible diagrams with appropriate external legs and impose 4-momentum 
conservation at each vertex. 


• Write down a factor of (— ig) at each vertex. 

• For each internal line, write down the propagator 

• Integrate over momentum k flowing through each loop f d 4 k/( 2tt) 4 . 

This last step deserves a short explanation. The diagrams we’ve computed so far have 
no loops. They are tree level diagrams. It’s not hard to convince yourself that in 
tree diagrams, momentum conservation at each vertex is sufficient to determine the 
momentum flowing through each internal line. For diagrams with loops, such as those 
shown in Figures 10 and 11, this is no longer the case. 
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Figure 12: The two lowest order Feynman diagrams for nucleon to meson scattering. 


Nucleon to Meson Scattering 


Let’s now look at the amplitude for a nucleon-anti-nucleon pair to annihilate into a 
pair of mesons: —> (j)(j). The simplest Feynman diagrams for this process are shown 

in Figure 12 where the virtual particle in these diagrams is now the nucleon i/j rather 
than the meson ef). This fact is reflected in the denominator of the amplitudes which 
are given by 


iA = ( -ig) 2 

As in (3.52), we’ve dropped 
vanishes. 


_(pi-p{) 2 -M 2 + ( Pl -p') 2 -M 2 _ ( 3 ' 58) 

the ie from the propagators as the denominator never 


Nucleon-Anti-Nucleon Scattering 



Figure 13: 


The two lowest order Feynman diagrams for nucleon-anti-nucleon scattering. 


For the scattering of a nucleon and an anti-nucleon, ipip —> ipij}, the Feynman 
diagrams are a little different. At lowest order, they are given by the diagrams of 
Figure 13. It is a simple matter to write down the amplitude using the Feynman rules, 


iA = {-igf 


+ 


(pi — p[) 2 — m 2 (pi + p 2 ) 2 — m 2 + ie 


(3.59) 
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Notice that the momentum dependence in the sec¬ 
ond term is different from that of nucleon-nucleon 
scattering (3.56), reflecting the different Feynman di¬ 
agram that contributes to the process. In the center 
of mass frame, p\ = — p 2 , the denominator of the sec¬ 
ond term is 4(M 2 + p 2 ) — m 2 . If m < 2M, then this 
term never vanishes and we may drop the ie. In con¬ 
trast, if m > 2 M, then the amplitude corresponding 
to the second diagram diverges at some value of p. 

In this case it turns out that we may also neglect the 
ie term, although for a different reason: the meson 
is unstable when m > 2 M, a result we derived in 
(3.30). When correctly treated, this instability adds Figure 14: 

a finite imaginary piece to the denominator which 

overwhelms the ie. Nonetheless, the increase in the scattering amplitude which we see 
in the second diagram when 4(M 2 +p 2 ) = m 2 is what allows us to discover new parti¬ 
cles: they appear as a resonance in the cross section. For example, the Figure 14 shows 
the cross-section (roughly the amplitude squared) plotted vertically for e + e _ —> p + 
scattering from the ALEPH experiment in CERN. The horizontal axis shows the center 
of mass energy. The curve rises sharply around 91 GeV, the mass of the Z-boson. 

Meson Scattering 

For ej)(j) —>■ 00, the simplest diagram we can write 
down has a single loop, and momentum conservation at 
each vertex is no longer sufficient to determine every 
momentum passing through the diagram. We choose 
to assign the single undetermined momentum k to the 
right-hand propagator. All other momenta are then de¬ 
termined. The amplitude corresponding to the diagram 
shown in the figure is 

\ 4 f d 4 k _1_ 

^ J (27r) 4 (k 2 — M 2 + ie)((k + p [) 2 — M 2 + ie) 

1 

((k + p[ — pi) 2 — M 2 + ie)((k — P 2) 2 — M 2 + ie) 

These integrals can be tricky. For large k, this integral goes as f cftk/k 8 , which is at 
least convergent as k —» 00 . But this won’t always be the case! 


>v -- 


p 2 
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3.5.1 Mandelstam Variables 


We see that in many of the amplitudes above — in particular those that include the 
exchange of just a single particle — the same combinations of momenta are appear¬ 
ing frequently in the denominators. There are standard names for various sums and 
differences of momenta: they are known as Mandelstam variables. They are 

s = (Pi + P 2) 2 — (p'i + P 2) 2 

t = (pi-p[) 2 = (P2-p 2 ) 2 (3.60) 

U = (pi~ P' 2 ) 2 = (P2 - AY 

where, as in the examples above, p\ and P 2 are the momenta of the two initial particles, 
and p\ and p' 2 are the momenta of the final two particles. We can define these variables 
whether the particles involved in the scattering are the same or different. To get a feel 
for what these variables mean, let’s assume all four particles are the same. We sit in 
the center of mass frame, so that the initial two particles have four-momenta 

p 1 = (E,0,0,p) and P 2 = (E, 0, 0, —p) (3.61) 

The particles then scatter at some angle 9 and leave with momenta 

p[ — (E, 0,psm9,pcos9) and p' 2 — (E, 0, — psin 9, — pcos 9) (3.62) 

Then from the above definitions, we have that 

s = AE 2 and t =—2p 2 (1 — cos9) and u = — 2p 2 (l + cos 9) (3.63) 

The variable s measures the total center of mass energy of the collision, while the 
variables t and u are measures of the momentum exchanged between particles. (They 
are basically equivalent, just with the outgoing particles swapped around). Now the 
amplitudes that involve exchange of a single particle can be written simply in terms of 
the Mandelstam variables. For example, for nucleon-nucleon scattering, the amplitude 
(3.56) is schematically A ~ (t — m 2 ) -1 + {u — m 2 )~ l . For the nucleon-anti-nucleon 
scattering, the amplitude (3.59) is A ~ (t — m 2 )^ 1 + (s — m 2 ) -1 . We say that the 
first case involves “t,-cha nn el” and “u-channel” diagrams. Meanwhile the nucleon-anti- 
nucleon scattering is said to involve “t-channel” and “s-channel” diagrams. (The first 
diagram indeed includes a vertex that looks like the letter “T”). 

Note that there is a relationship between the Mandelstam variables. When all the 
masses are the same we have s + t + u = AM 2 . When the masses of all 4 particles differ, 
this becomes s + t + u = Yhi M?. 
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3.5.2 The Yukawa Potential 

So far we’ve computed the quantum amplitudes for various scattering processes. But 
these quantities are a little abstract. In Section 3.6 below (and again in next term’s 
“Standard Model” course) we’ll see how to turn amplitudes into measurable quantities 
such as cross-sections, or the lifetimes of unstable particles. Here we’ll instead show 
how to translate the amplitude (3.52) for nucleon scattering into something familiar 
from Newtonian mechanics: a potential, or force, between the particles. 

Let’s start by asking a simple question in classical field theory that will turn out to 
be relevant. Suppose that we have a fixed ^-function source for a real scalar field 0, 
that persists for all time. What is the profile of 0(f)? To answer this, we must solve 
the static Klein-Gordon equation, 

-V 2 0 + m 2 0 = 5 ( 3 ) (f) (3.64) 


We can solve this using the Fourier transform, 

/d 3 k 

(3.65) 

Plugging this into (3.64) tells us that [k 2 + m 2 )0(/c) = 1, giving us the solution 


0(f) 


d 3 k e lk ' x 

( 2 vr ) 3 £2 _|_ m 2 


(3.66) 


Let’s now do this integral. Changing to polar coordinates, and writing k ■ x 
we have 


0 (f) 


1 r k 2 2 sin kr 

(27t) 2 J 0 k 2 + m 2 kr 

1 C + °° k sin kr 

(2tt) 2 r J _ 00 k 2 + m 2 

1 \ f +ao dk ke ikr 

-Re / -—=-- 

27rr L ./-00 27t* k z + m z 


kr cos 9 , 


(3.67) 


We compute this last integral by closing the contour in the upper half plane k —> +ioo , 
picking up the pole at k = +im. This gives 

0 (f) = J_ e ~ mr (3.68) 

' 47rr v ' 

The field dies off exponentially quickly at distances 1/m, the Compton wavelength of 
the meson. 
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Now we understand the profile of the 0 field, what does this have to do with the force 
between -0 particles? We do very similar calculations to that above in electrostatics 
where a charged particle acts as a 5-function source for the gauge potential: — V 2 A 0 = 
S^ 3 \x), which is solved by Aq = l/Anr. The profile for Aq then acts as the potential 
energy for another charged (test) particle moving in this background. Can we give the 
same interpretation to our scalar field? In other words, is there a classical limit of the 
scalar Yukawa theory where the 0 particles act as 5-function sources for 0, creating 
the profile (3.68)? And, if so, is this profile then felt as a static potential? The answer 
is essentially yes, at least in the limit M m. But the correct way to describe the 
potential felt by the 0 particles is not to talk about classical fields at all, but instead 
work directly with the quantum amplitudes. 


Our strategy is to compare the nucleon scattering amplitude (3.52) to the corre¬ 
sponding amplitude in non-relativistic quantum mechanics for two particles interacting 
through a potential. To make this comparison, we should first take the non-relativistic 
limit of (3.52). Let’s work in the center of mass frame, with p = p\ — —p 2 and 
p' = Pi = —p 2 - The non-relativistic limit means |p| <C M which, by momentum 
conservation, ensures that \p'\ <C M. In fact one can check that, for this particular 
example, this limit doesn’t change the scattering amplitude (3.52): it’s given by 


iA 


+ig- 


+ 


(p — p') 2 + m 2 {p + p') 2 + m 


(3.69) 


How do we compare this to scattering in quantum mechanics? Consider two particles, 
separated by a distance r, interacting through a potential U(r). In non-relativistic 
quantum mechanics, the amplitude for the particles to scatter from momentum states 
±p into momentum states ±p' can be computed in perturbation theory, using the 
techniques described in Section 3.1. To leading order, known in this context as the 
Born approximation, the amplitude is given by 


(p'| U(r) | p) 


j d 3 rU(r)e~ i{ P-^ ) ? 


(3.70) 


There’s a relative factor of (2M) 2 that arises in comparing the quantum field theory 
amplitude A to (p'\ U(r) \p), that can be traced to the relativistic normalization of the 
states \p \, p 2 ) • (If is also necessary to get the dimensions of the potential to work out 
correctly). Including this factor, and equating the expressions for the two amplitudes, 
we get 


d 3 r U(r) = 


-A 2 


(p — p') 2 + m 2 


(3.71) 







where we’ve introduced the dimensionless parameter A = g/2M. We can trivially invert 
this to find, 


U(r) = -A 2 


/ 


d 3 p e& ? 
(27t) 3 p 2 + m 2 


(3.72) 


But this is exactly the integral (3.66) we just did in the classical theory. We have 



(3.73) 


—mr 


e 


This is the Yukawa potential. The force has a range 1/m, the Compton wavelength of 
the exchanged particle. The minus sign tells us that the potential is attractive. 

Notice that quantum held theory has given us an entirely new perspective on the 
nature of forces between particles. Rather than being a fundamental concept, the force 
arises from the virtual exchange of other particles, in this case the meson. In Section 6 
of these lectures, we will see how the Coulomb force arises from quantum held theory 
due to the exchange of virtual photons. 

We could repeat the calculation for nucleon-anti-nucleon scattering. The amplitude 
from held theory is given in (3.59). The hrst term in this expression gives the same 
result as for nucleon-nucleon scattering with the same sign. The second term vanishes in 
the non-relativisitic limit (it is an example of an interaction that doesn’t have a simple 
Newtonian interpretation). There is no longer a factor of 1/2 in (3.70), because the 
incoming/outgoing particles are not identical, so we learn that the potential between 
a nucleon and anti-nucleon is again given by (3.73). This reveals a key feature of 
forces arising due to the exchange of scalars: they are universally attractive. Notice 
that this is different from forces due to the exchange of a spin 1 particle — such as 
electromagnetism — where the sign hips when we change the charge. However, for 
forces due to the exchange of a spin 2 particle — i.e. gravity — the force is again 
universally attractive. 

3.5.3 </> 4 Theory 

Let’s briehy look at the Feynman rules and scattering amplitudes for the interaction 
Hamiltonian 



4 


(3.74) 


The theory now has a single interaction vertex, which comes with a factor of (— iX), 
while the other Feynman rules remain the same. Note that we assign (— iX) to the 





vertex rather than (—iA/4!). To see why this is, we can look at (jxf> —> (jxj) scattering, 
which has its lowest contribution at order A, with the term 

(Pi,P 2 1 : <l>{x)<l>{x)<j)(x)<j)(x) : \pi,p 2 ) (3.75) 

Any one of the helds can do the job of annihilation or creation. This gives 4! different 
contractions, which cancels the 1/4! sitting out front. 


Feynman diagrams in the 0 4 theory sometimes come with extra 
combinatoric factors (typically 2 or 4) which are known as symmetry 
factors that one must take into account. For more details, see the book 
by Peskin and Schroeder. 
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Using the Feynman rules, the scattering amplitude for <jxj) —» (jxj) is Figure 16: 
simply i,A = —i\. Note that it doesn’t depend on the angle at which 
the outgoing particles emerge: in </> 4 theory the leading order two-particle scattering 
occurs with equal probability in all directions. Translating this into a potential between 
two mesons, we have 


U (r) 


A 

(2m) 2 


p+ip-r 

(2tt) 3 


A 

(2 m ) 2 


5 {3) (r) 


(3.76) 


So scattering in 0 4 theory is due to a ^-function potential. The particles don’t know 
what hit them until it’s over. 


3.5.4 Connected Diagrams and Amputated Diagrams 

We’ve seen how one can compute scattering amplitudes by writing down all Feynman 
diagrams and assigning integrals to them using the Feynman rules. In fact, there are 
a couple of caveats about what Feynman diagrams you should write down. Both of 
these caveats are related to the assumption we made earlier that “initial and final states 
are eigenstates of the free theory” which, as we mentioned at the time, is not strictly 
accurate. The two caveats which go some way towards ameliorating the problem are 
the following 

• We consider only connected Feynman diagrams, where every part of the diagram 
is connected to at least one external line. As we shall see shortly, this will be 
related to the fact that the vacuum |0) of the free theory is not the true vacuum 
112) of the interacting theory. An example of a diagram that is not connected is 
shown in Figure 17. 
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• We do not consider diagrams with loops on external lines, for example the diagram 
shown in the Figure 18. We will not explain how to take these into account in this 
course, but you will discuss them next term. They are related to the fact that the 
one-particle states of the free theory are not the same as the one-particle states of 
the interacting theory. In particular, correctly dealing with these diagrams will 
account for the fact that particles in interacting quantum field theories are never 
alone, but surrounded by a cloud of virtual particles. We will refer to diagrams 
in which all loops on external legs have been cut-off as “amputated”. 



Figure 17: A disconnected diagram. Figure 18: An un-amputated diagram 


3.6 What We Measure: Cross Sections and Decay Rates 

So far we’ve learnt to compute the quantum amplitudes for particles decaying or scat¬ 
tering. As usual in quantum theory, the probabilities for things to happen are the 
(modulus) square of the quantum amplitudes. In this section we will compute these 
probabilities, known as decay rates and cross sections. One small subtlety here is that 
the S-matrix elements (f\S — l |i) all come with a factor of (pp — pi), so we 

end up with the square of a delta-function. As we will now see, this comes from the 
fact that we’re working in an infinite space. 

3.6.1 Fermi’s Golden Rule 

Let’s start with something familiar and recall how to derive Fermi’s golden rule from 
Dyson’s formula. For two energy eigenstates | m) and | n), with E m ^ E n , we have to 
leading order in the interaction, 

(m\ U ( t ) | n) = — i (m\ f dt Hj(t) \n) 

Jo 

= —i (m\ H mt | n) f dt 1 e luJt ' 

Jo 

piut _ i 

= - {m\ H int | n) - (3.77) 
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where c o = E m — E n . This gives us the probability for the transition from | n) to | m) in 
time t, as 

P n _> m (t) = \(m\u ( t ) I n) 1 2 = 2| (m\ H int \ n) | 2 ( 

The function in brackets is plotted in Figure 19 for fixed t. 

We see that in time t, most transitions happen in a region 
between energy eigenstates separated by A E = 2n/t. As 
t —» oo, the function in the figure starts to approach a delta- 
function. To find the normalization, we can calculate 

— cos cot \ 

-2- = nt 

co z J 

—> 7rtS(uj) as t —> oo 

Consider now a transition to a cluster of states with density 
p(E). In the limit t —> oo, we get the transition probability 

Pn^m j dE m p^Efy^ 2| (m\ H[ n i | n) | 

->• 27T I (m\ H int \n) \ 2 p(E n )t (3.79) 

which gives a constant probability for the transition per unit time for states around 
the same energy E n E m — E. 





Figure 19: 


cos ut 


UJ 


(3.78) 


Pn—trn 2n\ (m\ H int | n) \ 2 p(E) (3.80) 

This is Fermi’s Golden Rule. 

In the above derivation, we were fairly careful with taking the limit as t —> oo. 
Suppose we were a little sloppier, and first chose to compute the amplitude for the 
state |n) at t —> —oo to transition to the state | m) at t —> +oo. Then we get 

rt =+oo 

—i (m\ / Hi(t ) |?i) = —i ( m\ H int \n) 2i r5(a;) (3.81) 

J t =—OO 

Now when squaring the amplitude to get the probability, we run into the problem of 
the square of the delta-function: P n ^>m = | (m\ H mt \n) \ 2 (2 tt) 2 5(u>) 2 . Tracking through 
the previous computations, we realize that the extra infinity is coming because P m ^ n 
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is the probability for the transition to happen in infinite time t —>• oo. We can write 
the delta-fnnctions as 

(27r) 2 5(cn) 2 = (2ir)5(u) T (3.82) 

where T is shorthand for t —> oo (we used a very similar trick when looking at the 
vacuum energy in (2.25)). We now divide out by this power of T to get the transition 
probability per unit time, 


P n ^ m = 2tt| (m \ H int | n) | 2 5(w) (3.83) 

which, after integrating over the density of final states, gives us back Fermi’s Golden 
rule. The reason that we’ve stressed this point is because, in our field theory calcula¬ 
tions, we’ve computed the amplitudes in the same way as (3.81), and the square of the 
(^-functions will just be re-interpreted as spacetime volume factors. 


3.6.2 Decay Rates 

Let’s now look at the probability for a single particle |i) of momentum pj (I=initial) 
to decay into some number of particles |/) with momentum p i: and total momentum 
Pf = YhiPi- This is given by 


P = 


(f\S\i) 


(/I/) <*l*> 

Our states obey the relativistic normalization formula (2.65), 

(i\i) = (2n) 3 2Ep I S^(0)=2E fl V 


(3.84) 


(3.85) 


where we have replaced <5®(0) by the volume of 3-space. Similarly, 

</i /> = n 2E <>y (3-86) 

final states 


If we place our initial particle at rest, so pj — 0 and Ep : = m, we get the probability 
for decay 


P 


2mV 


(2 OW-pUCT n ^7 

final states 


(3.87) 


where, as in the second derivation of Fermi’s Golden Rule, we’ve exchanged one of the 
delta-functions for the volume of spacetime: (27 t) 4 5^(0) = VT. The amplitudes Afi 
are, of course, exactly what we’ve been computing. (For example, in (3.30), we saw 
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that A = —g for a single meson decaying into two nucleons). We can now divide out 
by T to get the transition function per unit time. But we still have to worry about 
summing over all final states. There are two steps: the first is to integrate over all 
possible momenta of the final particles: V f d 3 pi/(2n) 3 . The factors of spatial volume 
V in this measure cancel those in (3.87), while the factors of 1/2 Ep i in (3.87) conspire 
to produce the Lorentz invariant measure for 3-momentum integrals. The result is an 
expression for the density of final states given by the Lorentz invariant measure 

dU = (2n)^ 4 \p F - Pl ) J] 44^7 (3-88) 

final states ' ' P' 

The second step is to sum over all final states with different numbers (and possibly 
types) of particles. This gives us our final expression for the decay probability per unit 
time, T = P. 

r=2_ 7 I' \A f ,\ 2 dU. (3.89) 

final states 

T is called the width of the particle. It is equal to the reciprocal of the half-life r = 1/T. 
3.6.3 Cross Sections 

Collide two beams of particles. Sometimes the particles will hit and bounce off each 
other; sometimes they will pass right through. The fraction of the time that they collide 
is called the cross section and is denoted by a. If the incoming flux F is defined to 
be the number of incoming particles per area per unit time, then the total number of 
scattering events N per unit time is given by, 

N = Fa (3.90) 

We would like to calculate a from quantum held theory. In fact, we can calculate a more 
sensitive quantity da known as the differential cross section which is the probability 
for a given scattering process to occur in the solid angle (#,</>). More precisely 

Differential Probability 1 1 , 2 

° Unit Time x Unit Flux 4E 1 E 2 V F fl 

where we’ve used the expression for probability per unit time that we computed in the 
previous subsection. E\ and E 2 are the energies of the incoming particles. We now 
need an expression for the unit flux. For simplicity, let’s sit in the center of mass frame 
of the collision. We’ve been considering just a single particle per spatial volume V, 
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meaning that the flux is given in terms of the 3-velocities hj as F = If?! — v 2 \/V. This 
then gives, 


da 


1 1 

AEiE 2 \vi — v 2 


\A fi \ 2 dU 


(3.92) 


If you want to write this in terms of momentum, then recall from your course on special 
relativity that the 3-velocities v t are related to the momenta by v = p/m\J 1 — v 2 = 

p/p°. 


Equation (3.92) is our final expression relating the S-matrix to the differential cross 
section. You may now take your favorite scattering amplitude, and compute the proba¬ 
bility for particles to fly out at your favorite angles. This will involve doing the integral 
over the phase space of final states, with measure dll. Notice that different scattering 
amplitudes have different momentum dependence and will result in different angular 
dependence in scattering amplitudes. For example, in </> 4 theory the amplitude for tree 
level scattering was simply A = —A. This results in isotropic scattering. In contrast, 
for nucleon-nucleon scattering we have schematically A ~ (t — m 2 ) -1 + (u — m 2 )~h 
This gives rise to angular dependence in the differential cross-section, which follows 
from the fact that, for example, t = —2|p| 2 (l — cos 9), where 9 is the angle between the 
incoming and outgoing particles. 


3.7 Green’s Functions 

So far we’ve learnt to compute scattering amplitudes. These are nice and physical (well 
- they’re directly related to cross-sections and decay rates which are physical) but there 
are many questions we want to ask in quantum field theory that aren’t directly related 
to scattering experiments. For example, we might want to compute the viscosity of 
the quark gluon plasma, or the optical conductivity in a tentative model of strange 
metals, or figure out the non-Gaussianity of density perturbations arising in the CMB 
from novel models of inflation. All of these questions are answered in the framework of 
quantum field theory by computing elementary objects known as correlation functions. 
In this section we will briefly define correlation functions, explain how to compute them 
using Feynman diagrams, and then relate them back to scattering amplitudes. We’ll 
leave the relationship to other physical phenomena to other courses. 

We’ll denote the true vacuum of the interacting theory as |fl). We’ll normalize El 
such that 


H |fi) = 0 


(3.93) 
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and (f2| 0) = 1. Note that this is different from the state we’ve called |0) which is the 
vacuum of the free theory and satisfies Hq |0) = 0. Define 

G {n \x 1 ,... ,x n ) = (QlT(/) H (x 1 )...^ H (x n )in) (3.94) 

where (pH is 0 in the Heisenberg picture of the full theory, rather than the interaction 
picture that we’ve been dealing with so far. The G are called correlation functions, 
or Green’s functions. There are a number of different ways of looking at these objects 
which tie together nicely. Let’s start by asking how to compute G^ using Feynman 
diagrams. We prove the following result 

Claim: We use the notation 0i = 0(aq), and write (pin to denote the held in the 
Heisenberg picture, and cpu to denote the held in the interaction picture. Then 

. ..,*„) = <Si| Kh in) = (3.95) 

where the operators on the right-hand side are evaluated on | 0 ), the vacuum of the free 
theory. 

Proof: Take ti > t 2 > • • • > t n . Then we can drop the T and write the numera¬ 
tor of the right-hand side as 

(0| [M+OO, fl)017 U (H, t 2 ) 027 ■ • • 0n7 Ul(t n , -OO) |0) 

We’ll use the factors of Ui(tk,tk+i) = T exp (—i f]' k+l Hj) to convert each of the 07 into 
077 and we choose operators in the two pictures to be equal at some arbitrary time to- 
Then we can write 

(0| Ui(+00, tl)0!7 U (ti, t 2 ) 027 •• • 0n7 U^tn, -OO ) |0) 

= (0| Uji + OO, t O )0l77 • • • 0n77 bh(t 0 , ~Oo) |0) 

Now let’s deal with the two remaining U (to, ±oo) at either end of the string of operators. 
Consider an arbitrary state T) and look at 

(T| Uj(t, -oo) | 0 ) = (T| U(t, -oo) | 0 ) (3.96) 

where f/(t, —oo) is the Schrodinger evolution operator, and the equality above follows 
because H 0 |0) = 0. Now insert a complete set of states, which we take to be energy 
eigenstates of H = H 0 + H - mt , 

(T| t/(t, — oo) | 0 ) = (T| f/(t, — oo) |fi) <fi| + Y, In) (n| | 0 ) 

n^O 

= (T|H)(D|0)+ lim VeW m n ) (n|0) (3.97) 

t'^ — CXj ' 

n^O 
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But the last term vanishes. This follows from the Ricmann-Lebesgue lemma which says 
that for any well-behaved function 

lim [ dxf(x)e ztix = 0 (3.98) 

^ 00 J a 

Why is this relevant? The point is that the Y2 n (3-97) is really an integral f dn , 
because all states are part of a continuum due to the momentum. (There is a caveat 
here: we want the vacuum |f2) to be special, so that it sits on its own, away from the 
continuum of the integral. This means that we must be working in a theory with a 
mass gap - i.e. with no massless particles). So the Riemann-Lebesgue lemma gives us 


lim (V\U(t,t')\0) = (V\n) (tt\0) (3.99) 

t'—>•—00 

(Notice that to derive this result, Peskin and Schroeder instead send t —> —00 in a 
slightly imaginary direction, which also does the job). We now apply the formula 
(3.99), to the top and bottom of the right-hand side of (3.95) to find 

(Q|n> (n|o) 

<o| n> <n| !J> (fi| 0) 1 J 

which, using the normalization (fi| £2) = 1, gives us the left-hand side, completing the 
proof. □. 


3.7.1 Connected Diagrams and Vacuum Bubbles 

We’re getting closer to our goal of computing the Green’s functions since we 
can compute both (0| T(j)j(xi)... 4>i(x n ) S |0) and (0|S'|0) using the same methods 
we developed for S-matrix elements; namely Dyson’s formula and Wick’s theorem or, 
alternatively, Feynman diagrams. But what about dividing one by the other? What’s 
that all about? In fact, it has a simple interpretation. For the following discussion, 
we will work in </> 4 theory. Since there is no ambiguity in the different types of line 
in Feynman diagrams, we will represent the (j) particles as solid lines, rather than the 
dashed lines that we used previously. Then we have the diagramatic expansion for 

(o| s |0). 

<oisio>=i+g+( 8 + ((j) + 8 8 ) + (3 - 101) 

These diagrams are called vacuum bubbles. The combinatoric factors (as well as the 
symmetry factors) associated with each diagram are such that the whole series sums 


77 



to an exponential, 


(0| S |0) = exp (8 + 8 + dD + -) (3.102) 

So the amplitude for the vacuum of the free theory to evolve into itself is (0| S |0) = 
exp(all distinct vacuum bubbles). A similar combinatoric simplification occurs for generic 
correlation functions. Remarkably, the vacuum diagrams all add up to give the same 
exponential. With a little thought one can show that 

(0| T(j>i... (f> n S |0) = connected diagrams j (0| S |0) (3.103) 

where “connected” means that every part of the diagram is connected to at least one 
of the external legs. The upshot of all this is that dividing by (0| S |0) has a very nice 
interpretation in terms of Feynman diagrams: we need only consider the connected 
Feynman diagrams, and don’t have to worry about the vacuum bubbles. Combining 
this with (3.95), we learn that the Green’s functions G^(x i... ,x n ) can be calculated 
by summing over all connected Feynman diagrams, 

(fl| T i) ... <j> H (x n ) |0) = E Connected Feynman Graphs (3.104) 

An Example: The Four-Point Correlator: (fi| T(j) H (x i) ... 0# (aq) |fl) 

As a simple example, let’s look at the four-point correlation function in </> 4 theory. The 
sum of connected Feynman diagrams is given by, 

X ' q _ q X 2 Xl o / 2 

+ 2 Similar + + 

x 3 ° °x 4 x 3 ° ^>x 4 


X| 


X 2 



+ 5 Similar + 


All of these are connected diagrams, even though they don’t look 
that connected! The point is that a connected diagram is defined 
by the requirement that every line is joined to an external leg. An 
example of a diagram that is not connected is shown in the figure. 
As we have seen, such diagrams are taken care of in shifting the 
vacuum from |0) to 10). 



Figure 20: 


Feynman Rules 

The Feynman diagrams that we need to calculate for the Green’s functions depend on 
Xi ,... ,x n . This is rather different than the Feynman diagrams that we calculated for 
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the S-matrix elements, where we were working primarily with momentum eigenstates, 
and ended up integrating over all of space. However, it’s rather simple to adapt the 
Feynman rules that we had earlier in momentum space to compute G^ n >{x \..., x n ). 
For 0 4 theory, we have 

• Draw n external points x±,... ,x n , connected by the usual propagators and ver¬ 
tices. Assign a spacetime position y to the end of each line. 

• For each line x -y from xtoy write down a factor of the Feynman propagator 

AfO - y). 

• For each vertex at position y , write down a factor of — iX J d 4 y. 


3.7.2 From Green’s Functions to S-Matrices 

Having described how to compute correlation functions using Feynman diagrams, let’s 
now relate them back to the S-matrix elements that we already calculated. The first 
step is to perform the Fourier transform, 




{ | (l 4 x.< 


3 - tPi'Xi 


i —1 


G^( Xl ,...,x n ) 


(3.105) 


These are very closely related to the S-matrix elements that we’ve computed above. The 
difference is that the Feynman rules for G^ n \x i,..., x n ), effectively include propagators 
Ap for the external legs, as well as the internal legs. A related fact is that the 4- 
momenta assigned to the external legs is arbitrary: they are not on-shell. Both of these 
problems are easily remedied to allow us to return to the S-matrix elements: we need 
to simply cancel off the propagators on the external legs, and place their momentum 
back on shell. We have 


(p \, • • • ,Pn’\ s ~ 1 \pi ■ ■ ■ ,p n ) = {-i) n+n ' ]J(Pi 2 - m 2 ) 


m 


(3.106) 


i= 1 


3 =1 


x G fn+n '\-p' 1 ,...,-p' n „p 1 ,...,p n ) 


Each of the factors ( p 2 —m 2 ) vanishes once the momenta are placed on-shell. This means 
that we only get a non-zero answer for diagrams contributing to G ( ' n \x i,..., x n ) which 
have propagators for each external leg. 

So what’s the point of all of this? We’ve understood that ignoring the unconnected 
diagrams is related to shifting to the true vacuum O). But other than that, introducing 
the Green’s functions seems like a lot of bother for little reward. The important point 
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is that this provides a framework in which to deal with the true particle states in 
the interacting theory through renormalization. Indeed, the formula (3.106), suitably 
interpreted, remains true even in the interacting theory, taking into account the swarm 
of virtual particles surrounding asymptotic states. This is the correct way to consider 
scattering. In this context, (3.106) is known as the LSZ reduction formula. You will 
derive it properly next term. 
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4. The Dirac Equation 

“A great deal more was hidden in the Dirac equation than the author had 
expected when he wrote it down in 1928. Dirac himself remarked in one of 
his talks that his equation was more intelligent than its author. It should 
be added, however, that it was Dirac who found most of the additional 
insights. Weisskopf on Dirac 

So far we’ve only discussed scalar fields such that under a Lorentz transformation 
x M —y (x') M = A^oA, the field transforms as 

4>(x) —> <t>'{x) = (f){A^ l x ) (4.1) 

We have seen that quantization of such fields gives rise to spin 0 particles. But most 
particles in Nature have an intrinsic angular momentum, or spin. These arise naturally 
in field theory by considering fields which themselves transform non-trivially under the 
Lorentz group. In this section we will describe the Dirac equation, whose quantization 
gives rise to fermionic spin 1/2 particles. To motivate the Dirac equation, we will start 
by studying the appropriate representation of the Lorentz group. 

A familiar example of a field which transforms non-trivially under the Lorentz group 
is the vector field A^{x) of electromagnetism, 

A^x) ->■ M x v A v {k~ 1 x) (4.2) 

We’ll deal with this in Section 6. (It comes with its own problems!). In general, a field 
can transform as 

r(x) -> D[A]“ t y(A-'i) (4.3) 

where the matrices D[ A] form a representation of the Lorentz group, meaning that 

D[Ax]D[A 2 ] = £>[A!A 2 ] (4.4) 

and D[A _1 ] = D[A] _1 and _D[1] = 1. How do we find the different representations? 
Typically, we look at infinitesimal transformations of the Lorentz group and study the 
resulting Lie algebra. If we write, 

K = 6» v + (4.5) 

for infinitesimal cu, then the condition for a Lorentz transformation A^A v p rf p = 
becomes the requirement that u ; is anti-symmetric: 

+ oj^ = 0 (4.6) 
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Note that an antisymmetric 4x4 matrix has 4 x 3/2 = 6 independent components, 
which agrees with the 6 transformations of the Lorentz group: 3 rotations and 3 boosts. 
It’s going to be useful to introduce a basis of these six 4x4 anti-symmetric matrices. 
We could call them with A = 1,..., 6. But in fact it’s better for us (although 

initially a little confusing) to replace the single index A with a pair of antisymmetric 
indices [pa], where p, a = 0,..., 3, so we call our matrices (Ai p<T ) p ' v . The antisymmetry 
on the p and a indices means that, for example, AA m = —Af 10 , etc, so that p and a 
again label six different matrices. Of course, the matrices are also antisymmetric on 
the pz/ indices because they are, after all, antisymmetric matrices. With this notation 
in place, we can write a basis of six 4x4 antisymmetric matrices as 

(M p<t Y v = r] pp rf v - rf p rf v (4.7) 

where the indices p and v are those of the 4x4 matrix, while p and a denote which 
basis element we’re dealing with. If we use these matrices for anything practical (for 
example, if we want to multiply them together, or act on some field) we will typically 
need to lower one index, so we have 

(M pa )Y = rj pp 8\ - r) ap 8 p u (4.8) 


Since we lowered the index with the Minkowski metric, we pick up various minus signs 
which means that when written in this form, the matrices are no longer necessarily 
antisymmetric. Two examples of these basis matrices are, 


(M°y u 



and (A4 12 K 



(4.9) 


The first, Af 01 , generates boosts in the x l direction. It is real and symmetric. The 
second, AI 12 , generates rotations in the (x 1 , x 2 )-plane. It is real and antisymmetric. 
We can now write any u^ u as a linear combination of the Af pcr , 

u,Y = ln p<J (M p T„ (4.io) 


where Vt pa are just six numbers (again antisymmetric in the indices) that tell us what 
Lorentz transformation we’re doing. The six basis matrices A4 pa are called the gen¬ 
erators of the Lorentz transformations. The generators obey the Lorentz Lie algebra 
relations, 


[M p Y M tu ] = g aT M pv - r) pT M° v + r) pv M ar - if v M pT (4.11) 
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where we have suppressed the matrix indices. A finite Lorentz transformation can then 
be expressed as the exponential 

A = exp (| Q p(J A4 pa ) (4-12) 

Let me stress again what each of these objects are: the Ai pa are six 4x4 basis elements 
of the Lorentz Lie algebra; the Cl pa are six numbers telling us what kind of Lorentz 
transformation we’re doing (for example, they say things like rotate by 9 = n/7 about 
the x 3 -direction and run at speed v — 0.2 in the x l direction). 

4.1 The Spinor Representation 

We’re interested in finding other matrices which satisfy the Lorentz algebra commuta¬ 
tion relations (4.11). We will construct the spinor representation. To do this, we start 
by defining something which, at first sight, has nothing to do with the Lorentz group. 
It is the Clifford algebra , 


{ 7 ^, Y} = YY + YY = 2t/ m, 1 (4.13) 

where 7 ^, with /a — 0,1, 2, 3, are a set of four matrices and the l on the right-hand side 
denotes the unit matrix. This means that we must find four matrices such that 


YY = ~YY when fj, v 


(4.14) 


and 


( 7 °) 2 = 1 , Yf = -1 * = 1,2,3 (4.15) 


It’s not hard to convince yourself that there are no representations of the Clifford 
algebra using 2 x 2 or 3 x 3 matrices. The simplest representation of the Clifford 
algebra is in terms of 4 x 4 matrices. There are many such examples of 4 x 4 matrices 
which obey (4.13). For example, we may take 


7 


0 





(4.16) 


where each element is itself a 2 x 2 matrix, with the a 1 the Pauli matrices 


cr 


1 





which themselves satisfy {a l ,Y} = 2 8 l f 



(4.17) 
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One can construct many other representations of the Clifford algebra by taking 
Py'W -1 for any invertible matrix V. However, up to this equivalence, it turns out 
that there is a unique irreducible representation of the Clifford algebra. The matrices 
(4.16) provide one example, known as the Weyl or chiral representation (for reasons 
that will soon become clear). We will soon restrict ourselves further, and consider only 
representations of the Clifford algebra that are related to the chiral representation by 
a unitary transformation V. 

So what does the Clifford algebra have to do with the Lorentz group? Consider the 
commutator of two y M , 



(4.18) 


Let’s see what properties these matrices have: 


Claim 4.1: [S pu , y p ] = — 7 v r\^ 1 


Proof: When p 7 ^ v we have 


\s^,Y] = \b“Y,Y\ 



{Y, Y }Y + \YY~t 


1 M, p V 


□ 


Claim 4.2: The matrices S pu form a representation of the Lorentz algebra (4.11), 
meaning 


ifiu gpo-j _ _ fjVPgvv _|_ g ,u p — 'q l " r S fip 


(4.19) 


Proof: Taking p 7 ^ < 7 , and using Claim 4.1 above, we have 



(4.20) 


Now using the expression (4.18) to write y^y 0 " = 2 S ptT + y MfT , we have 


[S'*", S P<T ] = S pr7 p up - S va p pp + S pp p U(T - S pv p ap 


(4.21) 


which is our desired expression. 


□ 
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4.1.1 Spinors 

The S pM are 4x4 matrices, because the 7 ^ are 4x4 matrices. So far we haven’t given 
an index name to the rows and columns of these matrices: we’re going to call them 
a, (3 = 1,2, 3,4. 

We need a field for the matrices ( S v ' l, ) a a to act upon. We introduce the Dirac spinor 
field ij; a (x), an object with four complex components labelled by a = 1,2, 3,4. Under 
Lorentz transformations, we have 

i/> a (x) -> S[A] a p ^(\~ l x) (4.22) 

where 

A = exp (1 fl pa J\4 pa ) (4.23) 

S'[A] = exp (1 Q pa S pa ) (4.24) 

Although the basis of generators A4 pa and S pa are different, we use the same six 
numbers Vt pa in both A and S[A]: this ensures that we’re doing the same Lorentz 
transformation on x and i/j. Note that we denote both the generator S pa and the full 
Lorentz transformation S[A] as “S”. To avoid confusion, the latter will always come 
with the square brackets [A], 

Both A and S[A] are 4x4 matrices. So how can we be sure that the spinor repre¬ 
sentation is something new, and isn’t equivalent to the familiar representation A^? To 
see that the two representations are truly different, let’s look at some specific transfor¬ 
mations. 

Rotations 


S ij 


if 0 / 0 cr j \ 

2 ^ —cd 0 ) \ -a j 0 ) 



(for i ^ j ) (4.25) 


If we write the rotation parameters as = —£i]kp k (meaning D 12 = — <p 3 , etc) then 
the rotation matrix becomes 

S[A] = exp (1 Sl„S<") = ( e+ ‘" /2 (4.26) 

where we need to remember that D 12 = — f^i = — <p 3 when following factors of 2 . 
Consider now a rotation by 2-k about, say, the x 3 -axis. This is achieved by <p = (0, 0, 2n), 
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and the spinor rotation matrix becomes, 


‘S'[A] 

Therefore under a 2tc rotation 



-1 


(4.27) 


il> a (x) ->• -ip a (x) (4.28) 

which is definitely not what happens to a vector! To check that we haven’t been 
cheating with factors of 2, let’s see how a vector would transform under a rotation by 
<p = (0, 0, <p 3 ). We have 

( 00 0 0 \ 

0 V o 3 0 ) ( 4 - 29 ) 

0 0 0 0 / 

So when we rotate a vector by p 3 = 2i r, we learn that A = 1 as you would expect. So 
S'[A] is dehnitely a different representation from the familiar vector representation A^ u . 


Boosts 


S« = ±( 01 \( 0 o\ 

2 ^ 1 0 / \ -cA 0 / 2^0 CT* / 

Writing the boost parameter as fho = — Ho* = Xi, we have 


‘S'[A] 


e + x-^/2 o \ 

0 J 


(4.30) 


(4.31) 


Representations of the Lorentz Group are not Unitary 

Note that for rotations given in (4.26), S[A] is unitary, satisfying S[A]tS[A] = 1. But 
for boosts given in (4.31), S[A] is not unitary. In fact, there are no finite dimensional 
unitary representations of the Lorentz group. We have demonstrated this explicitly for 
the spinor representation using the chiral representation (4.16) of the Clifford algebra. 
We can get a feel for why it is true for a spinor representation constructed from any 
representation of the Clifford algebra. Recall that 


‘S'[A] = exp (1 Q pa S pa ) (4.32) 

so the representation is unitary if S are anti-hermitian, i.e. (S 7 ^/ = —S^ u . But we 
have 



which can be anti-hcrmitian if all 7 ^ are hermitian or all are anti-hermitian. However, 
we can never arrange for this to happen since 

(7 0 ) 2 = 1 =>• Real Eigenvalues 

(Y) 2 = —1 => Imaginary Eigenvalues (4.34) 

So we could pick 7 0 to be hermitian, but we can only pick 7* to be anti-hermitian. 
Indeed, in the chiral representation (4.16), the matrices have this property: (7°)^ = 7 0 
and (7*)f = —7*. In general there is no way to pick 7^ such that S ^ are anti-hermitian. 

4.2 Constructing an Action 

We now have a new held to work with, the Dirac spinor i/j. We would like to construct 
a Lorentz invariant equation of motion. We do this by constructing a Lorentz invariant 
action. 

We will start in a naive way which won’t work, but will give us a clue how to proceed. 
Define 


= (0*) T (x) (4.35) 

which is the usual adjoint of a multi-component object. We could then try to form a 
Lorentz scalar by taking the product 040, with the spinor indices summed over. Let’s 
see how this transforms under Lorentz transformations, 

0(x) —» S'[A] 0(A -1 £) 

0 1 '(x) —* ^(A -1 #) S[A]1’ (4.36) 

So 0t(x)0(x) —> ^(A -1 x)S[A]1'S[A]'?/;(A -1 x). But, as we have seen, for some Lorentz 
transformation S[A]tS[A] Y 1 since the representation is not unitary. This means that 
0'0 isn’t going to do it for us: it doesn’t have any nice transformation under the 
Lorentz group, and certainly isn’t a scalar. But now we see why it fails, we can also see 
how to proceed. Let’s pick a representation of the Clifford algebra which, like the chiral 
representation (4.16), satisfies ( 7 °)^ = 7 0 and ( 7 *)^ = — 7 *. Then for all /i — 0,1,2,3 
we have 


7V7 0 = (7“) t 


(4.37) 


which, in turn, means that 

(sn* = j[(7") t . (7“)*] = -7°S'“'7° (4.38) 
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so that 

S[A]t = exp (i = 7° S[A]-V (4-39) 

With this in mind, we now dehne the Dirac adjoint 

ip(x) — ip\x) 7 0 (4.40) 

Let’s now see what Lorentz covariant objects we can form out of a Dirac spinor ip and 
its adjoint ip. 

Claim 4.3: ipip is a Lorentz scalar. 

Proof: Linder a Lorentz transformation, 

ip(x)ip(x) = ip'(x) 7 0 ip(x) 

ip^A^x) SjAj^y 0 S[A]-0(A _1 a;) 

= ip\A~ l x)^° ip{A~ l x) 

= ip(A~ 1 x)ip(A~ 1 x) (4-41) 

which is indeed the transformation law for a Lorentz scalar. □ 

Claim 4.4: ip^^ip is a Lorentz vector, which means that 

ip(x) 7 ^ ip(x) —> A p v ip(A~ l x) 7 " ip(A~ l x) (4.42) 

This equation means that we can treat the /j, — 0,1,2,3 index on the 7 Al matrices as 

a true vector index. In particular we can form Lorentz scalars by contracting it with 
other Lorentz indices. 

Proof: Suppressing the x argument, under a Lorentz transformation we have, 


iPYip ^^pS[A\~ l Y S[A\iP (4.43) 

If ip 7 m ip is to transform as a vector, we must have 

S'[A ] -1 Y S'[A] = A^Y (4.44) 

We’ll now show this. We work infinitesimally, so that 

A = exp (i n pa M pa ) ~ 1 + \D pa M pa + ... (4.45) 

S'[A] = exp (4 n pa S p ") « 1 + \D pa S p ° + ... (4.46) 
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so the requirement (4.44) becomes 


-[S pa ,^} = (M pry y v Y (4.47) 

where we’ve suppressed the cy/3 indices on 7 M and S pu , but otherwise left all other 
indices explicit. In fact equation (4.47) follows from Claim 4.1 where we showed that 
[S pcr , 7 m ] = 7 p r] ap — 7 a r] pjp . To see this, we write the right-hand side of (4.47) by 
expanding out A4, 


(M p °Tvl v 


(v p,1 K - V ° p 5 p )f 

T] pp 7 ff - rf V 


which means that the proof follows if we can show 

-[S pa , Y] =ri pp ^ -rf P Y 
which is exactly what we proved in Claim 4.1. 


(4.48) 


(4.49) 

□ 


Claim 4.5: 07^7^0 transforms as a Lorentz tensor. More precisely, the symmet¬ 
ric part is a Lorentz scalar, proportional to rf‘ v 00, while the antisymmetric part is a 
Lorentz tensor, proportional to 05 7 U ' / 0. 


Proof: As above. 


□ 


We are now armed with three bilinears of the Dirac field, -0-0, 0y M 0 and 07^700 
each of which transforms covariantly under the Lorentz group. We can try to build a 
Lorentz invariant action from these. In fact, we need only the first two. We choose 

S — J d 4 x jp(x) 07 m c0 — m) 0 (x) (4.50) 

This is the Dirac action. The factor of “i” is there to make the action real; upon complex 
conjugation, it cancels a minus sign that comes from integration by parts. (Said another 
way, it’s there for the same reason that the Hermitian momentum operator — iV in 
quantum mechanics has a factor i). As we will see in the next section, after quantization 
this theory describes particles and anti-particles of mass \m\ and spin 1/2. Notice that 
the Lagrangian is first order, rather than the second order Lagrangians we were working 
with for scalar fields. Also, the mass appears in the Lagrangian as m, which can be 
positive or negative. 
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4.3 The Dirac Equation 


The equation of motion follows from the action (4.50) by varying with respect to ^ and 
independently. Varying with respect to "0, we have 

(iy^d^ — m)i/j — 0 (4.51) 

This is the Dirac equation. It’s completely gorgeous. Varying with respect to 0 gives 
the conjugate equation 


id 00 Y + m0 = 0 (4.52) 

The Dirac equation is first order in derivatives, yet miraculously Lorentz invariant. If 
we tried to write down a first order equation of motion for a scalar field, it would look 
like = ..., which necessarily includes a privileged vector in spacetime and is 

not Lorentz invariant. However, for spinor fields, the magic of the matrices means 
that the Dirac Lagrangian is Lorentz invariant. 

The Dirac equation mixes up different components of 0 through the matrices 7V 
However, each individual component itself solves the Klcin-Gordon equation. To see 
this, write 


(iy u d u + m^iy^dn — tti) 0 = — + m2 ) 0 = 0 (4.53) 

But YY^Yv = \{YiY}d»d v = d so we get 

— ( d^ + m 2 )0 = 0 (4.54) 

where this last equation has no 7^ matrices, and so applies to each component 0“, with 
a = 1,2, 3,4. 

The Slash 

Let’s introduce some useful notation. We will often come across 4-vectors contracted 
with Y matrices. We write 


A^Y = A 


(4.55) 


so the Dirac equation reads 


(i0 — m)ip = 0 


(4.56) 
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4.4 Chiral Spinors 

When we’ve needed an explicit form of the 7^ matrices, we’ve used the chiral represen¬ 
tation 


7 


0 




(4.57) 


In this representation, the spinor rotation transformation 5[A rot ] and boost transfor¬ 
mation AfAboost] were computed in (4.26) and (4.31). Both are block diagonal, 


S[A rot ] = 


o +i0-ff/2 


0 


q e +i(p-9/2 


and S[A boost ] 


/ e +xW2 0 \ 

\ 0 e -X-3/2 J 


(4.58) 


This means that the Dirac spinor representation of the Lorentz group is reducible. It 
decomposes into two irreducible representations, acting only on two-component spinors 
u± which, in the chiral representation, are defined by 


= 



(4.59) 


The two-component objects u± are called Weyl spinors or chiral spinors. They trans¬ 
form in the same way under rotations, 

u± —>■ e*^'^ /2 u± (4.60) 

but oppositely under boosts, 

u± ->■ e ± * s/2 u± (4.61) 

In group theory language, u + is in the (|,0) representation of the Lorentz group, 
while U- is in the (0, |) representation. The Dirac spinor ^ lies in the (|,0) © (0, |) 
representation. (Strictly speaking, the spinor is a representation of the double cover of 
the Lorentz group SL( 2, C)). 


4.4.1 The Weyl Equation 

Let’s see what becomes of the Dirac Lagrangia.n under the decomposition (4.59) into 
Weyl spinors. We have 

C = i[){i 0 — m )^ = iuLa^d^u^ + iu\_o- ti d fl u + — m(u' + u^ + u^_u + ) = 0 (4.62) 
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where we have introduced some new notation for the Pauli matrices with a pc — 0,1, 2, 3 
index, 

cd 1 = (1 ,a l ) and = (1, — a 1 ) (4.63) 

From (4.62), we see that a massive fermion requires both u + and U-, since they couple 
through the mass term. However, a massless fermion can be described by u + (or uJ) 
alone, with the equation of motion 


ia^d^u + = 0 

or ia^d^U- = 0 (4.64) 

These are the Weyl equations. 

Degrees of Freedom 

Let me comment here on the degrees of freedom in a spinor. The Dirac fermion has 
4 complex components = 8 real components. How do we count degrees of freedom? 
In classical mechanics, the number of degrees of freedom of a system is equal to the 
dimension of the configuration space or, equivalently, half the dimension of the phase 
space. In held theory we have an infinite number of degrees of freedom, but it makes 
sense to count the number of degrees of freedom per spatial point: this should at least 
be finite. For example, in this sense a real scalar held 0 has a single degree of freedom. 
At the quantum level, this translates to the fact that it gives rise to a single type of 
particle. A classical complex scalar held has two degrees of freedom, corresponding to 
the particle and the anti-particle in the quantum theory. 

But what about a Dirac spinor? One might think that there are 8 degrees of freedom. 
But this isn’t right. Crucially, and in contrast to the scalar held, the equation of motion 
is hrst order rather than second order. In particular, for the Dirac Lagrangian, the 
momentum conjugate to the spinor 0 is given by 

7r^, = d£/dip = iift (4.65) 

It is not proportional to the time derivative of 0. This means that the phase space for 
a spinor is therefore parameterized by 0 and 0', while for a scalar it is parameterized 
by 0 and n = 0. So the phase space of the Dirac spinor 0 has 8 real dimensions and 
correspondingly the number of real degrees of freedom is 4. We will see in the next 
section that, in the quantum theory, this counting manifests itself as two degrees of 
freedom (spin up and down) for the particle, and a further two for the anti-particle. 

A similar counting for the Weyl fermion tells us that it has two degrees of freedom. 
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4.4.2 7 5 

The Lorentz group matrices 5[A] came out to be block diagonal in (4.58) because we 
chose the specific representation (4.57). In fact, this is why the representation (4.57) 
is called the chiral representation: it’s because the decomposition of the Dirac spinor 
0 is simply given by (4.59). But what happens if we choose a different representation 
7^ of the Clifford algebra, so that 


7 m —» U^U and 0 —» 1/0 ? 


(4.66) 


Now S'[A] will not be block diagonal. Is there an invariant way to define chiral spinors? 
We can do this by introducing the “fifth” gamma-matrix 


7 5 = —i7°7 1 7 2 7 3 


(4.67) 


You can check that this matrix satisfies 


{7 5 )7 m } = 0 and (7 5 ) 2 = +1 (4.68) 

The reason that this is called y 5 is because the set of matrices y A = (y^iy 5 ), with 
A = 0,1, 2, 3,4 satisfy the five-dimensional Clifford algebra {y^y 5 } = 2rj AB . (You 
might think that y 4 would be a better name! But y 5 is the one everyone chooses - it’s 
a more sensible name in Euclidean space, where A = 1,2, 3,4, 5). You can also check 
that [S^y 5 ] = 0, which means that y 5 is a scalar under rotations and boosts. Since 
(7 5 ) 2 = 1) this means we may form the Lorentz invariant projection operators 

P± = j(l ± 7 5 ) (4.69) 

such that P 2 = P + and P 2 = P_ and P+P~ = 0. One can check that for the chiral 
representation (4.57), 


7 


5 



(4.70) 


from which we see that the operators P± project onto the Weyl spinors u±. However, 
for an arbitrary representation of the Clifford algebra, we may use y 5 to define the 
chiral spinors, 


0 ± = P±^ 


(4.71) 


which form the irreducible representations of the Lorentz group. 0 + is often called a 
“left-handed” spinor, while 0_ is “right-handed”. The name comes from the way the 
spin precesses as a massless fermion moves: we’ll see this in Section 4.7.2. 
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4.4.3 Parity 

The spinors are related to each other by parity. Let’s pause to define this concept. 
The Lorentz group is defined by x p —> A^ v x v such that 

A^A p a rf ff = rf p (4.72) 

So far we have only considered transformations A which are continuously connected to 
the identity; these are the ones which have an infinitesimal form. However there are 
also two discrete symmetries which are part of the Lorentz group. They are 

Time Reversal T : x° —» — x° ; x l —* x l 

Parity P : x° —> x° ; x l —* —x l (4.73) 

We won’t discuss time reversal too much in this course. (It turns out to be represented 
by an anti-unitary transformation on states. See, for example the book by Peskin and 
Schroeder). But parity has an important role to play in the standard model and, in 
particular, the theory of the weak interaction. 

Under parity, the left and right-handed spinors are exchanged. This follows from the 
transformation of the spinors under the Lorentz group. In the chiral representation, we 
saw that the rotation (4.60) and boost (4.61) transformations for the Weyl spinors u± 
are 

U± ^ e i0.*/2 U± and u± Wt e ±x-ff/2 u± (4.74) 

Under parity, rotations don’t change sign. But boosts do flip sign. This confirms that 
parity exchanges right-handed and left-handed spinors, P : u± —> u T , or in the notation 
ip± — |(1 ± 7 5 )V , j we have 

P : Mx,t) ->■ (4.75) 

Using this knowledge of how chiral spinors transform, and the fact that P 2 = 1, we see 
that the action of parity on the Dirac spinor itself can be written as 

P : i/j(x, t) —> ^ 0r ip(—x,t) (4.76) 

Notice that if i^{x^t) satisfies the Dirac equation, then the parity transformed spinor 
7 O, 0(— x, t) also satisfies the Dirac equation, meaning 

(i^f 0 d t + i'fdi — m)7°'0(— x, t ) = 7°(i7°<9t — i^ l di — x, t) = 0 (4.77) 

where the extra minus sign from passing 7 0 through 7* is compensated by the derivative 
acting on — x instead of +x. 


94 


4.4.4 Chiral Interactions 


Let’s now look at how our interaction terms change under parity. We can look at each 
of our spinor bilinears from which we built the action, 

P : ipil;(x,t)ipip(—x,t) (4.78) 

which is the transformation of a scalar. For the vector yb 7 +/?, we can look at the 
temporal and spatial components separately, 

P : ip”/ 0 ip(x,t) — y 7/7 V(—a?, f) 

P : 'ip'y l 'ip(x,t) —* '07°7*7 O ^(—a:, t) = — 7/7 (4.79) 

which tells us that ^ 7 transforms as a vector, with the spatial part changing sign. 
You can also check that transforms as a suitable tensor. 

However, now we’ve discovered the existence of 7 s , we can form another Lorentz 
scalar and another Lorentz vector, 

'07 5, 0 and , 07 5 7 /i ^ (4.80) 


How do these transform under parity? We can check: 


P : t) —» '07°7 5 7° ip(—x,t) = — 

P : '07 5 7 M ^(a:, t) —* ^7°7 5 7 M 7 °^(—x,t) = 


^7 5 ^(-f, t) 
f — '07 5 7°a/; (— df, t) 
I +'07 5 7*^(— x, t) 


(4.81) 

H — 0 

H = i 


which means that transforms as a pseudoscalar, while transforms as an 

axial vector. To summarize, we have the following spinor bilinears, 



scalar 

'07 M, 0 

vector 


tensor 

'ip'y 5 ip 

pseudoscalar 

'07 5 7 m '0 

axial vector 


(4.82) 


The total number of bilinears is 1 + 4 +(4x3/2)+ 4 + 1 = 16 which is all we could 
hope for from a 4-component object. 
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We’re now armed with new terms involving y 5 that we can start to add to our 
Lagrangian to construct new theories. Typically such terms will break parity invariance 
of the theory, although this is not always true. (For example, the term 00y 5 0 doesn’t 
break parity if 0 is itself a pseudoscalar). Nature makes use of these parity violating 
interactions by using y 5 in the weak force. A theory which treats ip± on an equal 
footing is called a vector-like theory. A theory in which 0 + and 0_ appear differently 
is called a chiral theory. 

4.5 Majorana Fermions 

Our spinor 0 Q is a complex object. It has to be because the representation S'[A] 
is typically also complex. This means that if we were to try to make 0 real, for 
example by imposing 0 = 0*, then it wouldn’t stay that way once we make a Lorentz 
transformation. However, there is a way to impose a reality condition on the Dirac 
spinor 0. To motivate this possibility, it’s simplest to look at a novel basis for the 
Clifford algebra, known as the Majorana basis. 



These matrices satisfy the Clifford algebra. What is special about them is that they 
are all pure imaginary (y^)* = —yC This means that the generators of the Lorentz 
group S ,IV = 100, 7 "], and hence the matrices S'[A] are real. So with this basis of the 
Clifford algebra, we can work with a real spinor simply by imposing the condition, 

0 = 0* (4.83) 

which is preserved under Lorentz transformation. Such spinors are called Majorana 
spinors. 

So what’s the story if we use a general basis for the Clifford algebra? We’ll ask only 
that the basis satisfies (y°0 = 7 0 and (y*0 = — 7 *. We then define the charge conjugate 
of a Dirac spinor 0 as 

0 (c) = C0* (4.84) 

Here C is a 4 x 4 matrix satisfying 

&C = 1 and CVC = -( 7 **)* (4.85) 

Let’s firstly check that (4.84) is a good definition, meaning that 0^ transforms nicely 
under a Lorentz transformation. We have 

0 (c) C'S'[A]*0* = A[A]C0* = S'[A]0^ (4.86) 
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where we’ve made use of the properties (4.85) in taking the matrix C through 5[A]*. 
In fact, not only does ip^ transform nicely under the Lorentz group, but if ip satisfies 
the Dirac equation, then ip( c > does too. This follows from, 

(i 0 — m)ip = 0 =>• (—i 0 * — m)ip* = 0 

=>- C(—i 0* — m)ip* = (+i 0 — m)ip^ = 0 

Finally, we can now impose the Lorentz invariant reality condition on the Dirac spinor, 
to yield a Majorana spinor, 

ip {c) = Ip (4.87) 

After quantization, the Majorana spinor gives rise to a fermion that is its own anti¬ 
particle. This is exactly the same as in the case of scalar fields, where we’ve seen that 
a real scalar held gives rise to a spin 0 boson that is its own anti-particle. (Be aware: 
In many texts an extra factor of 7 0 is absorbed into the definition of C ). 


So what is this matrix C? Well, for a given representation of the Clifford algebra, it 
is something that we can find fairly easily. In the Majorana basis, where the gamma 
matrices are pure imaginary, we have simply CMaj = 1 and the Majorana condition 
ip = np( c ) becomes ip = ip*. In the chiral basis (4.16), only y 2 is imaginary, and we 
may take C Ai rai = iy 2 = ^ _° ct2 “ j. (The matrix ia 2 that appears here is simply the 

anti-symmetric matrix e Q/3 ). It is interesting to see how the Majorana condition (4.87) 
looks in terms of the decomposition into left and right handed Weyl spinors (4.59). 
Plugging in the various definitions, we find that u+ = ia 2 u*_ and U- = — ia 2 u* + . In 
other words, a Majorana spinor can be written in terms of Weyl spinors as 


ip = 


u+ 

—i<j 2 u\ 


(4.88) 


Notice that it’s not possible to impose the Majorana condition ip = ip^ c ' at the same 
time as the Weyl condition (u_ = 0 or u + = 0). Instead the Majorana condition relates 
U- and u + . 


An Aside: Spinors in Different Dimensions: The ability to impose Majorana 
or Weyl conditions on Dirac spinors depends on both the dimension and the signature of 
spacetime. One can always impose the Weyl condition on a spinor in even dimensional 
Minkowski space, basically because you can always build a suitable “y 5 ” projection 
matrix by multiplying together all the other y-matrices. The pattern for when the 
Majorana condition can be imposed is a little more sporadic. Interestingly, although 
the Majorana condition and Weyl condition cannot be imposed simultaneously in four 
dimensions, you can do this in Minowski spacetimes of dimension 2,10,18,.... 
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4.6 Symmetries and Conserved Currents 

The Dirac Lagrangian enjoys a number of symmetries. Here we list them and compute 
the associated conserved currents. 

Spacetime Translations 

Under spacetime translations the spinor transforms as 

Stjj = e p dY (4.89) 

The Lagrangian depends on d^, but not d^, so the standard formula (1.41) gives us 
the energy-momentum tensor 


T pv = i^Yd v 'ip - rf v L (4.90) 

Since a current is conserved only when the equations of motion are obeyed, we don’t lose 
anything by imposing the equations of motion already on T pv . In the case of a scalar 
field this didn’t really buy us anything because the equations of motion are second 
order in derivatives, while the energy-momentum is typically first order. However, for 
a spinor field the equations of motion are first order: (i 0 — m)i(j = 0. This means we 
can set £, = 0 in T pu , leaving 


T pu = (4.91) 

In particular, we have the total energy 

E = / A r«> = / A im = / *z + ■»)* (4.92) 

where, in the last equality, we have again used the equations of motion. 

Lorentz Transformations 

Under an infinitesimal Lorentz transformation, the Dirac spinor transforms as (4.22) 
which, in infinitesimal form, reads 

SY = -uYx u dY a + \tt pa (S pa YY P (4.93) 

where, following (4.10), we have u = |fi p<T (A / t pa )Y, and A4 P<7 are the generators of 

the Lorentz algebra given by (4.8) 

(M pa Y = V™ K ~ V afl K ( 4 -94) 
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which, after direct substitution, tells us that U) pu = Q pv . So we get 


sr = [x v a^ a - us^r^] 


(4.95) 


The conserved current arising from Lorentz transformations now follows from the same 
calculation we saw for the scalar held (1.54) with two differences: firstly, as we saw 
above, the spinor equations of motion set C = 0; secondly, we pick up an extra piece 
in the current from the second term in (4.95). We have 


(J p ) pa = x p T pa - x a T pp - ii])^S pa ip 


(4.96) 


After quantization, when { J {J ) pa is turned into an operator, this extra term will be 
responsible for providing the single particle states with internal angular momentum, 
telling us that the quantization of a Dirac spinor gives rise to a particle carrying spin 


1 / 2 . 


Internal Vector Symmetry 

The Dirac Lagrangian is invariant under rotating the phase of the spinor, 0 —>■ e _ *“0. 
This gives rise to the current 


3v = VD'V 


(4.97) 


where “V” stands for vector , reflecting the fact that the left and right-handed compo¬ 


nents 0± transform in the same way under this symmetry. We can easily check that 
jy is conserved under the equations of motion, 


d^ijy = (<9 M 0) 7^0 + 07 /i (c00) = im 00 — irrupip = 0 


(4.98) 


where, in the last equality, we have used the equations of motion = rml> and 
id ll 'i]j r Y p ’ = —mi]). The conserved quantity arising from this symmetry is 



(4.99) 


We will see shortly that this has the interpretation of electric charge, or particle number, 
for fermions. 

Axial Symmetry 

When rn — 0, the Dirac Lagrangian admits an extra internal symmetry which rotates 
left and right-handed fermions in opposite directions, 


0 ->■ e ia ^i\) and 0 ->■ i])e ia ^ 


(4.100) 
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Here the second transformation follows from the first after noting that e 7 0 
7 °e +ia75 . This gives the conserved current, 


3 a = VdVV’ 


(4.101) 


where A is for “axial” since j 1 ^ is an axial vector. This is conserved only when m = 0. 
Indeed, with the full Dirac Lagrangian we may compute 


d^A = + ' 07 m 7 5 ^^ = 


(4.102) 


which vanishes only for m = 0. However, in the quantum theory things become more 
interesting for the axial current. When the theory is coupled to gauge fields (in a 
manner we will discuss in Section 6 ), the axial transformation remains a symmetry 
of the classical Lagrangian. But it doesn’t survive the quantization process. It is the 
archetypal example of an anomaly, a symmetry of the classical theory that is not 
preserved in the quantum theory. 

4.7 Plane Wave Solutions 

Let’s now study the solutions to the Dirac equation 


(27 — m ) i /> = 0 


(4.103) 


We start by making a simple ansatz: 


xp = u(p) e~ ip ' x 


(4.104) 


where u{p) is a four-component spinor, independent of spacetime x which, as the no¬ 
tation suggests, can depend on the 3-momentum p. The Dirac equation then becomes 

(Y'Pi, ~ m)u(p) = f ~ m ) u(p) = 0 (4.105) 

V ~m J 

where we’re again using the definition, 


( 1 , a 1 ) and = ( 1 , -<7*) 


(4.106) 


Claim: The solution to (4.105) is 



(4.107) 
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for any 2-component spinor £ which we will normalize to £t£ = 1. 


Proof: Let’s write uip) 1 = (ui,it 2 ). Then equation (4.105) reads 

(p ■ a)u 2 = mui and (p ■ a)ui = mw 2 


(4.108) 


Either one of these equations implies the other, a fact which follows from the identity 
(p • cr)(p ■ a) = p o — PiPj<r l c = Pq — PiPjS u = p^p p = m 2 . To start with, let’s try 
the ansatz u\ — (p ■ cr)£' for some spinor £'. Then the second equation in (4.108) 
immediately tells us that w 2 = m C- So we learn that any spinor of the form 


u(p) = A 


ip ■ a ) £' A 
m? J 


(4.109) 


with constant A is a solution to (4.105). To make this more symmetric, we choose 
A = 1/m and £' = sjp ■ with constant £. Then ui — (p ■ cr)^/p ■ = m^/p ■ cr£. So 

we get the promised result (4.107) □ 


Negative Frequency Solutions 

We get further solutions to the Dirac equation from the ansatz 

= v(p) e +ip ' x (4.110) 

Solutions of the form (4.104), which oscillate in time as ij> ~ e~ lEt , are called positive 
frequency solutions. If we compute the energy of these solutions using (4.92), we find 
that it is positive. Those of the form (4.110), which oscillate as ip ~ e +lEt , are negative 
frequency solutions. Now if we compute the energy using (4.92), it is negative. 

The Dirac equation requires that the 4-component spinor v(p) satisfies 

{YPu + m)v(p) = I m _ Pl * a ] v(p) = 0 (4.111) 

\p^a p m J 


which is solved by 


v{p) = 


-yjp-ap 


(4.112) 


for some 2-component spinor rj which we take to be constant and normalized to rpr) = 1. 
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4.7.1 Some Examples 

Consider the positive frequency solution with mass m and 3-momentum p = 0, 


u(p) = \fm 



(4.113) 


where £ is any 2-component spinor. Spatial rotations of the held act on £ by (4.26), 


£ e+^' ?/2 £ 


(4.114) 


The 2-conrponent spinor £ dehnes the spin of the held. This should be familiar from 
quantum mechanics. A held with spin up (down) along a given direction is described 
by the eigenvector of the corresponding Pauli matrix with eigenvalue +1 (-1 respec¬ 
tively). For example, £ 3 = (1,0) describes a held with spin up along the z-axis. After 
quantization, this will become the spin of the associated particle. In the rest of this 
section, we’ll indulge in an abuse of terminology and refer to the classical solutions to 
the Dirac equations as “particles”, even though they have no such interpretation before 
quantization. 


Consider now boosting the particle with spin £ 2 = (1, 0) along the x 3 direction, with 
= (E,0,0,p). The solution to the Dirac equation becomes 


u(p) = 



V^P 3 

ywry 



(4.115) 


In fact, this expression also makes sense for a massless held, for which E = p 3 . (We 
picked the normalization (4.107) for the solutions so that this would be the case). For 
a massless particle we have 


u(p) = V2~E 



(4.116) 


Similarly, for a boosted solution of the spin down £ r = (0,1) held, we have 


u(p) = 


yp 7 ^ 

Vp 7 ^ 



JEEP 3 




(4.117) 
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4.7.2 Helicity 

The helicity operator is the projection of the angular momentum along the direction 
of momentum, 


h 


\€ijkP 


ojk 



(cr 1 0 
\ 0 n* 


(4.118) 


where S'*- 7 is the rotation generator given in (4.25). The massless field with spin f 1 = 
(1,0) in (4.116) has helicity h = 1/2: we say that it is right-handed. Meanwhile, the 
field (4.117) has helicity h = —1/2: it is left-handed. 


4.7.3 Some Useful Formulae: Inner and Outer Products 

There are a number of identities that will be very useful in the following section, 
regarding the inner (and outer) products of the spinors u{p) and v(p). It’s firstly 
convenient to introduce a basis f and r] s , s = 1, 2 for the two-component spinors such 
that 


ft = an d r f^ v s = S rs (4.119) 

for example, 

and (4.120) 

and similarly for r] s . Let’s deal first with the positive frequency plane waves. The two 
independent solutions are now written as 


u s (p) 


Vp 7 ^^ ) 


(4.121) 


We can take the inner product of four-component spinors in two different ways: either 
as vf -u, or as u-u. Of course, only the latter will be Lorentz invariant, but it turns out 
that the former is needed when we come to quantize the theory. Here we state both: 


u r \p) ■u s (p) = (f f Vp 7 ^, f t v / WW) ( 

Wp 7 ^^ 

= fV af-hfVbf = 2fW = 


2p 0 S rs 


(4.122) 


while the Lorentz invariant inner product is 


u r {p) ■ u s (p) = (C 7 Vp 7 ^, C ] Vp^) 



Vp 7 ^^ ) 


2mS rs (4.123) 
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We have analogous results for the negative frequency solutions, which we may write as 


with v r ^(p) ■ v s (p) = 2p 0 8 rs 
and v r {p) ■ v s (p) = —2m8 rs 

We can also compute the inner product between u and v. We have 


v s (p) = 


y/p • 17 V s 
— \/v ■ crr] s 


(4.124) 


u 


\ p ) ■ v s {p) = (r f v^, c^Vp 7 ^) 7° 


Vp^v 3 


-sJp ■ a rf 

= fV(p ■ &)(p ■ a )v s - CW(p ■ °)(p ■ a )v s = o 


(4.125) 


and similarly, v r {p) ■ u s (p) = 0. However, when we come to iP ■ v, it is a slightly 
different combination that has nice properties (and this same combination appears 
when we quantize the theory). We look at u r \p) ■ v s (—p), with the 3-momentum in 
the spinor v taking the opposite sign. Defining the 4-momentum (p') IJ ‘ = (p°, —p), we 
have 


u r] (p) ■ v s (-p) = (f rt v/p 7 ^, C^Vp 7 *) 


\/p' ■ (J rf 
—a/p 7 • a rj s 

= cW(p- <?)& ■ a )v s - cW {p ■ p){p' ■ p)p s 


(4.126) 


Now the terms under the square-root are given by (p-cr)(p'-a) = (po+Picd)(p 0 —pj<j*) = 
Pq — p 2 = m 2 . The same expression holds for (p • a)(p' • cr), and the two terms cancel. 
We learn 


u r \p) ■ v s (—p) = v r \p) ■ u s (—p) = 0 


(4.127) 


Outer Products 

There’s one last spinor identity that we need before we turn to the quantum theory. It is: 


Claim: 


2 

u s (p)u s (p) — jj + m 

S=1 


(4.128) 


where the two spinors are not now contracted, but instead placed back to back to give 
a 4 x 4 matrix. Also, 


2 

v s (p)v s (p) — ]j — m 

S= 1 


(4.129) 
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Proof: 




V u s (p)u s (p) = 


S=1 


But V. f ’ f ’ 1 = 1. the 2 x 2 unit matrix, which then gives us 

2 

Y1 u s (p)u s (p) = 

s =1 

which is the desired result. A similar proof works for Y2 S v s (p)v s (p). 


m p ■ cr 
p ■ a m 
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(4.130) 

(4.131) 

□ 



5. Quantizing the Dirac Field 

We would now like to quantize the Dirac Lagrangian, 

C = '(p(x) (i 0 — m) ip(x) (5.1) 

We will proceed naively and treat -0 as we did the scalar held. But we’ll see that things 
go wrong and we will have to reconsider how to quantize this theory. 


5.1 A Glimpse at the Spin-Statistics Theorem 

We start in the usual way and define the momentum, 

Fir 

7r = —- = i'l/j'y 0 = iifi (5.2) 

For the Dirac Lagrangian, the momentum conjugate to 0 is i0h It does not involve 
the time derivative of %This is as it should be for an equation of motion that is first 
order in time, rather than second order. This is because we need only specify 0 and 
0 t on an initial time slice to determine the full evolution. 

To quantize the theory, we promote the held 0 and its momentum 0t to operators, 
satisfying the canonical commutation relations, which read 

[0«(z),0/j(y)] = [^ a (x),^p(y)] = o 

s (3 \x-y) (5.3) 

It’s this step that we’ll soon have to reconsider. 


Since we’re dealing with a free theory, where any classical solution is a sum of plane 
waves, we may write the quantum operators as 


V’(^) = X! 


d 3 p 


= Y 


-7 (2vr) 3 
d 3 p 1 


Tt J (2vr) 3 ^2E p 


btu s (p)e +i P s + c s J v s {p)e- if;i 
b s J u s {p)^e-^ g + d\v a {pj*e + * i 


(5.4) 


where the operators bp create particles associated to the spinors u s (p), while ct' create 
particles associated to v s (p). As with the scalars, the commutation relations of the 
helds imply commutation relations for the annihilation and creation operators 
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Claim: The field commutation relations (5.3) are equivalent to 

1 ^,61'] = (2 jt )Vi< s >(p — if) 

[cp,c‘f] = -(2 x) 3 6 r ’6 m (p-q) (5.5) 

with all other commutators vanishing. 

Note the strange minus sign in the [c, c^\ term. This means that we can’t define the 
ground state |0) as something annihilated by ch |0) = 0, because then the excited states 
c S p |0) would have negative norm. To avoid this, we will have to flip the interpretation 
of c and cd, with the vacuum defined by cl' 10) = 0 and the excited states by c5|0). 
This, as we will see, will be our undoing. 


Proof: Let’s show that the [6, b 1] and [c, c^] commutators reproduce the field com¬ 
mutators (5.3), 


(50 ,*'m = v 


E 


d 3 p d 3 q 1 
( 2 vr )6 

d 3 p 1 




J&p-y-o) 


(2t r) 3 2 E p 


+ [$, 

U s (p)u s (p) 7 °e^ M + u s (p)u s (p)7°e-^ M )(5.6) 


At this stage we use the outer product formulae (4.128) and (4.129) which tell us 
'£2 s uS (p)u s (p) = tf + m and ^2 s v s (p)v s (p) = — m, so that 

= I ^2t) 3 2g- ^ Po7 ° + Pi7 ' + m - l7 ° + * Po7 ° ~ Pi7 ’ ~~ m ^ 7 °l e+0>i! ^’ 1 

where, in the second term, we’ve changed p —> — p under the integration sign. Now, 
using p 0 = Ep we have 

W*),^®] = f = 5< 3 ) (2- 7) (5.7) 

as promised. Notice that it’s a little tricky in the middle there, making sure that 
the piY terms cancel. This was the reason we needed the minus sign in the [c, c^] 
commutator terms in (5.5). □ 
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5.1.1 The Hamiltonian 

To proceed, let’s construct the Hamiltonian for the theory. Using the momentum 
7T = we have 

W. = 7 rip — £ = + m)7p (5.8) 


which means that H = f d 3 x7-L agrees with the conserved energy computed using 
Noether’s theorem (4.92). We now wish to turn the Hamiltonian into an operator. 
Let’s firstly look at 




d 3 p 1 
(2tt)3 v /2U^ 


b S p{~l l Vi + m)u s (p) e +lp ' x + ci t ( 7 ®p, + m)v s (p) e lft ' x 


where, for once we’ve left the sum over s = 1, 2 implicit. There’s a small subtlety with 
the minus signs in deriving this equation that arises from the use of the Minkowski 
metric in contracting indices, so that p-x = x l p % = —x l pi. Now we use the defining 
equations for the spinors u s (p) and v s (p) given in (4.105) and (4.111), to replace 


(—7 l Pi + m)u s (p) = r y°PoU s (p) and {YPi + m)v s (p) = — 7 °Po' yS (p) (5.9) 


so we can write 

{-i'y i d i + m)if> = J ^ 

We now use this to write the operator Hamiltonian 



btu s (p) e +i ^ 




v s (p) e 


-ip-x 


(5.10) 


H 


J d 3 x'i/j^'y 0 (—iyd i + m)-0 


d 3 xd 3 pd 3 q / E, 


(2tt) 
d 3 p 1 

(27^2 -I'* 


AE, 


p b r Ju r (q)^e~^ s + c}v r (q)* e +i<r " 

btu s (p)e +i P- x - c s J v s (p) e~* 3 

fyfsWipj* ■ v s (p)} 


b r Mu r tf)'-v 


-b r ~cl[u r (p)^ ■ v s {-p )] + c r pb s _ p [v r {p)^ ■ u s (-p)} 

where, in the last two terms we have relabelled p —y —p. We now use our inner product 
formulae (4.122), (4.124) and (4.127) which read 

u r (p ) t • u s (p) = v r {jrf ■ v s (p) = 2 p 0 5 rs and u r (p)^ ■ v s (—p) = v r (p) Jf ■ u s (—p) = 0 
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giving us 


H 



(5.11) 


(5.12) 


The <5 (3) term is familiar and easily dealt with by normal ordering. The tfb term is 
familiar and we can check that b1 create positive energy states as expected, 



The minus sign in front of the c^c term should make us nervous. If we think of c' as 


creation operators then there’s no problem since, using the commutation relation (5.5), 
we still find that cd creates positive energy states, 



However, as we noted after (5.5), these states have negative norm. To have a sensible 
Hilbert space, we need to interpret c as the creation operator. But then the Hamiltonian 
is not bounded below because 



This is a disaster. Taken seriously it would tell us that we could tumble to states of 
lower and lower energy by continually producing c particles. As the English would say, 
it’s all gone a bit Pete Tong. (No relation). 

Since the above calculation was a little tricky, you might think that it’s possible to 
rescue the theory to get the minus signs to work out right. You can play around with 
different things, but you’ll always find this minus sign cropping up somewhere. And, 
in fact, it’s telling us something important that we missed. 

5.2 Fermionic Quantization 

The key piece of physics that we missed is that spin 1/2 particles are fermions, meaning 
that they obey Fermi-Dirac statistics with the quantum state picking up a minus sign 
upon the interchange of any two particles. This fact is embedded into the structure 
of relativistic quantum field theory: the spin-statistics theorem says that integer spin 
fields must be quantized as bosons, while half-integer spin fields must be quantized 
as fermions. Any attempt to do otherwise will lead to an inconsistency, such as the 
unbounded Hamiltonian we saw in (5.12). 
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So how do we go about quantizing a field as a fermion? Recall that when we quantized 
the scalar field, the resulting particles obeyed bosonic statistics because the creation 
and annihilation operators satisfied the commutation relations, 

[ 4,41 = ° =c <44|0> = |p,9> = |?,p> (5.13) 

To have states obeying fermionic statistics, we need anti-commutation relations, {A, B} = 
AB + BA. Rather than (5.3), we will ask that the spinor fields satisfy 

{>«(£), Vv»(£)} = = 0 

{^a(x),^p(y)} = S a pS^(x-y) (5.14) 

We still have the expansion (5.4) of and ijj' in terms of b,b',c and cfi But now the 
same proof that led us to (5.5) tells us that 

{ 63 , 6 *.'}= (2r) 3 S"i< 3 >(p-f) 

{c5,cl'}= (2T) 3 i r ‘^>(p-i~) (5.15) 

with all other anti-commutators vanishing, 

{63,6*-} = {4,4} = {63,4'} = {63,4-} = ... =0 (5.16) 


The calculation of the Hamiltonian proceeds as before, all the way through to the 
penultimate line (5.11). At that stage, we get 


H 


d 3 p 
(2?r) 3 
d 3 p 
(27r) 3 


Ef, 


b s J bp - c s pC s J 


Ei 5 


b s Jb%+r s ^ 


C -> C -> — 
'p "p ' p p 


(2ir) 3 6 <3) (0) 


(5.17) 


The anti-commutators have saved us from the indignity of an unbounded Hamiltonian. 
Note that when normal ordering the Hamiltonian we now throw away a negative contri¬ 
bution — (27r) 3 h^ 3 ' ) (0). In principle, this could partially cancel the positive contribution 
from bosonic fields. Cosmological constant problem anyone?! 


5.2.1 Fermi-Dirac Statistics 

Just as in the bosonic case, we define the vacuum |0) to satisfy, 

63|0) = 4-|0) = 0 (5.18) 
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Although b and c obey anti-commutation relations, the Hamiltonian (5.17) has nice 
commutation relations with them. You can check that 


\H,by l = -E fS b} and [ H,bO] = E p -bO 

[ff.Cp-] = and [H,cO] = E f cO (5.19) 

This means that we can again construct a tower of energy eigenstates by acting on the 
vacuum by b r J and cb* to create particles and antiparticles, just as in the bosonic case. 
For example, we have the one-particle states 

W,r) = b r J |0) (5.20) 


The two particle states now satisfy 

\pi,ri,p 2 ,r 2 ) = 6g + 6g f |0) = - \p 2 , r 2 ;pi, n) (5.21) 

confirming that the particles do indeed obey Fermi-Dirac statistics. In particular, we 
have the Pauli-Exclusion principle | p,r;p,r) = 0. Finally, if we wanted to be sure 
about the spin of the particle, we could act with the angular momentum operator 
(4.96) to confirm that a stationary particle \p = 0, r) does indeed carry intrinsic angular 
momentum 1/2 as expected. 

5.3 Dirac’s Hole Interpretation 

“In this attempt, the success seems to have been on the side of Dirac rather 
than logic” 

Pauli on Dirac 

Let’s pause our discussion to make a small historical detour. Dirac originally viewed 
his equation as a relativistic version of the Schrodinger equation, with ip interpreted 
as the wavefunction for a single particle with spin. To reinforce this interpretation, he 
wrote (i 0 — m)ip = 0 as 

d ip -> 

i— = —ia ■ Wip + mphp = Hpj (5.22) 

where a = — 7°7 and /3 — 7 0 . Here the operator H is interpreted as the one-particle 
Hamiltonian. This is a very different viewpoint from the one we now have, where ip 
is a classical field that should be quantized. In Dirac’s view, the Hamiltonian of the 
system is H defined above, while for us the Hamiltonian is the field operator (5.17). 
Let’s see where Dirac’s viewpoint leads. 


Ill 


With the interpretation of as a single-particle wavefunction, the plane-wave solu¬ 
tions (4.104) and (4.110) to the Dirac equation are thought of as energy eigenstates, 
with 

-0 = u (p) e~ ip ' x => = E p-0 

i/j = v(p) e +ip ' x => i d ^ = -E^ (5.23) 

which look like positive and negative energy solutions. The spectrum is once again 
unbounded below; there are states v(p) with arbitrarily low energy —Eg. At first 
glance this is disastrous, just like the unbounded held theory Hamiltonian (5.12). Dirac 
postulated an ingenious solution to this problem: since the electrons are fermions (a 
fact which is put in by hand to Dirac’s theory) they obey the Pauli-exclusion principle. 
So we could simply stipulate that in the true vacuum of the universe, all the negative 
energy states are filled. Only the positive energy states are accessible. These filled 
negative energy states are referred to as the Dirac sea. Although you might worry about 
the infinite negative charge of the vacuum, Dirac argued that only charge differences 
would be observable (a trick reminiscent of the normal ordering prescription we used 
for held operators). 

Having avoided disaster by boating on an infinite sea comprised of occupied negative 
energy states, Dirac realized that his theory made a shocking prediction. Suppose that 
a negative energy state is excited to a positive energy state, leaving behind a hole. 
The hole would have all the properties of the electron, except it would carry positive 
charge. After flirting with the idea that it may be the proton, Dirac finally concluded 
that the hole is a new particle: the positron. Moreover, when a positron comes across 
an electron, the two can annihilate. Dirac had predicted anti-matter, one of the greatest 
achievements of theoretical physics. It took only a couple of years before the positron 
was discovered experimentally in 1932. 

Although Dirac’s physical insight led him to the right answer, we now understand 
that the interpretation of the Dirac spinor as a single-particle wavefunction is not really 
correct. For example, Dirac’s argument for anti-matter relics crucially on the particles 
being fermions while, as we have seen already in this course, anti-particles exist for 
both fermions and bosons. What we really learn from Dirac’s analysis is that there is 
no consistent way to interpret the Dirac equation as describing a single particle. It is 
instead to be thought of as a classical held which has only positive energy solutions 
because the Hamiltonian (4.92) is positive definite. Quantization of this held then gives 
rise to both particle and anti-particle excitations. 
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This from Julian Schwinger: 


“Until now, everyone thought that the Dirac equation referred directly to 
physical particles. Now, in field theory, we recognize that the equations 
refer to a sublevel. Experimentally we are concerned with particles, yet the 
old equations describe fields.... When you begin with field equations, you 
operate on a level where the particles are not there from the start. It is 
when you solve the field equations that you see the emergence of particles. ” 


5.4 Propagators 

Let’s now move to the Heisenberg picture. We define the spinors if(x, t) at every point 
in spacetime such that they satisfy the operator equation 


We solve this by the expansion 


~di 




2 

x ) = X] 

s= 1 
2 

^\ x ) = 5Z 

S=1 


d 3 p 1 
(2tt)3 v /2E> L 

d 3 p 1 
(2tt)3 v /2 E> L' 


blu s (p)e~ ip - x + c s Jv s (p)e +ipc 


b s Ju s (p)^e +ipx + Av s {p)'e- ip " 


(5.24) 


(5.25) 


Let’s now look at the anti-commutators of these fields. We define the fermionic prop¬ 
agator to be 


iS af 3 = {^ a {x),iJjfi{y)} 


(5.26) 


In what follows we will often drop the indices and simply write iS(x—y ) = {fi>(x), (y)}, 
but you should remember that S(x — y ) is a 4 x 4 matrix. Inserting the expansion (5.25), 
we have 


iS(x-y) = 


d 3 p d 3 q 


(2vr) 6 v /4£p; 
d 3 p 1 




- i(p-x-q-y ) 


,+i(p-x-q-y) 


(2vr) 3 2 Ep 
d 3 p 1 
(27t) 3 2 Ep 


+ { C p, Cr q} vS (P) Vr (^ e 

[u s (p)u s (p)e~ ip < x - y) + v s (p)v s (p)e +ip < x - y) ] 

+ m)e- ip ' {x - y) + [fi- m)e +ip - (x ~ y) ] 


(5.27) 
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where to reach the final line we have used the outer product formulae (4.128) and 
(4.129). We can then write 

iS(x-y) = 

in terms of the propagator for a 
as (2.90) 

D(x — 

Some comments: 

• For spacelike separated points (x — y) 2 < 0, we have already seen that D(x — y ) — 
D(y — x) = 0. In the bosonic theory, we made a big deal of this since it ensured 
that 

[<j>(x),<j>{y)\= 0 (x — y) 2 < 0 (5.30) 

outside the lightcone, which we trumpeted as proof that our theory was causal. 
However, for fermions we now have 

{*/> < x{x),il>p(y)} = 0 (x~y) 2 < 0 (5.31) 

outside the lightcone. What happened to our precious causality? The best that 
we can say is that all our observables are bilinear in fermions, for example the 
Hamiltonian (5.17). These still commute outside the lightcone. The theory re¬ 
mains causal as long as fermionic operators are not observable. If you think this is 
a little weak, remember that no one has ever seen a physical measuring apparatus 
come back to minus itself when you rotate by 360 degrees! 

• At least away from singularities, the propagator satisfies 

(i0 x — m)S(x — y) = 0 (5.32) 

which follows from the fact that (0 2 x + m 2 )D(x — y) — 0 using the mass shell 
condition p 2 = m 2 . 

5.5 The Feynman Propagator 

By a similar calculation to that above, we can determine the vacuum expectation value, 
(0| V’oMV’/jfo) |0) = j % ^ E J J + TO )<# e- p( — v) 

(o| Mv)Mx) |0) = / J^fW^~ mUe+ ‘ P<x ~’ ) (5 ' 33) 


(i0 x + m)(D(x-y) - D(y - x)) (5.28) 

real scalar held D(x — y) which, recall, can be written 

») = / < 5 - 29 > 
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We now define the Feynman propagator S F (x — y), which is again a 4 x 4 matrix, as 
the time ordered product, 


S F (x -y) = (0| Tip(x)^(y) |0) 


( 0 | ip(x)ip(y) | 0 ) x° > y° 
( 0 | — 0(?/)0(;r) | 0 ) y° > x° 


Notice the minus sign! It is necessary for Lorentz invariance. When (x—y) 2 < 0, there is 
no invariant way to determine whether x° > y° or y° > x°. In this case the minus sign is 
necessary to make the two definitions agree since {i/>(x),i/j(y)} = 0 outside the lightcone. 
We have the 4-momentum integral representation for the Feynman propagator, 


S F (x - y) 


dA P _ip.(x- v ) 7 ' P + m 

(27r) 4 p 2 — m 2 + ie 


(5.35) 


which satishes {i$ x — m)S F (x — y) = i8^(x — y), so that S F is a Green’s function for 
the Dirac operator. 


The minus sign that we see in (5.34) also occurs for any string of operators inside 
a time ordered product T(...). While bosonic operators commute inside T, fermionic 
operators anti-commute. We have this same behaviour for normal ordered products as 
well, with fermionic operators obeying : 0102 : 0 2 0i With the understanding 

that all fermionic operators anti-commute inside T and Wick’s theorem proceeds 
just as in the bosonic case. We define the contraction 

V>( x)i>{y) = T(^(x)^(y)) - : ^{x)^{y) : = S F (x - y) (5.36) 

5.6 Yukawa Theory 

The interaction between a Dirac fermion of mass m and a real scalar held of mass y is 
governed by the Yukawa theory, 

C = — \y 2 (f) 2 + — m)-0 — (5.37) 

which is the proper version of the baby scalar Yukawa theory we looked at in Section 3. 
Couplings of this type appear in the standard model, between fermions and the Higgs 
boson. In that context, the fermions can be leptons (such as the electron) or quarks. 


Yukawa originally proposed an interaction of this type as an effective theory of nuclear 
forces. With an eye to this, we will again refer to the 0 particles as mesons, and the 
-0 particles as nucleons. Except, this time, the nucleons have spin. (This is still not 
a particularly realistic theory of nucleon interactions, not least because we’re omitting 
isospin. Moreover, in Nature the relevant mesons are pions which are pseudoscalars, so 
a coupling of the form 00y 5 0 would be more appropriate. We’ll turn to this briefly in 
Section 5.7.3). 
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Note the dimensions of the various fields. We still have [0] = 1, but the kinetic 
terms require that [0] = 3/2. Thus, unlike in the case with only scalars, the coupling 
is dimensionless: [A] = 0. 

We’ll proceed as we did in Section 3, firstly computing the amplitude of a particular 
scattering process then, with that calculation as a guide, writing down the Feynman 
rules for the theory. We start with: 

5.6.1 An Example: Putting Spin on Nucleon Scattering 

Let’s study 00 ~^ 00 scattering. This is the same calculation we performed in Section 
(3.3.3) except now the fermions have spin. Our initial and final states are 

\i) = yjAEpEfbp* bf 1 | 0 ) = | p,s\q,r) 

|/) = ^4 Ep'Efbp b\p |0) = \p', s'; q’, r') (5.38) 

We need to be a little cautious about minus signs, because the 0’ s now anti-commute. 
In particular, we should be careful when we take the adjoint. We have 

{/I = \/Spf <0| b^.bi, (5.39) 

We want to calculate the order A 2 terms from the S-matrix element (f\S — 1 |i). 

- y~ J d 4 x 1 cl 4 x 2 T (rjj{xi)ip(x 1 )<j>(x 1 )^(x 2 )^(x 2 )</)(x 2 )) (5.40) 

where, as usual, all helds are in the interaction picture. Just as in the bosonic calcula¬ 
tion, the contribution to nucleon scattering comes from the contraction 

: 0(zi)0(zi)0(z 2 )0(a; 2 ) : (p(x 1 )4>(x 2 ) (5.41) 

We just have to be careful about how the spinor indices are contracted. Let’s start 
by looking at how the fermionic operators act on |i). We expand out the 0 helds, 
leaving the 0 helds alone for now. We may ignore the <0 pieces in 0 since they give no 
contribution at order A 2 . We have 


: 0(xi)0(xi)0(x 2 )0(:r 2 ) : b s J b r J | 0 ) 


d kl [0(ah) ' u m (h)} [0(x 2 ) • u n (k 2 )} 


(2vr) £ 


0 —iki-x\—ik2-X2 




b l b l b p b i 1 °) < 5 ’ 42 ) 
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where we’ve used square brackets [•] to show how the spinor indices are contracted. The 
minus sign that sits out front came from moving Y( x i) past ^(x 2 ). Now anti-commuting 
the 6 ’s past the tf’s, we get 


0 fjT-rr ■u r (q)][^(x 2 )-u s {p}}e lp ' X2 lq ' Xl 

A yj hjphjq 

- ${ Xl ) ■ u s (p}] $(x 2 ) ■ u r (q)] e - ip - Xl - iqX2 ) |0) (5.43) 


Note, in particular, the relative minus sign that appears between these two terms. Now 
let’s see what happens when we hit this with (/1. We look at 


p +ip'-xi+iq'-X2 

<01 b\? hi [Hi) • u r (q)\ [ip(x 2 ) ■ U s (p)} | 0 ) = — -7== [u s \v ) ' u r (q)] [u r '(q) ■ u s (p)} 


2 

p+ip'-X 2 +iq' 

2 \J EfiEqi 


u r '(q')-u r (q)} [u s \p')-u s (p)] 


— S 1 / „.S! 


The [4>(xi) ■ u s (p)\ [i)(x 2 ) • u r (q)] term in (5.43) doubles up with this, cancelling the 
factor of 1/2 in front of (5.40). Meanwhile, the 1 /\[E terms cancel the relativistic 
state normalization. Putting everything together, we have the following expression for 

</|S-l|i> 


HA ) 2 


d 4 xid 4 x 2 d‘ l k ie lk '( Xl X2 ' ) 
(27r) 4 k 2 — p 2 + ie 


([ii s V') • u s (p )] [u r \q') ■ u r (q)\e +ixi < q '- q)+iX2 < p '- p) 

- [u s \p') ■ u r (q)] [u r \q’) ■ u s {p)]e ixv{p '- q)+ixoAq, - p) ^ 


where we’ve put the </> propagator back in. Performing the integrals over x\ and x 2 , 
this becomes, 

[ d A k ^^ ^ ^ ([ u s \p’) ■ u s (p)\ [u r \q ') ■ u r (q)\5^\q' -q + k)5 (4 \p' - p - k) 

J k z — p- + ie V 

— [u s \p') ■ u r (q)\ [u r '(q') ■ u s (p)\8^\p' — q + k)5 <yi \q' — p — k)^j 

And we’re almost there! Finally, writing the S-matrix element in terms of the amplitude 
in the usual way, (f \ S — 1 |z) = iA{2'K) 4 5 (y4 Yp + q — p' — </), we have 

= , 2 f [u s \p') • u s (p )] [u r \q') • u r (q)\ _ [u s \p') • u r (q)\ [.u r \q') • u s (p)\ 

\ {p' — p) 2 — p 2 -Vie ( q' — p) 2 — p 2 -Vie 

which is our final answer for the amplitude. 
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5.7 Feynman Rules for Fermions 

It’s important to bear in mind that the calculation we just did kind of blows. Thankfully 
the Feynman rules will once again encapsulate the combinatoric complexities and make 
life easier for us. The rules to compute amplitudes are the following 

• To each incoming fermion with momentum p and spin r, we associate a spinor 
u r (p). For outgoing fermions we associate u r (p). 


u r (p) 




Figure 21: An incoming fermion 




u r Cp) 


Figure 22: An outgoing fermion 


• To each incoming anti-fermion with momentum p and spin r, we associate a spinor 
v r (p). For outgoing anti-fermions we associate v r {p). 



Figure 23: An incoming anti-fermion Figure 24: An outgoing anti-fermion 


• Each vertex gets a factor of —iX. 


• Each internal line gets a factor of the relevant propagator. 


p 


p 


> 


-> 


i 

p 2 — /i 2 + ie 
i(]j + m) 

p2 _ m 2 _|_ 


for scalars 
for fermions 


(5.44) 


The arrows on the fermion lines must flow consistently through the diagram (this 
ensures fermion number conservation). Note that the fermionic propagator is a 
4x4 matrix. The matrix indices are contracted at each vertex, either with further 
propagators, or with external spinors u, u, v or v. 


• Impose momentum conservation at each vertex, and integrate over undetermined 
loop momenta. 

• Add extra minus signs for statistics. Some examples will be given below. 
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5.7.1 Examples 

Let’s run through the same examples we did for the scalar Yukawa theory. Firstly, we 
have 

Nucleon Scattering 

For the example we worked out previously, the two lowest order Feynman diagrams are 
shown in Figure 25. We’ve drawn the second Feynman diagram with the legs crossed 




Figure 25: The two Feynman diagrams for nucleon scattering 


to emphasize the fact that it picks up a minus sign due to statistics. (Note that the 
way the legs point in the Feynman diagram doesn’t tell us the direction in which the 
particles leave the scattering event: the momentum label does that. The two diagrams 
above are different because the incoming legs are attached to different outgoing legs). 
Using the Feynman rules we can read off the amplitude. 

. 2 / [ u s \p') • u s (p )] [u r \g') • u r (g)} [u s \p') • u r (q)} [u r ' (<f) • u s (p )] 

1 ' V {p-p'Y-p 2 {p-q'f-V 2 



The denominators in each term are due to the meson propagator, with the momen¬ 
tum determined by conservation at each vertex. This agrees with the amplitude we 
computed earlier using Wick’s theorem. 


Nucleon to Meson Scattering 

Let’s now look at 00 — > 00. The two lowest order Feynman diagrams are shown in 
Figure 26. Applying the Feynman rules, we have 

_ . 2 ( v\q) [7~ P'p) + m}u s (p) v r {g)[Y(p» ~ gjj + rn]u s (p) \ 

\ (p — p ') 2 — m 2 (p — q ') 2 — m 2 ) 

Since the internal line is now a fermion, the propagator contains 7 n(P/j,—p' IJ )+'i r n factors. 
This is a 4 x 4 matrix which sits on the top, sandwiched between the two external 
spinors. Now the exchange statistics applies to the final meson states. These are 
bosons and, correspondingly, there is no relative minus sign between the two diagrams. 
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Figure 26: The two Feynman diagrams for nucleon to meson scattering 


Nucleon-Anti-Nucleon Scattering 

For —> 'i/j'ij), the two lowest order Feynman diagrams are of two distinct types, just 

like in the bosonic case. They are shown in Figure 27. 

The corresponding amplitude is given by, 


A 


( [u a '(p') -u s {p)} [v r {g) ■v r \q’)} [v r (q) ■ u s (p )] [u B> jf) • v r '{q')} \ , , 

V (p-p') 2 -h 2 (p + q) 2 -li 2 + ie ) [ - ’ 


As in the bosonic diagrams, there is again the difference in the momentum dependence 
in the denominator. But now the difference in the diagrams is also reflected in the 
spinor contractions in the numerator. 


More subtle are the minus signs. The fermionic statistics mean that the first diagram 
has an extra minus sign relative to the scattering of Figure 25. Since this minus sign 
will be important when we come to figure out whether the Yukawa force is attractive 
or repulsive, let’s go back to basics and see where it comes from. The initial and final 
states for this scattering process are 

I*) = V 4E pEqbp c r A |0) = I p,s-,q,r) 

|/) = yj-bp cp |0) = |p", s '; q\r’) (5.47) 

The ordering of and A in these states is crucial and reflects the scattering r/j-i/j —y ipi/j, 
as opposed to , 0'0 '0'0 which would differ by a minus sign. The first diagram in 

Figure 27 comes from the term in the perturbative expansion, 

(/I : A(xi)A{xi)A(x 2 )A(x 2 ) : b s J c r J |0) ~ (/1 [b m (£i) • ^(xi)] [^(x 2 ) • u n c r J |0) 

where we’ve neglected a bunch of objects in this equation like f d A ki and exponential 
factors because we only want to keep track of the minus signs. Moving the annihilation 
operators past the creation operators, we have 

+ (/I [v r (q) ■ A(xi)} $(x 2 ) • u s (p )] |0) (5.48) 
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Figure 27: The two Feynman diagrams for nucleon-anti-nucleon scattering 


Repeating the process by expanding out the i/j(x i) and ^(^ 2 ) fields and moving them 
to the left to annihilate (/|, we have 


4bi,c™' 6” f [v r (q] ■ v m (h)\ [S”(f 2 ) . u‘(p ,|] |0) ~ ~{v r (q) ■ V'(<f)] 


u 


where the minus sign has appeared from anti-commuting cH past b s L. This is the 
overall minus sign found in (5.46). One can also follow similar contractions to compute 
the second diagram in Figure 27. 


Meson Scattering 


Finally, we can also compute the scattering of 00 —> 00 
which, as in the bosonic case, picks up its leading contribution 
at one-loop. The amplitude for the diagram shown in the figure 

is 



lA 


-HA ) 4 


d A k tf + m tf+ A + m 

(27t) 4 (k 2 — m 2 + ie) {{k + p'0 2 — m 2 + ie) Figure 28: 

A + m - A + m 

((k + p[ — Pi ) 2 — rn 2 + ie) (( k — p '^) 2 — m 2 + ie) 


Notice that the high momentum limit of the integral is f d^k/k 4 , which is no longer 
finite, but diverges logarithmically. You will have to wait until next term to make sense 
of this integral. 


There’s an overall minus sign sitting in front of this amplitude. This is a generic 
feature of diagrams with fermions running in loops: each fermionic loop in a diagram 
gives rise to an extra minus sign. We can see this rather simply in the diagram 
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which involves the expression 







fpa(x)ip a (x)ipy(y)^y(y) 


~^y{y)^a{x) ip a (x)ipp(y) 
-Tr (S F (y-x)S F (x-y)) 


After passing the fermion fields through each other, a minus sign appears, sitting in 
front of the two propagators. 


5.7.2 The Yukawa Potential Revisited 

We saw in Section 3.5.2, that the exchange of a real scalar particle gives rise to a 
universally attractive Yukawa potential between two spin zero particles. Does the 
same hold for the spin 1/2 particles? 

Recall that the strategy to compute the potential is to take the non-relativistic limit 
of the scattering amplitude, and compare with the analogous result from quantum 
mechanics. Our new amplitude now also includes the spinor degrees of freedom u(p) 
and v(p). In the non-relativistic limit, p —> ( m,p ), and 


u(p) = 


v(p) = 


V Vp 7 ^^ J 

- Vp ^ J 



(5.49) 


In this limit, the spinor contractions in the amplitude for 00 —> scattering (5.45) 

become u s ' ■ u s = 2 m8 ss ' and the amplitude is 



— i{—i\) 2 (2m) 


8 ss 8 r 


8 sr 8 r 


(. p-p r )+y 2 (p-^+p 2 


(5.50) 


The 8 symbols tell us that spin is conserved in the non-relativistic limit, while the 
momentum dependence is the same as in the bosonic case, telling us that once again 
the particles feel an attractive Yukawa potential, 

A 2 e -Mr 

U(r) = (5.51) 


Repeating the calculation for 'ijj'ijj —> -0^, there are two minus signs which cancel each 
other. The first is the extra overall minus sign in the scattering amplitude (5.46), 
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due to the fermionic nature of the particles. The second minus sign comes from the 
non-relativistic limit of the spinor contraction for anti-particles in (5.46), which is 
v s ' ■ v s = —2 mS ss '. These two signs cancel, giving us once again an attractive Yukawa 
potential (5.51). 

5.7.3 Pseudo-Scalar Coupling 

Rather than the standard Yukawa coupling, we could instead consider 

£y u k = -A0^7 5 ^ (5.52) 

This still preserves parity if 0 is a pseudoscalar, i.e. 

P : <j)(x, t) —> —(f>(—x,t) (5.53) 


We can compute in this theory very simply: the Feynman rule for the interaction vertex 
is now changed to a factor of — iXj 5 . For example, the Feynman diagrams for ip'ip —» 'ijj'i/j 
scattering are again given by Figure 25, with the amplitude now 

2 f[u s ’(f)'fu a (p)] [u r '{q')^u r {q)\ [u s '{p')^u r {q)\ [u r \q')^u s {p)} 

[ 1 ’ V (p-p') 2 -p 2 (p-q'Y-p 2 

We could again try to take the non-relativistic limit for this amplitude. But this 
time, things work a little differently. Using the expressions for the spinors (5.49), we 
have u s '^f 5 u s —> 0 in the non-relativistic limit. To find the non-relativistic amplitude, 
we must go to next to leading order. One can easily check that u s ' (p r ) r 'f 5 u s (p) —» 
m£ s ' T (p — p') •cr£ s . So, in the non-relativistic limit, the leading order amplitude arising 
from pseudoscalar exchange is given by a spin-spin coupling, 



[z s ' T (p- f) • [c' T (p- f) • 

ip-p ') 2 + p 2 


(5.54) 
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6. Quantum Electrodynamics 


In this section we finally get to quantum electrodynamics (QED), the theory of light 
interacting with charged matter. Our path to quantization will be as before: we start 
with the free theory of the electromagnetic field and see how the quantum theory gives 
rise to a photon with two polarization states. We then describe how to couple the 
photon to fermions and to bosons. 


6.1 Maxwell’s Equations 

The Lagrangian for Maxwell’s equations in the absence of any sources is simply 

£ = — | F iiv F' w (6.1) 

where the field strength is defined by 


F„v = S r K - d„A„ 


( 6 . 2 ) 


The equations of motion which follow from this Lagrangian are 




= 0 


(6.3) 


Meanwhile, from the definition of the field strength also satisfies the Bianchi 
identity 


dx Ffiv + d/iFux + dyFx,,, — 0 


(6.4) 


To make contact with the form of Maxwell’s equations you learn about in high school, 
we need some 3-vector notation. If we define A M = (</>, A), then the electric field E and 
magnetic field B are defined by 


_ Qj\ 

E = —V(p —— and B = V x A 
at 


which, in terms of F ;i „, becomes 


0 Ey E z 

-E x 0 —B z By 


F = 

A**' I -Ey B z 0 -B x 


— E z —By B x 0 


The Bianchi identity (6.4) then gives two of Maxwell’s equations, 


d B 

V • B = 0 and —— = —V x F 

dt 


(6.5) 


( 6 . 6 ) 


(6.7) 
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These remain true even in the presence of electric sources. Meanwhile, the equations 
of motion give the remaining two Maxwell equations, 

_ d E ^ 

V • E = 0 and = V x B (6.8) 

dt y J 

As we will see shortly, in the presence of charged matter these equations pick up extra 
terms on the right-hand side. 


6.1.1 Gauge Symmetry 

The massless vector field A M has 4 components, which would naively seem to tell us that 
the gauge field has 4 degrees of freedom. Yet we know that the photon has only two 
degrees of freedom which we call its polarization states. How are we going to resolve 
this discrepancy? There are two related comments which will ensure that quantizing 
the gauge field A M gives rise to 2 degrees of freedom, rather than 4. 

• The field Aq has no kinetic term Aq in the Lagrangian: it is not dynamical. This 
means that if we are given some initial data A{ and Ai at a time to, then the field 
Aq is fully determined by the equation of motion V ■ E = 0 which, expanding out, 
reads 

BA 

V 2 A 0 + V ■ — = 0 (6.9) 

This has the solution 


A 0 (x) 


3 (V -dA/dt)(x') 
47t|x — x'\ 


( 6 . 10 ) 


So A 0 is not independent: we don’t get to specify A 0 on the initial time slice. It 
looks like we have only 3 degrees of freedom in A^ rather than 4. But this is still 
one too many. 


• The Lagrangian (6.3) has a very large symmetry group, acting on the vector 
potential as 


A^x) ->• A^x) + d^X(x) (6.11) 

for any function \(x). We’ll ask only that X(x) dies off suitably quickly at spatial 
x —» oo. We call this a gauge symmetry. The field strength is invariant under the 
gauge symmetry: 


Ffj, u d^{A v + d u X) — d u (A^ + d tl X) — F fW (6.12) 
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So what are we to make of this? We have a theory with an infinite number of 
symmetries, one for each function A (a;). Previously we only encountered symme¬ 
tries which act the same at all points in spacetime, for example ^ —> e ia ty for a 
complex scalar field. Noether’s theorem told us that these symmetries give rise 
to conservation laws. Do we now have an infinite number of conservation laws? 

The answer is no! Gauge symmetries have a very different interpretation than 
the global symmetries that we make use of in Noether’s theorem. While the 
latter take a physical state to another physical state with the same properties, 
the gauge symmetry is to be viewed as a redundancy in our description. That is, 
two states related by a gauge symmetry are to be identified: they are the same 
physical state. (There is a small caveat to this statement which is explained in 
Section 6.3.1). One way to see that this interpretation is necessary is to notice 
that Maxwell’s equations are not sufficient to specify the evolution of A The 
equations read, 


[n,w(d P d P ) - <9,A] A v = 0 (6.13) 

But the operator [r]^i / {d p dp) — d^d l ] is not invertible: it annihilates any function of 
the form d^X. This means that given any initial data, we have no way to uniquely 
determine A M at a later time since we can’t distinguish between A M and A /t + <9 /t A. 
This would be problematic if we thought that A M is a physical object. However, 
if we’re happy to identify A M and A M + d^X as corresponding to the same physical 
state, then our problems disappear. 

Since gauge invariance is a redundancy of the system, 
we might try to formulate the theory purely in terms of 
the local, physical, gauge invariant objects E and B. This 
is fine for the free classical theory: Maxwell’s equations 
were, after all, first written in terms of E and B. But it is 
not possible to describe certain quantum phenomena, such 
as the Aharonov-Bohm effect, without using the gauge 
potential A M . We will see shortly that we also require the Figure 29- 

gauge potential to describe classically charged fields. To 

describe Nature, it appears that we have to introduce quantities A M that we can never 
measure. 

The picture that emerges for the theory of electromagnetism is of an enlarged phase 
space, foliated by gauge orbits as shown in the figure. All states that he along a given 
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line can be reached by a gauge transformation and are identified. To make progress, 
we pick a representative from each gauge orbit. It doesn’t matter which representative 
we pick — after all, they’re all physically equivalent. But we should make sure that we 
pick a “good” gauge, in which we cut the orbits. 

Different representative configurations of a physical state are called different gauges. 
There are many possibilities, some of which will be more useful in different situations. 
Picking a gauge is rather like picking coordinates that are adapted to a particular 
problem. Moreover, different gauges often reveal slightly different aspects of a problem. 
Here we’ll look at two different gauges: 

• Lorentz Gauge: = 0 

To see that we can always pick a representative configuration satisfying d^A^ = 0, 
suppose that we’re handed a gauge held A' satisfying d^A'Y = fix). Then we 
choose An = A' + d^X, where 


d,d' l \ = ~f (6.14) 

This equation always has a solution. In fact this condition doesn’t pick a unique 
representative from the gauge orbit. We’re always free to make further gauge 
transformations with <9 /t <9 /J A = 0, which also has non-trivial solutions. As the 
name suggests, the Lorentz gauge 3 has the advantage that it is Lorentz invariant. 

• Coulomb Gauge: V ■ A = 0 

We can make use of the residual gauge transformations in Lorentz gauge to pick 
V ■ A = 0. (The argument is the same as before). Since A 0 is fixed by (6.10), we 
have as a consequence 


A 0 = 0 (6.15) 

(This equation will no longer hold in Coulomb gauge in the presence of charged 
matter). Coulomb gauge breaks Lorentz invariance, so may not be ideal for some 
purposes. However, it is very useful to exhibit the physical degrees of freedom: 
the 3 components of A satisfy a single constraint: V • A = 0, leaving behind just 
2 degrees of freedom. These will be identified with the two polarization states of 
the photon. Coulomb gauge is sometimes called radiation gauge. 

3 Named after Lorenz who had the misfortune to be one letter away from greatness. 
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6.2 The Quantization of the Electromagnetic Field 

In the following we shall quantize free Maxwell theory twice: once in Coulomb gauge, 
and again in Lorentz gauge. We’ll ultimately get the same answers and, along the way, 
see that each method comes with its own subtleties. 


The first of these subtleties is common to both methods and comes when computing 
the momentum 7 conjugate to A^, 


n 


7r 


8C 

d A 0 
DC 

dli 


0 

_ F 0i = E l 


(6.16) 


so the momentum 7T° conjugate to To vanishes. This is the mathematical consequence of 
the statement we made above: A 0 is not a dynamical held. Meanwhile, the momentum 
conjugate to Tj is our old friend, the electric held. We can compute the Hamiltonian, 


H — I d 3 x 7T i Ai - C 

= j d 3 x ^E ■ E + \B ■ B — T 0 (V ■ E) (6.17) 

So A 0 acts as a Lagrange multiplier which imposes Gauss’ law 

V-E = 0 (6.18) 

which is now a constraint on the system in which A are the physical degrees of freedom. 
Let’s now see how to treat this system using different gauge hxing conditions. 


6.2.1 Coulomb Gauge 

In Coulomb gauge, the equation of motion for A is 

8^1=0 (6.19) 


which we can solve in the usual way, 



with pi = |p) 2 . The constraint V • A = 0 tells us that £ must satisfy 

Z-p = o 


( 6 . 20 ) 


( 6 . 21 ) 
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which means that £ is perpendicular to the direction of motion p. We can pick £ (p) to 
be a linear combination of two orthonormal vectors e r , r = 1,2, each of which satisfies 
e r (p) ■ p — 0 and 

trip) -tsip) = &rs r,s = 1,2 (6.22) 


These two vectors correspond to the two polarization states of the photon. It’s worth 
pointing out that you can’t consistently pick a continuous basis of polarization vectors 
for every value of p because you can’t comb the hair on a sphere. But this topological 
fact doesn’t cause any complications in computing QED scattering processes. 

To quantize we turn the Poisson brackets into commutators. Naively we would write 

[A i (x),A j (y)\ = lE i (x),E>(y)} = 0 
[4(f), EP(y)\ = iSj h (3) (f - y) (6.23) 


But this can’t quite be right, because it’s not consistent with the constraints. We 
still want to have V • A = V ■ E = 0, now imposed on the operators. But from the 
commutator relations above, we see 

[V ■ A(x),V ■ E(y)] =iV 2 5 {3 \x-y)^0 (6.24) 


What’s going on? In imposing the commutator relations (6.23) we haven’t correctly 
taken into account the constraints. In fact, this is a problem already in the classical 
theory, where the Poisson bracket structure is already altered 4 . The correct Poisson 
bracket structure leads to an alteration of the last commutation relation, 


[Ai(x),Ej(y)\ = i (da - ) h (3) (f - y) 


(6.25) 


To see that this is now consistent with the constraints, we can rewrite the right-hand 
side of the commutator in momentum space, 



(6.26) 


which is now consistent with the constraints, for example 


[diAi(x), Ej(y)\ 



d 3 p 

(2vr) 3 


Sij 


PiPj 

\P\ 2 


i Pi e= 0 


(6.27) 


4 For a nice discussion of the classical and quantum dynamics of constrained systems, see the small 
book by Paul Dirac, “Lectures on Quantum Mechanics ” 
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We now write A in the usual mode expansion, 


A(x) = 

a- 

CO 

2 

J2?r(p) + 

r —1 

(2vr ) 3 v ^f 

E{x) = ^ 


y X>W> [&** - a'/e-** 

r=1 


where, as before, the polarization vectors satisfy 


(6.28) 


e r (p) ■ P = 0 and e r (p) ■ e s (p) = S rs 


(6.29) 


It is not hard to show that the commutation relations (6.25) are equivalent to the usual 
commutation relations for the creation and annihilation operators, 

a ?-] = ] = 0 

[ah, aJ] = {2ir) 3 S rs 6®(p- q) (6.30) 

where, in deriving this, we need the completeness relation for the polarization vectors, 

E4(p>;'(p) = ^-^ (6.3i) 

You can easily check that this equation is true by acting on both sides with a basis of 
vectors (ei(p), e 2 (p),p). 

We derive the Hamiltonian by substituting (6.28) into (6.17). The last term vanishes 
in Coulomb gauge. After normal ordering, and playing around with e r polarization 
vectors, we get the simple expression 



The Coulomb gauge has the advantage that the physical degrees of freedom are man¬ 
ifest. However, we’ve lost all semblance of Lorentz invariance. One place where this 
manifests itself is in the propagator for the fields A^x) (in the Heisenberg picture). In 
Coulomb gauge the propagator reads 

D%(x - y) = <0| TMx)My) 0) = / (s„ - e~^ (6.33) 

The tr superscript on the propagator refers to the “transverse” part of the photon. 
When we turn to the interacting theory, we will have to fight to massage this propagator 
into something a little nicer. 
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6.2.2 Lorentz Gauge 

We could try to work in a Lorentz invariant fashion by imposing the Lorentz gauge 
condition d^A' 1 = 0. The equations of motion that follow from the action are then 

d il &‘A v = 0 (6.34) 

Our approach to implementing Lorentz gauge will be a little different from the method 
we used in Coulomb gauge. We choose to change the theory so that (6.34) arises directly 
through the equations of motion. We can achieve this by taking the Lagrangian 

C = - 1 -F, V F^ (6.35) 

The equations of motion coming from this action are 


d l ,F IJ,v + d v {d^) = d lt & J A v = 0 


(6.36) 


(In fact, we could be a little more general than this, and consider the Lagrangian 

C = -\F, U F^ - ^(<V^) 2 (6.37) 

with arbitrary a and reach similar conclusions. The quantization of the theory is 
independent of a and, rather confusingly, different choices of a are sometimes also 
referred to as different “gauges”. We will use a = 1, which is called “Feynman gauge”. 
The other common choice, a = 0, is called “Landau gauge”.) 


Our plan will be to quantize the theory (6.36), and only later impose the constraint 
= 0 in a suitable manner on the Hilbert space of the theory. As we’ll see, we will 
also have to deal with the residual gauge symmetry of this theory which will prove a 
little tricky. At first, we can proceed very easily, because both 7T° and 7i l are dynamical: 


7r 


7r 


dC 

DAq 

dC 

dAi 


-d,A> J ' 
<TA° - A 1 


(6.38) 


Turning these classical fields into operators, we can simply impose the usual commu¬ 
tation relations, 


[A^x),A v {y)\ = [t ^{x),Tr v {y)\ = 0 
[Ap{x), tg(2/)] = - y) (6.39) 
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and we can make the usual expansion in terms of creation and annihilation operators 
and 4 polarization vectors (e M ) A , with A = 0,1, 2,3. 


A n(v) = 


■kHx) = 


d 3 p 


(2ir) 3 p2\p\ ^ 


J24(p> ' i12 


% M (+<) 


(2tt 


A=0 


a±e ip - s -afe-^"- 


(6.40) 


Note that the momentum 7 T m comes with a factor of (+*), rather than the familiar (— i) 
that we’ve seen so far. This can be traced to the fact that the momentum (6.38) for the 
classical fields takes the form 7 t m = — + .... In the Heisenberg picture, it becomes 

clear that this descends to (+i) in the definition of momentum. 

There are now four polarization 4-vectors e A (p), instead of the two polarization 3- 
vectors that we met in the Coulomb gauge. Of these four 4-vectors, we pick e° to be 
timelike, while e 1,2,3 are spacelike. We pick the normalization 


e 


A 


■ e 


A' 



(6.41) 


which also means that 


(cO A (e„) x ' r) aa' = rj ^ (6.42) 

The polarization vectors depend on the photon 4-momentum p = (\p\,p)- Of the two 
spacclike polarizations, we will choose e 1 and e 2 to lie transverse to the momentum: 

e 1 • p = e 2 ■ p = 0 (6.43) 

The third vector e 3 is the longitudinal polarization. For example, if the momentum lies 
along the x 3 direction, so p ~ (1, 0, 0,1), then 



For other 4-momenta, the polarization vectors are the appropriate Lorentz transforma¬ 
tions of these vectors, since (6.43) are Lorentz invariant. 

We do our usual trick, and translate the field commutation relations (6.39) into those 
for creation and annihilation operators. We find [ap,a~] = [aV,a^ *] = 0 and 

= -V XX ’ (2v r) 3 6®(p- q) (6.45) 
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The minus signs here are odd to say the least! For spacclikc A = 1,2,3, everything 
looks fine, 

[4, < f ] = <5 AV (2vr) 3 5&(p - q) A, A' = 1, 2, 3 (6.46) 

But for the timelike annihilation and creation operators, we have 

[4’ a ?] = -( 27 0 3 <5 (3 ) (p- q) (6-47) 

This is very odd! To see just how strange this is, we take the Lorentz invariant vacuum 
|0) defined by 


4 |0> = 0 (6.48) 

Then we can create one-particle states in the usual way, 

|p,A) = a At |0) (6.49) 

For spacelike polarization states, A = 1,2,3, all seems well. But for the timelike 
polarization A = 0, the state | p, 0) has negative norm, 

(p, 0| q, 0) = (0| 4 aV |0) = -(2tt) 3 5^\p - q) (6.50) 

Wtf? That’s very very strange. A Hilbert space with negative norm means negative 
probabilities which makes no sense at all. We can trace this negative norm back to the 
wrong sign of the kinetic term for A 0 in our original Lagrangian: C = + \A 2 — + .... 

At this point we should remember our constraint equation, d^A^ = 0, which, until 
now, we’ve not imposed on our theory. This is going to come to our rescue. We will see 
that it will remove the timelike, negative norm states, and cut the physical polarizations 
down to two. We work in the Heisenberg picture, so that 

d^A 1 - 1 = 0 (6.51) 

makes sense as an operator equation. Then we could try implementing the constraint 
in the quantum theory in a number of different ways. Let’s look at a number of 
increasingly weak ways to do this 

• We could ask that = 0 is imposed as an equation on operators. But this 

can’t possibly work because the commutation relations (6.39) won’t be obeyed 
for 7r° = —d^A^. We need some weaker condition. 
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• We could try to impose the condition on the Hilbert space instead of directly 
on the operators. After all, that’s where the trouble lies! We could imagine that 
there’s some way to split the Hilbert space up into good states |\&) and bad states 
that somehow decouple from the system. With luck, our bad states will include 
the weird negative norm states that we’re so disgusted by. But how can we define 
the good states? One idea is to impose 


|t> = o 


(6.52) 


on all good, physical states |\P). But this can’t work either! Again, the condition 
is too strong. For example, suppose we decompose A^ix) = A+(x) + A~(x) with 


Ai(x) = 

A ~(x) = 


d 3 p 


— y 

(2ir) 3 y 2 p 


44 ^ 


d 3 p 


— y 

pir ) 3 


44 ' e+ip '‘ 


(6.53) 


Then, on the vacuum A+ |0) = 0 automatically, but d^A^ |0) ^ 0. So not even 
the vacuum is a physical state if we use (6.52) as our constraint 


• Our final attempt will be the correct one. In order to keep the vacuum as a good 
physical state, we can ask that physical states T) are defined by 


<9A4+|T) = 0 


(6.54) 


This ensures that 


= 0 (6.55) 

so that the operator has vanishing matrix elements between physical states. 
Equation (6.54) is known as the Gupta-Bleuler condition. The linearity of the 
constraint means that the physical states T) span a physical Hilbert space "Kphys- 

So what does the physical Hilbert space "Kphys look like? And, in particular, have we 
rid ourselves of those nasty negative norm states so that "Hphys has a positive definite 
inner product defined on it? The answer is actually no, but almost! 

Let’s consider a basis of states for the Fock space. We can decompose any element 
of this basis as |T) = |^ r ) 10), where l^r) contains only transverse photons, created by 
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ah 2 \ while |0) contains the timelike photons created by aV and longitudinal photons 
created by a|d. The Gupta-Blculer condition (6.54) requires 

(4-4)|0> = O (6.56) 

This means that the physical states must contain combinations of timelike and longi¬ 
tudinal photons. Whenever the state contains a timelike photon of momentum p, it 
must also contain a longitudinal photon with the same momentum. In general |0) will 
be a linear combination of states |0 n ) containing n pairs of timelike and longitudinal 
photons, which we can write as 

OO 

|0) =Y J C n\(t>n) (6.57) 

71=0 

where |0o) = |0) is simply the vacuum. It’s not hard to show that although the condition 
(6.56) does indeed decouple the negative norm states, all the remaining states involving 
timelike and longitudinal photons have zero norm 

(0m| 0ri) ^nO^mO (6.58) 


This means that the inner product on "Hphys is positive semi-definite. Which is an 
improvement. But we still need to deal with all these zero norm states. 

The way we cope with the zero norm states is to treat them as gauge equivalent 
to the vacuum. Two states that differ only in their timelike and longitudinal photon 
content, |0 n ) with n > 1 are said to be physically equivalent. We can think of the gauge 
symmetry of the classical theory as descending to the Hilbert space of the quantum 
theory. Of course, we can’t just stipulate that two states are physically identical unless 
they give the same expectation value for all physical observables. We can check that 
this is true for the Hamiltonian, which can be easily computed to be 



But the condition (6.56) ensures that (T| o'pai T) = (T| aVa^ |T) so that the contri¬ 
butions from the timelike and longitudinal photons cancel amongst themselves in the 
Hamiltonian. This also renders the Hamiltonian positive definite, leaving us just with 
the contribution from the transverse photons as we would expect. 

In general, one can show that the expectation values of all gauge invariant operators 
evaluated on physical states are independent of the coefficients C n in (6.57). 
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Propagators 

Finally, it’s a simple matter to compute the propagator in Lorentz gauge. It is given 
by 

( 0 | TA„(x)A„(y) | 0 ) = / ^ (6.60) 

This is a lot nicer than the propagator we found in Coulomb gauge: in particular, it’s 
Lorentz invariant. We could also return to the Lagrangian (6.37). Had we pushed 
through the calculation with arbitrary coefficient a, we would find the propagator, 

(0\TA,(x)A„(y)\0) = J ^ ^ („„„ + ( a - 1 )^) e-V-*) (6.61) 

6.3 Coupling to Matter 

Let’s now build an interacting theory of light and matter. We want to write down 
a Lagrangian which couples A^ to some matter fields, either scalars or spinors. For 
example, we could write something like 

C = -\F IU ,F' U ' (6.62) 

where is some function of the matter fields. The equations of motion read 

d IJ F' 11 ' = f (6.63) 

so, for consistency, we require 

d tl f = 0 (6.64) 

In other words, must be a conserved current. But we’ve got lots of those! Let’s look 
at how we can couple two of them to electromagnetism. 

6.3.1 Coupling to Fermions 

The Dirac Lagrangian 


C = ifi(i $ — m)*p (6.65) 

has an internal symmetry ip — y e~ ia, ip and Pp — y e +lo "ip, with a £ R. This gives rise to 
the conserved current jy = ^ 7 ^ip. So we could look at the theory of electromagnetism 
coupled to fermions, with the Lagrangian, 

C = -\F^ U F^ U + ip{i0 - m)ip - exp^A^ip (6.66) 
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where we’ve introduced a coupling constant e. For the free Maxwell theory, we have 
seen that the existence of a gauge symmetry was crucial in order to cut down the 
physical degrees of freedom to the requisite 2. Does our interacting theory above still 
have a gauge symmetry? The answer is yes. To see this, let’s rewrite the Lagrangian 
as 


C = -\F lu/ F lu ' + ^(i]p-m)iP (6.67) 

where D /t -0 = d^i/j + ieA^ is called the covariant derivative. This Lagrangian is 
invariant under gauge transformations which act as 

Afj, — » + d^X and ^ —» e~ ieX i\) (6.68) 

for an arbitrary function X(x). The tricky term is the derivative acting on %!>, since this 
will also hit the e _ieA piece after the transformation. To see that all is well, let’s look 
at how the covariant derivative transforms. We have 


D = d^ + ieA^ 

-»■ <9 M (e _ie V) + ie{Ap + d ft X)(e~ ieX ^) 

= e~ ieX D^ (6.69) 

so the covariant derivative has the nice property that it merely picks up a phase under 
the gauge transformation, with the derivative of e~ ieX cancelling the transformation 
of the gauge held. This ensures that the whole Lagrangian is invariant, since —y 


Electric Charge 

The coupling e has the interpretation of the electric charge of the ijj particle. This 
follows from the equations of motion of classical electromagnetism d tJ F ,llJ = j v \ we 
know that the j° component is the charge density. We therefore have the total charge 
Q given by 


Q = e 


d 3 x^(x) , y 0 'ip(x) 


After treating this as a quantum equation, we have 


Q = e 




(6.70) 


(6.71) 


which is the number of particles, minus the number of antiparticles. Note that the 
particle and the anti-particle are required by the formalism to have opposite electric 
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charge. For QED, the theory of light interacting with electrons, the electric charge 
is usually written in terms of the dimensionless ratio cc, known as the fine structure 
constant 


a = 



47 rhc 137 


Setting h — c— 1, we have e = \Ja ~ 0.3. 


(6.72) 


There’s a small subtlety here that’s worth elaborating on. I stressed that there’s a 
radical difference between the interpretation of a global symmetry and a gauge symme¬ 
try. The former takes you from one physical state to another with the same properties 
and results in a conserved current through Noether’s theorem. The latter is a redun¬ 
dancy in our description of the system. Yet in electromagnetism, the gauge symmetry 
ij) —» e+* eA ( a: )-0 seems to lead to a conservation law, namely the conservation of electric 
charge. This is because among the infinite number of gauge symmetries parameterized 
by a function A(x), there is also a single global symmetry: that with A(x) = constant. 
This is a true symmetry of the system, meaning that it takes us to another physical 
state. More generally, the subset of global symmetries from among the gauge symme¬ 
tries are those for which A(x) —> a = constant as a: —y oo. These take us from one 
physical state to another. 

Finally, let’s check that the 4x4 matrix C that we introduced in Section 4.5 really 
deserves the name “charge conjugation matrix”. If we take the complex conjugation of 
the Dirac equation, we have 

(«7^ - eq^A* - m)ip = 0 => - e( 7 #1 )M #t - m)^* = 0 

Now using the defining equation = — (q M )*, and the definition t/A c ) = Cip*, we 

see that the charge conjugate spinor satisfies 

(i^d^ + = 0 (6.73) 

So we see that the charge conjugate spinor satisfies the Dirac equation, but with 
charge —e instead of +e. 


6.3.2 Coupling to Scalars 

For a real scalar field, we have no suitable conserved current. This means that we can’t 
couple a real scalar field to a gauge field. 
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Let’s now consider a complex scalar field p. (For this section, I’ll depart from our 
previous notation and call the scalar held p to avoid confusing it with the spinor). We 
have a symmetry p —» e~ ia p. We could try to couple the associated current to the 
gauge held, 

Ant = -i^d^p^p - p*d fi p)A tl (6.74) 

But this doesn’t work because 

• The theory is no longer gauge invariant 

• The current j^ that we coupled to A M depends on d^p. This means that if we 
try to compute the current associated to the symmetry, it will now pick up a 
contribution from the j^A^ term. So the whole procedure wasn’t consistent. 

We solve both of these problems simultaneously by remembering the covariant deriva¬ 
tive. In this scalar theory, the combination 

V^p = d^p + ieA^p (6.75) 

again transforms as T>^p —> e~ ieX T> iL p under a gauge transformation + d^X 

and p —> e~ ieX p. This means that we can construct a gauge invariant action for a 
charged scalar held coupled to a photon simply by promoting all derivatives to covariant 
derivatives 

C = + V^V^p - m 2 \p\ 2 (6.76) 

In general, this trick works for any theory. If we have a U(l) symmetry that we wish to 
couple to a gauge held, we may do so by replacing all derivatives by suitable covariant 
derivatives. This procedure is known as minimal coupling. 

6.4 QED 

Let’s now work out the Feynman rules for the full theory of quantum electrodynamics 
(QED) - the theory of electrons interacting with light. The Lagrangian is 

C = ~ l -F^ v + $(ilp - m)V> (6.77) 

where D M + ieA 

The route we take now depends on the gauge choice. If we worked in Lorentz gauge 
previously, then we can jump straight to Section 6.5 where the Feynman rules for QED 
are written down. If, however, we worked in Coulomb gauge, then we still have a bit of 
work in front of us in order to massage the photon propagator into something Lorentz 
invariant. We will now do that. 
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In Coulomb gauge V ■ A = 0, the equation of motion arising from varying A 0 is now 


-V 2 A 0 = e^ij; = ej° 


(6.78) 


which has the solution 


A 0 (x, t) — e d 3 x 




A-k\x — x’\ 

In Coulomb gauge we can rewrite the Maxwell part of the Lagrangian as 


(6.79) 


^Maxwell — 


d 3 x \E 2 - \B 2 


= J d 3 x l(A + VA 0 ) 2 - ±B 2 
= [ d 3 x \A 2 + i(V. A 0 ) 2 - \B 2 


(6.80) 


where the cross-term has vanished using V ■ A = 0. After integrating the second term 
by parts and inserting the equation for A 0 , we have 


^Maxwell 


= / d 3 x 


1 1 D 2 




/ jo(x)jo(x') 


±B Z + - d 6 x' A _ . . 

2 2 J 4tt\x — x'\ 


(6.81) 


We find ourselves with a nonlocal term in the action. This is exactly the type of 
interaction that we boasted in Section 1.1.4 never arises in Nature! It appears here as 
an artifact of working in Coulomb gauge: it does not mean that the theory of QED is 
nonlocal. For example, it wouldn’t appear if we worked in Lorentz gauge. 

We now compute the Hamiltonian. Changing notation slightly from previous chap¬ 
ters, we have the conjugate momenta, 

n = A = i 


dA 

dC -/t 

** = e$ = '* 


(6.82) 


which gives us the Hamiltonian 


H = d 3 x 


\A + ^B + + m)ip - ej ■ A + — / d x‘ 


47t|x — x’\ 


where j = and j° = 


140 









6.4.1 Naive Feynman Rules 

We want to determine the Feynman rules for this theory. For fermions, the rules are 
the same as those given in Section 5. The new pieces are: 

• We denote the photon by a wavy line. Each end of the line comes with an i,j = 
1, 2, 3 index telling us the component of A. We calculated the transverse photon 

propagator in (6.33): it is v/wwv and contributes Df- = —- ( <5 i? - — ( ) 

J p 2 + ie\ \p\ 2 ) 

• The vertex contributes — ie r f. The index on 7 * contracts with the 

index on the photon line. 


• The non-local interaction which, in position space, is given by 

contributes a factor of ~ ^ 

4:7t\x — y | 

These Feynman rules are rather messy. This is the price we’ve paid for working in 
Coulomb gauge. We’ll now show that we can massage these expressions into something 
much more simple and Lorentz invariant. Let’s start with the offending instantaneous 
interaction. Since it comes from the Aq component of the gauge field, we could try to 
redefine the propagator to include a D 00 piece which will capture this term. In fact, it 
fits quite nicely in this form: if we look in momentum space, we have 



5(x° — y°) f d A p e ip '( x yS> 
An\x-y\ J (27t) 4 \p\ 2 


(6.83) 


so we can combine the non-local interaction with the transverse photon propagator by 
defining a new photon propagator 


D^{p) 


p 2 + ie 


PiPj 


IpT 

s,i ~ IpT 

0 


H,u — 0 
/i = i^0, z/ = j^0 
otherwise 


(6.84) 


With this propagator, the wavy photon line now carries a /x, v = 0,1, 2,3 index, with 
the extra /1 = 0 component taking care of the instantaneous interaction. We now need 
to change our vertex slightly: the —ie 7 * above gets replaced by —ie 7 ^ which correctly 
accounts for the (e 7°) 2 piece in the instantaneous interaction. 
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The D 00 piece of the propagator doesn’t look a whole lot different from the transverse 
photon propagator. But wouldn’t it be nice if they were both part of something more 
symmetric! In fact, they are. We have the following: 

Claim: We can replace the propagator D^ v (p) with the simpler, Lorentz invariant 
propagator 


DM = (6.85) 

pZ 

Proof: There is a general proof using current conservation. Here we’ll be more pedes¬ 
trian and show that we can do this for certain Feynman diagrams. In particular, we 
focus on a particular tree-level diagram that contributes to e~e~ —>■ e~e~ scattering, 



e 2 [u(p')Yu(p)\ D^ik) [u(q)Yu(q)] 


( 6 . 86 ) 


where k = p — p' = q' — q. Recall that u(p) satisfies the equation 

(rf — m)u(p) = 0 (6.87) 

Let’s define the spinor contractions a M = u(p')'y fJ ’u(p) and f3 v = u(q')^ v u(q). Then 
since k = p — p' = q' — q, we have 


k^a^ = u{p')(]f — /)u(p) = u(p')(m — m)u(p) = 0 


( 6 . 88 ) 


and, similarly, k u /3 u = 0. Using this fact, the diagram can be written as 


1 a ■ /3 (a ■ k)(f3 ■ k) + a°/3° 


ofDwP’ = 


= l 


= l 


k 2 k 2 \k\ 2 \k\ 2 

a-/3 kla 0 f3 0 a 0 /3 0 


k 2 

a ■ /3 




k 2 \k\ 2 
1 

k 2 \k\ 2 

= n*( "k 


+ \k\ 2 

(k 2 0 - k 2 ) a 0 /3 0 

P 


k 2 


(6.89) 
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which is the claimed result. You can similarly check that the same substitution is legal 
in the diagram 

~ e 2 [v(^Yu(p)}D^(k)[u(p')Yv(q')] (6.90) 

In fact, although we won’t show it here, it’s a general fact that in every Feynman dia¬ 
gram we may use the very nice, Lorentz invariant propagator = —irj^/p 2 . □ 



Note: This is the propagator we found when quantizing in Lorentz gauge (using the 
Feynman gauge parameter). In general, quantizing the Lagrangian (6.37) in Lorentz 
gauge, we have the propagator 

D n» = + (« - !)^) (6-91) 

Using similar arguments to those given above, you can show that the p^p^/p 2 term 
cancels in all diagrams. For example, in the following diagrams the p^p v piece of the 
propagator contributes as 



~ u(p')Yu{p) k = u(p')(tf- / )u(p) = 0 

>-~ v(p)^u(q) kn = u(p)(rf + i')u{q) = 0 

6.5 Feynman Rules 

Finally, we have the Feynman rules for QED. For vertices and internal lines, we write 


(6.92) 



Photon Propagator: v/X/X/X/X/V 




Fermion Propagator: 


> 


p 2 + ie 
i( 4 + m) 


p 2 — m 2 + ie 
For external lines in the diagram, we attach 

• Photons: We add a polarization vector e^ n /e^ ut for incoming/outgoing photons. 
In Coulomb gauge, e° = 0 and e-p — Q. 

• Fermions: We add a spinor u r (p)/u r (p) for incoming/out going fermions. We add 
a spinor v r (p)/v r (p) for incoming/outgoing anti-fermions. 


143 






6.5.1 Charged Scalars 

“Pauli asked me to calculate the cross section for pair creation of scalar 
particles by photons. It was only a short time after Bethe and Heitler 
had solved the same problem for electrons and positrons. I met Bethe 
in Copenhagen at a conference and asked him to tell me how he did the 
calculations. I also inquired how long it would take to perform this task; 
he answered, “It would take me three days, but you will need about three 
weeks.” He was right, as usual; furthermore, the published cross sections 
were wrong by a factor of four.” 

Viki Weisskopf 


The interaction terms in the Lagrangian for charged scalars come from the covariant 
derivative terms, 

C = V^ Vd^ip — ieA^^d^ip — _|_ e 2 A^A^ip^tp (6.93) 

This gives rise to two interaction vertices. But the cubic vertex is something we haven’t 
seen before: it contains kinetic terms. How do these appear in the Feynman rules? 
After a Fourier transform, the derivative term means that the interaction is stronger 
for fermions with higher momentum, so we include a momentum factor in the Feynman 
rule. There is also a second, “seagull” graph. The two Feynman rules are 



The factor of two in the seagull diagram arises because of the two identical particles 
appearing in the vertex. (It’s the same reason that the 1/4! didn’t appear in the 
Feynman rules for 0 4 theory). 


6.6 Scattering in QED 

Let’s now calculate some amplitudes for various processes in quantum electrodynamics, 
with a photon coupled to a single fermion. We will consider the analogous set of 
processes that we saw in Section 3 and Section 5. We have 
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Electron Scattering 


Electron scattering e e —>■ e e is described by the two leading order Feynman dia¬ 
grams, given by 

2 f {u s \p')Yu s {p)} [u r \q'h^u r (g)\ 

V (p 1 - p ) 2 

_ [u s \f)Yu r {q)\ [u f '(f)7/(p)] \ 
(p - q ') 2 ) 

The overall — % comes from the —iq^ in the propagator, which contract the indices on 
the 7 -matrices (remember that it’s really positive for / 1 , v = 1,2, 3). 



Electron Positron Annihilation 


Let’s now look at e e + —> 2y, two gamma rays. The two lowest order Feynman 
diagrams are, 



i(—ie) 2 v\q) 


f 7m(^~ rf' + m) 7 ^ 

\ (p — p') 2 — m 2 


l' + rn) 7/A 

(p — q’) 2 — m 2 J 


u s {m{p’K{q') 


Electron Positron Scattering 

For e _ e + —> e~e + scattering (sometimes known as Bhabha scattering) the two lowest 
order Feynman diagrams are 




u s \p')Yu s (p)\ [v r (q)inv r '(q')_ 

(p — p') 2 

[v r (q)yu s (p)} [u s '(p') r r fJ ,v r '(q')} ’ 

(p + q ) 2 


Compton Scattering 

The scattering of photons (in particular x-rays) off electrons e _ 7 —> e _ 7 is known 
as Compton scattering. Historically, the change in wavelength of the photon in the 
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scattering process was one of the conclusive pieces of evidence that light could behave 
as a particle. The amplitude is given by, 



This amplitude vanishes for longitudinal photons. For example, suppose e m ~ q. Then, 
using momentum conservation p + q = p' + </, we may write the amplitude as 


iA = i(—ie) 2 u r '(p') 


= i(-ie)V'(p') 4 mit u s (p) 


/out(^+ 4 + m) 4 4( / - 4 + m) 

(p + q ) 2 — m 2 (p' — q ) 2 — m 2 

( 2 pa t 2 p'-q 

(p + q ) 2 — m 2 (p' — q ) 2 — m 2 


u s (p) 


(6.94) 


where, in going to the second line, we’ve performed some 7 -matrix manipulations, 
and invoked the fact that q is null, so 4 4 = 0, together with the spinor equations 
{rf-m)u{p) and u(p')(rf' — m) = 0. We now recall the fact that q is a null vector, while 
p 2 = (p ') 2 = m 2 since the external legs are on mass-shell. This means that the two 
denominators in the amplitude read (p+g) 2 — m 2 = 2 p-q and (p'—q)—m 2 = —2 p'-q. This 
ensures that the combined amplitude vanishes for longitudinal photons as promised. A 
similar result holds when / out ~ q’. 


Photon Scattering 

In QED, photons no longer pass through each other unimpeded, y 

At one-loop, there is a diagram which leads to photon scattering. b—<—T 

Although naively logarithmically divergent, the diagram is actually x/ /x 
rendered finite by gauge invariance. ^ LLla H 

Adding Muons Figure 30: 

Adding a second fermion into the mix, which we could identify as a 
muon, new processes become possible. For example, we can now have processes such 
as e~p~ —> e~ p~ scattering, and e + e~ annihilation into a muon anti-muon pair. Using 
our standard notation of p and q for incoming momenta, and p' and q 1 for outgoing 
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momenta, we have the amplitudes given by 



and 



6.6.1 The Coulomb Potential 

We’ve come a long way. We’ve understood how to compute quantum amplitudes in 
a large array of field theories. To end this course, we use our newfound knowledge to 
rederive a result you learnt in kindergarten: Coulomb’s law. 

To do this, we repeat our calculation that led us to the Yukawa force in Sections 

3.5.2 and 5.7.2. We start by looking at e _ e _ —> e~e~ scattering. We have 



v ^ [u{p')Y u (p)\ W(q'hMv)} 

H le ) „\2 


(6.96) 



means that the 7 0 piece of the interaction gives a term u s (p)Yu r (q) —> 2 mS rs , while 
the spatial 7®, i = 1,2,3 pieces vanish in the non-relativistic limit: u s (p)^ l u r (q) —>■ 0. 
Comparing the scattering amplitude in this limit to that of non-relativistic quantum 
mechanics, we have the effective potential between two electrons given by, 



(6.97) 


We find the familiar repulsive Coulomb potential. We can trace the minus sign that 
gives a repulsive potential to the fact that only the Aq component of the intermediate 
propagator ~ —if]^ contributes in the non-relativistic limit. 

For e _ e + —>■ e _ e + scattering, the amplitude is 



+i (_ ie ) 2 Hp')Yu(p)} [v(q)iMq')} 


( 6 . 98 ) 









The overall + sign comes from treating the fermions correctly: we saw the same minus 
sign when studying scattering in Yukawa theory. The difference now comes from looking 
at the non-relativistic limit. We have v^/°v —y 2m, giving us the potential between 
opposite charges, 


U(r) = — e 2 


d 3 p 

(27r) 3 \p \ 2 


47rr 


(6.99) 


Reassuringly, we find an attractive force between an electron and positron. The differ¬ 
ence from the calculation of the Yukawa force comes again from the zeroth component 
of the gauge field, this time in the guise of the 7 0 sandwiched between v^°v —» 2m, 
rather than the vv —y —2m that we saw in the Yukawa case. 


The Coulomb Potential for Scalars 

There are many minus signs in the above calculation which somewhat obscure the 
crucial one which gives rise to the repulsive force. A careful study reveals the offending 
sign to be that which sits in front of the A 0 piece of the photon propagator —irj llL ,/p 2 . 

Note that with our signature (H-), the propagating A* have the correct sign, while 

A 0 comes with the wrong sign. This is simpler to see in the case of scalar QED, where 
we don’t have to worry about the gamma matrices. From the Feynman rules of Section 
6.5.1, we have the non-relativistic limit of scalar e _ e _ scattering, 



where the non-relativistic limit in the numerator involves (p+p') ■ {q+q') ~ ( p+p')°(q + 
q ')0 « ( 2 m ) 2 and is responsible for selecting the A 0 part of the photon propagator rather 
than the A* piece. This shows that the Coulomb potential for spin 0 particles of the 
same charge is again repulsive, just as it is for fermions. For e _ e + scattering, the 
amplitude picks up an extra minus sign because the arrows on the legs of the Feynman 
rules in Section 6.5.1 are correlated with the momentum arrows. Flipping the arrows 
on one pair of legs in the amplitude introduces the relevant minus sign to ensure that 
the non-relativistic potential between e _ e + is attractive as expected. 
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6.7 Afterword 

In this course, we have laid the foundational framework for quantum field theory. Most 
of the developments that we’ve seen were already in place by the middle of the 1930s, 
pioneered by people such as Jordan, Dirac, Heisenberg, Pauli and Weisskopf 5 . 

Yet by the end of the 1930s, physicists were ready to give up on quantum held theory. 
The difficulty lies in the next terms in perturbation theory. These are the terms that 
correspond to Feynamn diagrams with loops in them, which we have scrupulously 
avoided computing in this course. The reason we’ve avoided them is because they 
typically give infinity! And, after ten years of trying, and failing, to make sense of this, 
the general feeling was that one should do something else. This from Dirac in 1937, 

Because of its extreme complexity, most physicists will be glad to see the 
end of QED 

But the leading figures of the day gave up too soon. It took a new generation of postwar 
physicists — people like Schwinger, Feynman, Tomonaga and Dyson — to return to 
quantum held theory and tame the infinities. The story of how they did that will be 
told in next term’s course. 


5 For more details on the history of quantum field theory, see the excellent book “QED and the 
Men who Made it” by Sam Schweber. 
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